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NOTICE TO AUTHORS 


1. Communications.—Papers must be communicated to the Society by a Fellow. They 
should be accompanied by a summary at the beginning of the paper conveying briefly the 
content of the paper, and drawing attention to important new information and to the main 
conclusions. ‘The summary should be intelligible in itself, without reference to the paper, 
to a reader with some knowledge of the subject; it should not normally exceed 200 words 
in length. Authors are requested to submit MSS. in duplicate. These should be 
be ae Pp me sec apna margin of not less than one inch on the 

hand side. Corrections to the MSS. should be made in the text and not in the 

margin. Unless a paper reaches the Secretaries more than seven days before a Council 

poem it will not normally be considered at that meeting. By Council decision, MSS. of 

accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation.—Authors are allowed considerable latitude, but they are requested to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespective of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively :— 

A. Corlin, Zs. f. Astrophys., 15, 239, 1938. 
H. Jeffreys, Theory of Probability, and edn., section 5.45, p. 258, Oxford, 1948. 


3. Notation.—For technical astronomical terms, authors should conform closely to the 
recommendations of Commission 3 of the International Astronomical Union (Trans. 
I.A.U., Vol. VI, p. 345, 1938). Council has decided to adopt the I.A.U. 3-letter abbrevi- 
ations for constellations where contraction is desirable (Vol. IV, ms 221, 1932). In general 
matters, authors should follow the recommendations in Symbols, Si igns and Abbreviations 
(London: Royal Society, 1951) except where these conflict with I. AU. U. practice. 


4. Diagrams.—These should be designed to appear upright on the page, drawn 
about twice the size required in print and prepared for direct photographic 
reproduction except for the lettering, which should be inserted in pencil. 
Legends should be given in the manuscript indicating where in the text the 
figure should appear. Blocks are retained by the Society for 10 years; unless the author 
requires them before the end of this period they are then destroyed. 


5. Tables—These should be arranged so that they can be printed upright on 
the page. 


6. Proofs.—Costs of alterations exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society’s interests to seek the maximum conciseness and simplifi- 
cation of symbols and equations consistent with clarity. 


7. Revised Manuscripts —When ee are submitted in revised form it is especially 


requested that they be accompanied by the original MSS. 


Reading of Papers at Meetings 


8. When submitting papers authors are requested to indicate whether they will be 
willing and able to read the paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 


g. Postcards giving the programme of each meeting are issued some days before the 
ing concerned. Fellows wishing to receive such cards whether for i 


Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 
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MEETING OF 1954 JANUARY 8 
Dr J. Jackson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


*Lewis Grenville Bentley, 27 Gracefield Gardens, Streatham, London, S.W.16 
(proposed by F. M. Holborn); 

Robert John Bray, University Observatory, Oxford (proposed by H. H. 
Plaskett) ; 

Ronald Edgar Edwards, 14 Leaside Way, Southampton (proposed by E. W. 
Mason) ; 

*Mavis McIntyre Smith Gibson, H.M. Nautical Almanac Office, Royal 
Greenwich Observatory, Herstmonceux Castle, Hailsham, Sussex (pro- 
posed by D. H. Sadler) ; 

*Raymond Hide, 22 Castle Hills Road, Scawthorpe, near Doncaster, Yorkshire 
(proposed by S. Chandrasekhar) ; 

*John Noel Hodgson, 123 Cherrydown Avenue, Chingford, London, E.4 
(proposed by R. O. Redman); 

*Francis John Lowes, Department of Geodesy and Geophysics, Downing 
Place, Cambridge (proposed by M. N. Hill); 

William Scott Matthews, 81 Albert Road, Colne, Lancashire (proposed by 
W. M. Mason) ; 

*Gilbert Elliott Satterthwaite, 48 Upper Abbey Road, Belvedere, Kent 
(proposed by A. F. Alexander) ; 

*Gordon Ingram Thompson, Dunsink Observatory, Co. Dublin, Eire (proposed 
by H. A. Briick); and 

William Henry Vale, 11 Edgar Street, Brighton, Victoria, Australia (proposed 
by J. Boyle). 
The election by the Council of the following Junior Member was duly 
confirmed :— 
Philip Henry Vince, Wyeross, 17 Stanley Hill Avenue, Amersham, Bucking- 
hamshire (proposed by A. W. Vince). 


Ninety-seven presents were announced as having been received since the last 
meeting. 

The President announced that the Council had awarded the Gold Medal of the 
Society to Dr Walter Baade, of Mount Wilson Observatory, for his observational 
work on galactic and extra-galactic objects. 

* Transferred from Junior Membership. 





ON THE INTENSITIES OF INTERLOCKED MULTIPLET 
LINES IN THE MILNE-EDDINGTON MODEL 


I. W. Busbridge and D. W. N. Stibbs 
(Communicated by the Director, University Observatory, Oxford) 


(Received 1953 July 11) 


Summary 


The equations of transfer in the Milne—Eddington model for multiplet 
lines with a common upper state are solved exactly. It is assumed that no 
redistribution in frequency takes place other than that due to interlocking, 
that the ratios of the line absorption coefficients to the continuous absorption 
coefficient are independent of depth, and that the lower states of the lines are 
sharp. A linear approximation for the Planck function is used. The solution is 
found to depend on H-functions of a fairly simple type. The case of a 
doublet is considered in detail. Values for the H-functions are computed 
and the residual intensities of the two components are evaluated as functions of 
»(=cos 8) for three distances from the centres of the lines. These are 
compared with the residual intensities for doublet lines when there is no 
interlocking. It is shown that interlocking causes a small reduction in the 
difference between the residual intensities of the two lines, the reduction 
being greater towards the limb. A solution for the doublet intensities based 
on Eddington’s approximation is also considered and it is shown that the 
errors in the approximate solution are satisfactorily small for 4 >0-2. 





1. Introduction.—The theory of the formation of absorption lines by coherent 
scattering may be applied to multiplets with a common lower state if the lower 
state itself is sharp and if there is no redistribution in frequency due to collision 
broadening of the upper states. Under these conditions an atom, once excited 
to any one of the upper states, can only return to the common lower state. When, 
however, there is a common upper state, the re-emission will be governed by the 
transition probabilities corresponding to the various lines and the atom may 
return to a different state from the one in which it started. Examples of such 
interlocked lines are 


, 
A3961 A rz 2 P59 — 2S. 
I 3991 ad Gal 4172A lel . 1/2 
39444 4032A *Pyj2—* Si 
A5184A A6162A 8P, —3S, 
MgI< 5173A and Cal; 6122A °P, —%S, 
| 5167A | 6102A 3Py —8S, 

The equations of transfer in the Milne-Eddington model for interlocking 
without redistribution have been discussed by Woolley and Stibbs in (5), where 
a clear statement of the problem will be found. Taking the Planck function to 
be linear, they have obtained a solution by means of Eddington’s approximation 
and calculated the residual intensities and the total absorption in the emergent 


flux for doublet and triplet lines. Although their equations are set out for 
triplets, they clearly extend to multiplets containing any number of lines with a 
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common upper state and these equations, suitably generalized, are reproduced 
in Section 2. 

The method of solution (Sections 3-5) is a slight modification of one which 
has already been applied successfully in (x) and (4) to problems of non-coherent 
scattering. It easily gives the exact solution for k interlocked lines in terms of 
H-functions of a fairly simple type. 

When k=1, these functions are those already tabulated by Chandrasekhar 
and Breen (see (2) Table XI). In Section 6.1 we describe the method which 
has been used to compute them when k= 2 and we give a table of values (Table I). 
These values are used in Section 6.2 to evaluate the residual intensities for 
doublet lines at three distances from their centres. The residual intensities are 
calculated as functions of 4, where cos! is the angle which the direction of 
emergence makes with the outward normal to the surface of the atmosphere. 
The results are shown in Table II and Fig. 1. 

In order to show the effect of interlocking on the centre-to-limb variation of 
the residual intensity, the corresponding residual intensities are evaluated for a 
doublet without interlocking (‘Table III and Fig. 2). It is shown that inter- 
locking reduces the differences between the residual intensities of the two lines 
for all values of 1. The reduction near the centre of the disk (4 = 1) is very small, 
but it increases towards the limb. The effect is shown clearly in Fig. 3. 

The labour involved in computing the H-functions is considerable and it 
would be even greater for a triplet (k=3). For this reason we consider, in 
Section 6.3, the accuracy of an approximate solution of the equations of transfer 
for a doublet based on Eddington’s approximation, and we show that, in the 
cases examined, the proportional errors in the residual intensities lie within 
+1 per cent for values of 4 >0-2. 

2. The equations of transfer.—Woolley and Stibbs have made certain assump- 
tions in (§) which we shall also make here. These are 

(i) that the lines are so close together that variations of the continuous absorp- 
tion coefficient « and of the Planck function B(v, T) with wave-length may be 
neglected. This also means that the lower states are nearly equal in excitation 
potential and that they have the same classical damping constant. Then the 
values of y;, . .. , »,(the ratios of the line absorption coefficients tox) are proportional 
to the transition probabilities for spontaneous emission from the upper state to 
the respective lower states ; 


(ii) that y,,...,7, are independent of depth; 
(iii) that Biv, T)=a+br, (2.1) 


where 7 = | xp dx, x being the depth below the surface of the atmosphere: By (i) 
0 


a, 6 and 7 are independent of v; 
(iv) that the coefficient e, which is introduced to allow for thermal emission 
associated with the line absorption, is independent of both frequency and depth. 
With these assumptions, the equations of transfer are 


ote) =(1+7,)(7,4)— (+e, a +br) 


k 1 
—(t-)a,24np{ Iy(rse')du! (r= 1,009) 
p=1 -1 


where t= -/(M1 +++ +x); 
so that Ay +agt... +a, =I. 
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In the equations (2.2) the subscript r denotes the quantity corresponding to the 
line of frequency v,. Later it will be convenient to use J,(7, 4) and », when we 
are considering any frequency v, but on putting y= v,, we shall replace them by 
T(z, ») and »,. 

The equations (2.2) have to be solved subject to the boundary conditions 

I(0,-p')=0  (0<p' <I, r=I,...,h), (2.5) 

together with a condition limiting J,(7,u) for large +. We shall assume that 
I,{7, 4) is at most linear in 7 for large r. 

Exact solutions of (2.2) are easily found, but they do not satisfy (2.5). These 


Irs) =0+ (7+ +) (r 


are 


as is easily verified with the help of (2.3). 
Write 
) +1,*(7,u) (r=1,...,). (2.6) 


Then /,*(7, «) satisfies the equations of transfer 


*(7 a “ 
dl, “et =(1+7,)I,*(7, u)—(1—«€)a, = bnp| I, *(7, p’)dp’ (r= ae 
p=1 i 2.7 


Ifo,n)=a+b (r+ 


I+” 


r 


together with the boundary conditions 


I,*(0, -p) =e —a@ (0<p’ <I, r=I,...,h). (2.8) 


Moreover J,*(7,) must be at most linear in + as tr 00. 

We have here the problem of a scattering atmosphere subject to external 
radiation, whose intensity is given by (2.8). We want to find the emergent 
intensity J,*(0,) of frequency v,. This will be the intensity of the diffusely 
reflected radiation; it can be calculated when the appropriate scattering function 
is known. 

In the present problem the scattering function splits up into k? functions 


Sfp’) (r=t,...,k; s=I,...,h), 
but it is convenient to reunite them temporarily in the function 
Pl», v’)S(v, v’ > Bs Be’), 
where vr is any one of 14, v9,..., Vp» 
. . 
Hr»’)=a, Z3(vy—v') (2.9) 
6 denoting Dirac’s 5-function, and 
S(v,, Vs> By B’')=S,(e, Bw’). (2.10) 
Then the law of diffuse reflection for the atmosphere can be written (cf. (4), 
equation (15)) 
ra 1 
1(0,n)= > | ol v')dv’ | S(v,5 wou’ (0, —n')du’- (2.22) 
* 0 0 


, 


The equivalent form in terms of the functions S,,(u, 4’) is 


' I fa i ’ , 
1,™(0,n)=2,3 = f° Sip(as’Up™(0, —H') de’ (2.12) 
p=1-h lg: 


a form to which we shall return in Section 5. 
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With the aid of (2.9), (2.7) can be written in the form 
BOC) tn )l.e(rs0)—ME—6) | nePlse'dav’ [Leesa !)de (2.23) 


This is equation (12) of (4), to which we now refer for the general equation for 


the scattering function. 
3. The scattering functions.—In (4), Stibbs has found the following general 


equation for the scattering function S(v, v’; ju, 4’): 


: 
[+ EE] EM — 6, 

ae 
| 2h" 
ia al 


Sv’, ¥5 BS Ma) a 





1 
+n | P(», v’)p(v’, v;) dv’ Fh cide 5 pts 


+ 2 nvP(», v’)p(v’, v;) dv’ 


+ |p P(r, v") dv” 4 S(v, v"; ',v;) dv’ 


, 


du 
S(v',%5 by Mas: I 
xf (V4 HB Mi). (3-1) 


In this v;and y, refer to the incident radiation and v, » to the emergent radiation. 
On substituting the value (2.9) for p(v, v’) and dividing by «,, this equation 
becomes 
I+yn, I+ne]e Sv, v5 by Hy 
[ : , = a. y,) Slits Hs bi) =7,,25(v,—%) 
be i 8 I—e€ 8 


i ‘ du 
+1, 28(%—m) | | 2vEB(Up— v') dy’ I S(v, v's use") 5 
8 ¢ p 





0 d , 
+ X8(v,—v,) | %yLO(V_ — v’) dv’ 2 S(v’, v5 o's Mi) 4 
8 J0 qa 


+23(¥,-%) j | yEB( vy —v") do" f S(v, v”; me 
F ;, — 
x f ay v2 P(g —v')dv’ f S(v', 33 Pe 


On using the properties “ the 5-function we get 


I+yM, Iti le S(v, 455 Hy Mi) 
[ “ae i |=. we eee 





. wa » dp’ 
+,,28(v,— 44) Dat, J SQ %5 HET 
3 p 


> . dw 
+ 28(v,— V;)U%—Ng J ‘ink 43 By 2 
8 q bh 


+33(v,— nay | Suisi wy 
8 Pp 0 lied 


. ¥ » \ dp 
x LaeNg | S(vq Yj3 B's bi) 2u" 
q 0 lod 


When v,;=»v, the 5-function becomes indefinitely large; hence this equation will 
be satisfied only if the coefficient of 5(v,— 1) is zero for y,=v,. Hence 
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|; +n, is I+n, | S(v, v3 Ms Ps) 


B Vy I-—e 





iy 1 4h 

=n n2ay | S(v, Vy> Bs an 

du du’ 

+Za Sv, V3 H+ Bx, | S(r.¥93 Wt ee 

Sata |. (Yq Mss HoH) o (%%3 HoH) D 
x Dagte | S( Vg Vs5 heats 
q 0 2p 


Finally replacing v by v,, using (2.10), and factorizing, we have 





I+, . 1+, | S,(Hs 1) 2 1 » dp’ 
[ + ] a eee fa. |, Silo eo 


du’ 
x {n+ 3 = tata. Soe’ H;) ae 


lad Bj 


and this holds for r=1,...,k, s=1,...,k. 
From the theory of coherent scattering we know that 


S,(u, My) = S,(Mis b). (3-3) 


On putting s=r in (3.2), we see that this will be so if 


pS plt’s H)= 1S rp(Hs #’), (3-4) 


a relation which includes (3.3). In the Appendix we prove. conversely that 
(3-2) and (3.3) together imply (3.4). From (3.2) and (3.4) we have 


I+%, I+, Sy o(My Hi) = {r+ . ’ al 
[ oo al wi- z ‘ |, Slee) 2 


14+ 5 ['s, “nth. ' 
| = 0 oh Ou (3-5) 





Write 
A,(u)=1+ = Xp {. Srp(Hs oH = 2m? (3.6) 
so that 





|; +, + I =1e] S(Hs Hi) 


be by; nA(I— €) =H Au)A, (144), (3-7) 


and let* 


I +e, I 


tie coven - = . 8 
A, ite.’ n, +7, (3-8) 





Then 


nm, I+M%s 
cre (3-9) 


1~A, ai den J 
I+”, 
and (3.7) can be written 


Sre(Hs B’) =n{I— A eon A (p)H (p' ). (3-10) 


a 





* This definition of Ay agrees with that adopted in (2), Section 84. If e=o, then mr=Ar. 
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On substituting from (3.10) into (3.6), we have 


H{u)=1+4 E ay(1—2, null (u) | AD gy, 
; ont ‘ 0M, + Nye 


Ae’) 4 


ree wide 


1 k 
a 3 RA %p(I—A,)n,p joes 
Consider also 


—4 3 ag(—Aya fe ay 


Be 
Hu) 

From these it at once follows that 
H,(x/n,)=H,(x/n,)=H(x) (say). 


H,(u)=H(n,p), (3-13) 


Hence 


where, from (3.11), 


Hn) =X + In uHng) © %,(I—A,) [ie ' 


; I 
ot +n (3.14) 
Change the variables in this from yu and y’ to x and x’ by writing 2,u =x, np’ =x’. 
Then 

"» H(x') 


0 x+x' 


H(x)=1+ $xH(x) z %,(I1—A,)n,* | ; dx", 
1 F(x')H(2’) 


H(x)=1+xH(x) J. os 
where, assuming that 


> Ne>+++> Nk» 
so that 


oO ey ++ <M, SI, 
%,(I—A,)/n, (0 <x’ <n), 


¥(x')= 4 p(I—Ap)/My (n, <x" <m,41), 


(n, <x’ <1). 








{. F(x") dx’ = = $ay(I--A,) = (x- > a 4 >») <, (3-19) 


‘Y'(x’) satisfies the conditions of (2), Chap. V, and the solution of (3.15) can at 
once be found. 
From (3.10) and (3.13) we have 


S,s(H5 #’) =7,(1 —A,.) ———, H(n,»)H(n’). (3.20) 


<n 


This is the final form for the scattering function from frequency v, and direction 
—p’ into frequency v, and direction p. 
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4. The H-function.—Before completing the solution of the problem in question, 
it is convenient to obtain some properties of the H-function which is the solution 
of (3.15). Inthis section all references are to Chap. V of (2) and we shall therefore 
quote (2) followed by the number of the relevant formula or by the theorem 
required. 

By (2), Theorem 1, and by (3.19), 


k “n j 1/2 
D 4a,(1—A,)m,* | ” H(x’) dx’ =1—- ( >> ayy) 
p=1 0 p=1 


k rl & 1/2 
D> dop( Ay) | H(n,p’) du’ =1— (= x ay) (4.1) 
p=1 0 p=1 


By means of (4.1), (3.14) can be re-written in the form (cf. (2), Theorem 1, 
Corollary 2) 
oe np H(n uw’), , 
H(n,,1) = (3 2a ty) +43 a E ay ( 1 _ ry) { eT a (4.2) 


0 2+ Nye 


By (2), (66) and (67), 


sedygl In T(w) a 


where 





1P(u) du = ) (" 
T(w)=1-2 a) = %p vB ("x 
(») 0 w*— 3 ve oe. 0 ss 
This gives 
A,),. win 


k om 
Te)=1-je © I) 8 t% | (4.4) 
p=1 = 


'p w-—n, 





where the principal values of the logarithms have to be taken. 

It is simplest now to concentrate on the strongest line of the multiplet and to 
calculate H(n,). ‘The changes required for any other line are obvious. In 
(4.3) and (4.4) put x=n,u, w=in, cot@ (0<@<47). Then (4.4) simplifies to 


k 
T(in, cot @) =1 cot {ay( ~),)0+ ¥ a,(1—-A,)— cot (> cot 8) b, (4.5) 
p=2 Ny Ny 


and (4.3) becomes 
How) e0{—" 


7 





(4-6) 


When »,..-,, and ¢ are known, «,, A, and n, (r=1,...,) can be calculated 
by means of (2.3) and (3.8) and then H(n,u) can be evaluated from (4.5) and (4.6) 
by numerical integration. Had the strongest line of the multiplet been formed by 
coherent scattering (without interlocking) the H-function involved would have 
been given by (4.6) with 


T(in, cot @) = 1 —(1 —A,)0 cot 0. (4.7) 


By comparing this with (4.5) we see the changes introduced by interlocking. 
5. The emergent eae (2.8), (2.12) and (3.8), we have 


a2 In T(in, cot 6) 
9 cos?6 +p? sin? 6 ? 


1,*(0,u)= as f Syp(vts 2’ \(bmyp’ — a) du 
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Substituting from (3.20) we get 


k rl , b en ) 
* - . me B ( Ny —a D ’ 
1,*(0, 4) = }%,H(n,) & n,(t-2y) ag Se rma tot’) dy 


k 1 (bn,u+a 
eo > = , i: Peete 
= $a,H(n,) KA n,{I—A,) i; {ou Nf + Ny’ 


= }o,H(n,1) a E (2—Ap) |’) de 
wH(nye’) ‘. 


, 
0 2, + Npyft 


+ Hyp )dp’ 


—nfa+bngu) © (x-A my) [ 
From (2.6), 


(5-1) 


I,(0, 4) =a + bn,p +I1,*(0, 1). (5-2) 
By making use of (4.2) multiplied by H(n,), this can be written in the form 


1,(0,u)=(a+ bmg) H(nu){ (> isd +4 E oy(r—a,) [CoD ay'h 


0 nb + Ny’ 
1 w'H (nyu') 


tte rere ea ree 
= r Pp 


rl 


+ 4ba.n, Hn) © E (1 ~h) p(y’) du’, 
and thus 


1,(0, «) = (a+ bn,u)H(n,1) { (= ) typ), 2 


ad “1 "H (n ’) 
s = a 6 es et 
* p=1 (pty — ate —Ay) oN, + Nye . } 


k rl 
+ Jba,H (nu) & (t—Ay) | mHtpe') du’. (5.3) 


This is the final form for the emergent intensity in the rth line. 

We shall use the notation H(w; ») for the H-function tabulated by 
Chandrasekhar and Breen (see (2), Table XI). Then H(1—A; yz) is the 
H-function for coherent scattering. When k=1, (4.5) reduces to (4.7), 
A(ny) = H(1—A,; ») and (5.3) becomes Chandrasekhar’s solution for coherent 
scattering ((2), Section 84). This solution is alsoobtained when yn, = y2.=...=7,- 


6. The calculations for a doublet 
6.1. The H-functions.—For a doublet, ; 
°42=2: I and a= %, Oe = 4. (6.1) 
Let ¢(0 <¢ < 4) be defined in terms of 6(0 <8 < 47) by the equation 
n, cot? =n, cot ¢. (6.2) 
Then, when k=2, (4.5) and (4.6) give 
H(mps)=exp {= er In[1 —a,(1—A,) @ cot é—a,(1—A,) ¢ cot d| ao} 


0 cos? 6 +? sin? @ 
and we have the relation 





» (6.3) 


H(my:)=H(m."n). (6.4) 
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The integral (6.3) could be evaluated by numerical integration, but it was 
decided to calculate the values of the H-functions by means of sequence 
iterations based on (4.2). For this, good initial approximations for H(n,) 
and H(n.) were needed and these were obtained from (6.2)-(6.4) as follows: 

Since «a, +a,=1, (6.3) can be written 


ri | —[1I—«,A, —%A.—k(8)] Ocoté 
Hm) exp} — * {" nee one rns. ah, (6.5) 





Tt 
where 


h(0) =09(1 —Aq){1 — (d cot $)/(8 cot 4)}. (6.6) 


As @ increases from 0 to $7, k(@) increases from 0 to «(I —A,)(I — ,/mg), which is 
less than 4 in all cases. As a first approximation k(@) can be replaced by 


hy = }a9(1 —Ay)(I —m/Mg), (6.7) 


H( nyu) = H(1— ayy — ayy 1): (6.8) 


and then (6.5) gives 


Thus approximate values for H(n,u), when 0<p <1, can be found from the 
tables for H(w; ,) in (2), interpolating linearly between the values for successive 
@ when necessary. 

In order to evaluate H(n.) from (6.4), it is necessary to know values for 
H(n,) not only in the interval (0, 1) but also for I<p<a,/n,. Since, by (4.2), 


Hm) (% A, +%Ag) 4? as p> 00, (6.9) 
fairly reliable values can be obtained by an extrapolation formula of the form 
H( my) = (4 Ay + %_Aq)*? — (Ce + C9), (6.10) 


¢, and c, being determined from the approximate values of H(n,u) for »=0-95 
and 1-00. In this way approximate values were also found for H(n.) (0 <p <1). 
Three cases were considered : 


A=", =%, Ag=ng=4; | 
A, =n, =}, da=me=45 } (6.11) 
A, =m, =}, Ag = Me = §. J 


In Cases I and II the initial values found by the method described above differed 
from the final values in the second decimal place; in Case III the difference 
only occurred in the third decimal place. 

The equations (4.2), on which the iterations were based, can be written 


BAL, (u’) dp’ 
bth’ 


oa ie _ 1/2 7 t wy (u’) dp’ 
Hap) (cyAq + %gAq)? + $a5(I OY yas erg 


I 


Fray =a ae) + deat) [ tw Halo!) de 


0 My /Ng+p' 


’ 


+ $09(I — Ag) | 
(tH, (u') du’ 


p+ hag(t—Aq) | Ea, 


(6.12) 


where H,(u)=H(n,u) and H,(u)=H(ngu). At each stage of the iterations the 
approximate values of the H-functions at ~=0, 0-1, 0-2, 0-4, 0°6, 0-8 and 1-0 
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were represented by polynomials of degree six in x. This enabled the iterations 
to be carried out with the help of tables of the integral 


1 , m , 
p- (ee 


=I,2,... 
m 9 Pxtp’ (m I, , »7) 


for x=1,m,/n, and n,/n,. These were constructed from the recurrence relation 


I I+px 
In= = —Hlm1y Ty=In——. 


Let H, and H,™ be the nth iterates of the functions H,() and H,() and let 


A.) =H, - HH, Af» =H - Hy, 

so that these are the corrections to the (m—1)th iterates. After several iterations 
it was found that A,”/A,» and A,”/A,-», though alternating in sign, were 
practically independent of both m and ». When this state of affairs was reached, 
the remaining errors in the H-functions were almost completely removed by 
summing a series of extrapolated corrections as follows: 

Let B,= |A,/A,"? |, B.=|A,”/A,"-» |, and assume that f, and fy are 
independent of m. Then the equations 


H,(u) = Hy +A — By + BY... )= Ay + Ay/(1 + By), 


H,(u) = He + Ag — By + Ba? —...)= Hy” + AX” /(I+ ft in 


will give the H-functions exactly. Since, however, 8, and f, were not entirely 
independent of n, the equations (6.13) only gave much improved values for the 
H-functions. These values were re-iterated until further iteration made no 
alteration in the fifth decimal place. The extrapolated corrections were so 
effective that only one or two such iterations proved to be necessary. The final 
values of the H-functions are given in Table I. 


TABLE I 
Values of the H-functions for a doublet 





I: 9,=10, 72=5 II: =4, 72=2 III: »,=1, 73=4 


Ai,(p) Aiy(p) Hu) Ai,() A{u) 








1 ‘0000 0000 
1*1440 2278 
1°2429 ‘3745 
1°3975 ‘5893 


“0000 10000 
“1140 1°1647 0598 0724 
‘1869 1°2620 0935 ‘1108 


10000 I 
I I 
I I 
2936 1°3944 1°1382 I‘I59Q1 
I I 
I I 
I I 


“0000 


I°5197 7462 3722 1°4844 ‘1679 -1896 
1°6211 8684 4341 1°5509 -1896 ‘2111 
1°7075 ‘9670 1°4844 1‘6020 2064 ‘2271 


I 
I 
I 
I 
I 
I 





























6.2. The residual intensities—The emergent intensities for the lines of the 
doublet are given by equation (5.3), and in the continuous spectrum outside 
the lines we have [“™ (0, u)=a+bu. The residual intensities, expressed as 
percentages of the intensity in the continuous spectrum, are 


r{u) = 1001,(0, n)/T™(0,n) (i=1,2). (6.14) 
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Values of r,(u) and r,(u) were calculated for Cases I, II and III for a region 
of the spectrum where b=$a. The first moments of the H-functions and certain 
integrals required in the calculations had already been evaluated as part of the 
iteration process described above. The results are given in Table II and are 
shown by the curves of Fig. 1, in which the residual intensities are plotted as 
functions of . Since, in the wings of the lines, 7 is proportional to (AA)-?, where 
AA is the distance from the line centre, the cases I-III of (6.11) lead, respectively, 
to cross-sections of the line profiles at increasing distances from the centres of 
the lines. This representation of the results is similar to that adopted by Houtgast 
in (3), where cross-section curves for multiplets in the solar spectrum are given. 

In Table III we give the corresponding exact values for a doublet without 
interlocking. ‘These values are shown by the curves in Fig. 2. They have been 
calculated from the exact solution in (2), Section 84.5, equation (69) and from 
Tables XI and XXXIII of (2). Appropriate interpolations were used when 
necessary. 

Tas_e II 
Residual intensities for a doublet with interlocking 





I: 9, =10, 72=5 Il: 1=4, 92=2 III: y,=1, n2=3$ 





r1(u) ro(4) ri (4) re(y2) ri() Yo(j2) 





’ 38°57 43°69 56-83 64°02 83-60 89:06 
39°09 46°39 56°88 65°74 82°58 88-80 
38°16 46°35 55°17 64°76 80°24 87°28 
35°86 44°73 51°49 61°81 75°97 84°39 
33°70 42°70 48:28 58-95 72°60 82:07 
31°82 40°75 45°60 56-50 69°96 80°26 
30°18 39°00 43°37 54°41 67-86 78-80 





























Taste III 
Residual intensities for a doublet without interlocking 





I: 9,=10, ne=5 II: 7,=4, 72=2 III: n,;=1, n2=$ 





r(“) ro(H) r(#) ro() ry() ro(}2) 





35°34 48°74 53°55 68-79 82°03 90°97 
36°58 49°71 54°36 68-90 81°44 90°08 
36°12 48°71 53°15 67°07 79°34 88-24 
34°51 46°03 50°12 63°18 75°36 85-00 
32°77 43°45 47°30 59°85 72°16 82°49 
31°17 41°21 44°88 57°12 69°63 80°56 
29°72 39°26 42°81 54°86 67°56 79°03 





























From Tables II and III and from Figs. 1 and 2 it will be seen that the effect of 
interlocking is to increase r,() and to decrease 7,() for all ». Let d(u) denote the 
difference r,() — 7,(4) in the interlocked case, and let d,() denote this difference 
when there is no interlocking. Then d(u)<d,(u) for all ». The values of 
d(u) and d,() in Cases I-III are shown by the full and broken curves respectively 
in Fig. 3. Near the centre of the disk d)() —d(y) is very small, but it increases 
towards the limb in every case and for the sections nearest to the centres of the 
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lines (Cases I and II) it rises to more than 4 per cent of the continuous spectrum 
at w=0°2. 

An interesting feature of the curves in Figs. 1 and 2 is the drop in 7,(u) and 
r() near to.=0in CasesI and II. This occurs both when there is interlocking 
and when there is no interlocking. This result is, however, only of limited 
interest, since the plane-parallel representation of a stellar atmosphere only applies 
for values of ~ >0-2. 








ry) 
90 


30 30 
l N l | I | N 

lO 08 O68 04 O2 OO lO O8 O06 04 O2 OO 
Centre kK Limb Centre ‘. Limb 




















Fic. 1.—The residual intensities for a Fic. 2.—The residual intensities for a 
doublet with interlocking. doublet without interlocking. 


Both figures show cross-sections of the line profiles at three distances from their centres, the full 
curves representing the stronger component and the broken curves the weaker component. 


6.3. Eddington’s approximation.—In this section we consider the accuracy 
of an approximate solution for the residual intensities for interlocked doublet 
lines based on Eddington’s approximation. ‘The details of the solution for the 
residual intensities in the emergent flux are given in (5). ‘The residual intensities 
in any direction and in a region of the spectrum where b= a are given by 


Ay, +%,4 %,Ai.+%.A4 
r(u)=100(2-+ uy 4x43 Be +n( ce Rts mn Aen eee a =), 
B + ee + teAe, aAre + “st) } 
T+ *\I+m+Hp, | +7 t+Hp2 J” | 








ru) =100(1 + $x)? {x + : 
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In these expressions p, and p, are, respectively, the larger and smaller positive roots 


of the equation 


P*— (Gir? + Ga2")P? + 9117420” — 912741” = 0; (6.16) 
where 


Gra” = 393(I +71)%2, Gar” = 392(I + 2)%4, (6.17) 
ux? = 3(1+m){I+7(I—%)}, Goao®=3(I +2){I + 7a(I — H2)}. 























7 





0:8 ; 
Centre Limb 
Fic. 3.—The variation of d(u)=re(u)—1r,(u) from centre to limb. The full curves represent 
the variations for a doublet with interlocking and the broken curves those for a doublet without 


interlocking. 


The A’s are given by the simultaneous equations 


(6.18) 


(I++ $f) Ant (2+ + $P2)412 + =9, 
(I++ $P))Ao, + (1+ No+ $P2)Acg+2=0, 


re 22° — p" A — u—Ps* 4 6 
ir mee 22= > “712 (6.19) 
q21 Tie 


The residual intensities for Cases I-III, calculated by means of equations 
(6.15)+(6.19), are given in Table IV. Comparing these approximate values with 
the exact values given in Table II, the proportional errors in the residual intensities, 
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although large at the extreme limb, are found to lie within +1 per cent for values 
of ~>0o-2. The error in d(u)=7,(u)—7,(u) is likewise small for ~>0-2, and 
when »=0-2 it only amounts to 0-5 per cent of the intensity in the continuous 
spectrum. Thus the errors of the approximate solution are satisfactorily 
small in the range of values for » for which the plane-parallel model applies, and 
they are sufficiently small to predict the limb effect of Fig. 3. 


TABLE IV 
Residual intensities for interlocked doublet lines according to Eddington’s approximation 





I: ,=10, 72=5 II: 7,;=4, n2=2 III: 4,=1, n2=$ 
r1(4) ro(14) r,(4) ro(4) r1() ro(1) 


42°35 47°16 61°53 . 67°94 88-06 92°04 
40°38 46°97 58°43 66-50 84:21 89-69 
38-64 46°27 55°77 64°82 81:06 87-67 
35°74 44°31 51°42 61°53 76-21 84°47 
33°40 42°26 48°04 58-65 72°66 82:07 
31°46 40°36 45°32 56-22 69°96 80°23 
29°84 3865 43°09 54°18 67°83 78°77 



































7. Conclusions.—The theory and calculations show that there is a limb 
effect (Fig. 3) in the inner wings of interlocked doublet lines, the difference 
d(u) =r.(%)—1,(u) between the residual intensities of the two lines at the same 
distance from their centres being smaller than for doublet lines without inter- 


locking. The limb effect for 1 >0-2 is also predicted by an approximate solution 
of the equations of transfer by Eddington’s approximation. This suggests that 
results obtained in this way for the centre-to-limb effects in more complicated 
multiplets are probably reliable for »>0-2 and that conclusions drawn from the 
approximate solution are likely to remain true for the exact solution. 





APPENDIX 


8. The principle of reciprocity for the scattering functions.—In Section 3 we 
have shown that the scattering functions S,,(u, y;) (r, s=1I,..., %) satisfy the 
equations (3.2), viz. 


I+”, I+”, Sy o(H, Hi) ze : [ ' dy’ 
| pb + Bi; ] I-e mn Se Bn! Xp role Bo 
: : , dp 
x {n+ x ata | So(u » Hy) ae (8.1) 
q=1 0 1d 
where r,s=1,...,, and it is known from the theory of coherent scattering that 


Si (Ms Hi)=Sr(on He) (r=1,.-+,h). (8.2) 
We shall now deduce (3.4) from (8.1) and (8.2).* 
* In (1), the proof of the principle of reciprocity, which was omitted, is very similar to that 


given here, but it ends with an appeal to the theory of homogeneous integral equations instead of 
to the theory of homogeneous linear equations. 
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Put s=r in (8.1) and interchange p and yw; Then it follows from (8.2) that 


k ‘1 dp’ ~ fs du’ 
y he oll av , a, 
( ' ont 7 J, Srolts #’) “ { ” ont “aq 0 Sarl's Hi) a 


t ny dp’ ¥ en 8 
- {x ~ © ap | Seale 4} {n,+ baad aang | Sor('s ) ae 
This is of the form 


fees) =Aedg(), 


Sv )/g(+) =fls)/g(+s) =a constant independent of . 


The constant does, however, depend on r and we therefore write 
~ A » Gp -. - ae 
Net fad %oNq | P Solu MoT =C, {x a * Ly J. Syy(Hs oe ye} . (8.3) 
From (8.1) and (8.3) we have 
I+, | I+ me] S,6(Hs 4) { hee 'g » e'7 
+ =C,41+ = | Sele, pt’) 3 
L Me My I—e Poa 0 ae 
S ri du’ 
z ae 
x {x + 2% \, So(Hir # ra , 
and interchanging 7 and s, » and p,;, we get 


I+, 1 <2] Ss (His #) { . f j =} 
+ =C,jI+ XT a S ‘9 ae 
| by im  — r ae fe ep(His & ) iL 


k 4 du" 
x {r+ Xap | Seals’) \. 


> , 
P - 


whence 








From these we have 
C,S,,(u, Hi) wa CYS 5A Mir lt). 
By (8.4), (8.3) can be written (changing p into q) 


: 4 ‘t , \ ap’ 
(n,-—C,)+ — (ng—Cq)%q | Sale ey =0, 
q= . é 


where r=1,...,&. This gives k homogeneous equations for the constants 
(n,—C,) (r=1,...,k). The determinant formed by the coefficients is a function 
of » and is not identically zero; hence the only solution independent of y is 


"e—-C,=0 (r=I,...,R). 
Thus (8.4) becomes 

ArSr(Ms Hi) = NS or(Hir )s 
and this is (3.4). 


St. Hugh’s College, University Observatory, 
Oxford: Oxford. 
1953 July 9. 
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MOTIONS IN THE SUN AT THE PHOTOSPHERIC LEVEL 


IV. THE EQUATORIAL ROTATION AND POSSIBLE VELOCITY FIELDS 
IN THE PHOTOSPHERE 


A. B. Hart 
(Communicated by the Director, University Observatory, Oxford) 


(Received 1953 November 28) 


Summary 


From spectra of the equatorial photosphere taken near the limb of the Sun 
between 1950 October and 1952 June, 480 determinations of line-of-sight 
velocity are made. A careful investigation shows, contrary to an earlier 
suggestion, that the effect on velocities, of light scattered in the Earth’s atmo- 
sphere and in the telescope, is negligible. In order that each day on which 
observations are made should carry equal weight in the final solution, only 
334 of the determinations are used to find the mean equatorial velocity. 
After the elimination of limb effect, the velocity due to solar rotation is 
found to be 1-91 km sec™! for the mean epoch 1951°5. The individual 
determinations have a large dispersion about this mean and a statistical 
discussion shows that the chance is less than one in a thousand that the 
dispersion is due to errors of measurement. It is therefore highly probable 
that the velocity varies from point to point in the photosphere. The extreme 
range in measured velocity is approximately 0-3 km sec™! and successive 
maxima or minima of velocity appear to be separated by a distance of the 
order of 75 000 km. Some evidence is found for bands of constant velocity 
occurring perpendicular to the solar equator, but this is less certain. The 
existence of such velocity changes may account for the discrepant rates of 
rotation found by earlier observers. 





Introduction.—F rom visual measures of Doppler displacements, Dunér (§) in 
1888 made the first spectroscopic determination of solar rotation. The potential 
accuracy of the method and its applicability to all solar latitudes gave grounds for 
the belief that the determination of rotation by sunspots would soon be superseded. 
Discrepant results were, however, obtained by photographic spectroscopy 
between 1906 and 1913, and since 1915 investigations have been restricted to a 
determination of equatorial velocity ‘in the hope of elucidating the causes of 
these discrepancies. 

Adams (2) at Mt. Wilson between 1906 and 1908 made the first photographic 
measures and laid the foundations of a method used by many later investigators. 
Working in the wave-length region 4250 A witha solar image 170 mm in diameter 
and a dispersion of 0-71 or 0-56 Amm, he used two small prisms to bring light 
from opposite limbs of the Sun into juxtaposition on the slit of the spectroscope. 
Half the displacement between the same spectral line on the two limbs gave the 
value of the rotation. ‘The spectra were measured with an ordinary cross-wire 
micrometer. 

This procedure was improved by J. S. Plaskett (19), at Ottawa in 1911, who 
used a triple-prism system. It produced spectra arranged so that the spectrum 


2 
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from one limb lay between the two spectra from the other, thus removing spurious 
displacements caused by any non-parallelism of the measuring wire and spectral 
lines. A series of measures in 1911, 1912 and 1913 was carried out in this way (18) 
and a precisely similar method was adopted at Allegheny by Schlesinger (20). 
It was clear, however, that the small-prism method was still subject to systematic 
error in that displacements of the prisms might cause them to illuminate the 
collimator differently and so produce spurious shifts in the spectral lines. By 
returning to Dunér’s original method of using water-vapour lines as standards of 
zero velocity, H. H. Plaskett (13) eliminated this source of error which, unless 
such a precaution is taken, vitiates any determination of velocity. 

With these improvements, in 1915, a variation in the equatorial velocity with 
amplitude + 0-15 kmsec"! was shown to exist. Of the reality of the measured 
change in velocity there can be little doubt, but the change itself was very differently 
interpreted. On the one hand, it was ascribed by H. H. Plaskett to an actual 
change in velocity of the photospheric matter and, on the other hand, by R. E. 
DeLury (3) to the superposition of a variable amount of sky spectrum on the limb 
spectrum. 











| | 


INO 4920 





Fic. 1.—The rotational velocity at the equator. 
O W. S. Adams and F. B. Lasby, Mt. Wilson. @ Mt. Wilson. 
C) ¥. S. Plaskett and R. E. DeLury, Ottawa. @ H. H. Plaskett, Ottawa. 
A F. Schlesinger, Allegheny. @ Edinburgh. 


For many years systematic observations were made on a few days in each year at 
Mt. Wilson (12) and at Edinburgh (9, 21, 22), using atmospheric lines as standards, 
in an effort to obtain a good mean value for the equatorial velocity. The values 
obtained up to 1935 are shown in Fig. 1. Results from Mt. Wilson and Edin- 
burgh are averaged over the five year periods beginning in 1914, I919, 1924 and 
1929. Determinations using atmospheric standards are shown as filled symbols 
(the only ones shown after 1914) and the values obtained by some of the earlier 
workers are shown as unfilled ones. It will be noted that once atmospheric lines 
were used as standards, the measured equatorial velocity settled down to 2 mean 
of 1-95 kmsec~! with a dispersion +0-05kmsec-4. This dispersion has been 
variously attributed to the effect of superposed sky spectrum, to the vertical 
apparatus function of the spectroscope and to actual changes in velocity, but no 
analysis has been made of the relative sizes of these effects. 





No. I, 1954 Motions in the Sun at the photospheric level 19 


Consequently, new observations have been made at Oxford with the intention, 
(1) of obtaining a value of the equatorial velocity together with a critical estimate of 
its reliability, (2) of determining precisely the effect of superposed sky spectrum 
on the measured velocities, and (3) of making some investigation of any irregu- 
larities in velocity which may occur. The following account of an investigation 
carried out between 1950 and 1952 is divided into four sections. The first section 
concerns the observations and the method of measurement and the second the 
effect of superposed light on measured velocities, a discussion of which is necessary 
before measured velocities can be converted into velocities on the Sun. The third 
section contains the details of the calculation of the equatorial velocity and a value 
of the limb effect, and the last section is a discussion of the residuals in the measured 
velocities. 


1. Observations and measurements 

1.1. Observations.—The observations were made with the Oxford solar tele- 
scope (14). The telescope is of the vertical Cassegrain type and produces at its 
focal plane an image 180 mm in diameter giving a linear scale 1 mm=10"'4 arc. 
The image is formed in the plane of the spectroscope slit (15). The slit lies very 
nearly in the East-West direction on the image and the precise direction is 
determined for each series of observations. The position on the image at which 
the observations are made can be found and accurately maintained throughout an: 
exposure by means of the rocking plate described by H. H. Plaskett.* The 
spectra were taken in the wave-length region 5900 A where there is a number of 
suitable unblended solar and water-vapour lines. ‘The dispersion is 1-5 A mm“. 

Spectra of the two limbs were photographed separately, any one plate carrying 
two (or in the third series, six) exposures of the same limb. The height of the 
spectroscope slit was limited to 12-5 mm but the solar image did not occupy the 
whole height. Each spectrum was only 6 to 10 mm high so that the limb proper 
could be identified on the plate. A metal diaphragm just in front of the dark slide 
of the plate prevented the light distributed in the vertical apparatus function of the 
spectroscope during one exposure falling on any other spectrum. 

Three series of limb spectra were taken. The first, an exploratory series, 
consisted of repeated exposures along that chord at each limb, which passes through 
the equator atthe limb proper. Inthe second series, successive exposures covered 
an appreciable area of the disk and followed each other as rapidly as possible. The 
maximum time which elapsed between the first and last of a set of six exposures on 
any one limb was 10 minutes. Adjacent spectra were separated by 1° latitude at 
the limb and were arranged so that approximately equal areas North and South 
of the equator were covered. Fig. 2 shows the arrangement of the slits in this 
case. The third series was similar to the first. For high resolution on the solar 
image it is necessary to have good seeing. ‘The second and third series had, in 
general, seeing I or 2 on a scale of 0 to 5 (5 the best, see Paper I), but series I was 
taken during a bad season and is inferior in this respect. Of the numerous spectra 
so obtained, three sets made in the best available conditions were selected for 
measurement. Details are given in Table I. A spectrum of the centre of the 
Sun, 12-5 mm in height, was taken for use as the standard spectrum in the 
measuring machine. Ilford Rapid Process Panchromatic plates were used for all 


* This is the first paper in the series, Motions in the Sun at the photospheric level. Papers 
in this series will be referred to by their serial numbers, I (16), III (11). 
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the work; they were developed in Kodak Press Contrast Developer or Parkhurst 
two-solution developer. 

In Series I and II, on each day that a set of limb spectra was obtained, at least 
one “‘ sky brightness ”’ plate was taken to find the ratio of the intensity of the light 
at a fixed distance outside the limb to that at a fixed distance inside. Each plate 
consisted of one exposure at the centre of the Sun through a Zeiss calibrated step- 
wedge and one without it to show any variation in illumination along the slit; 
then narrow spectra, 3 mm high, were taken of the light from points 40” arc radially 
inside and outside the limb. The exposure time was the same for each of the four 
spectra and suitable densities on the plate were obtained by using N.P.L. calibrated 
filters. The plates were not developed until 24 hours after they had been taken, so 








SOLAR EQUATOR, w 


SOLAR AXIS. 














Fic. 2.— Arrangement of observations in Series II. (The diagram is not drawn to scale.) 


that the initial fading of the photographic images did not introduce any spurious 
differences in density on the plate. Table I, col. (8) shows the ratio of the bright- 
ness of the sky spectrum as a percentage of the brightness of the limb spectrum, 
as found in this way. 

1.2. Measurements.—Spectra are marked perpendicular to the dispersion, at 
either end of the spectral region, by a needle attached to a Hilger measuring 
machine. This determines to approximately 0-01 mm the position on the spectrum 
and hence the position on the Sun at which measurement is to be made. The 
height of spectrum measured is defined by a limiting diaphragm of variable 
aperture in the eye-piece of the Evershed machine(Paperl). Apart from personal 
equation, error in measurement is chiefly produced by irregularities in the 
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distribution of silver grains. An investigation showed, as one would expect, that 
the consistency of measurement increases as the height of the limiting diaphragm 
increases, up to a certain limit. This limit is approached asymptotically and 
occurs at about a height of rmm. __Increase of diaphragm height decreases, 
however, the resolution on the solar image and in Series II and III the spectra were 
measured with a diaphragm aperture of 0-75 mm at intervals of 0-75 mm so that 
the entire spectrum was measured between the first and last marked points. In 
these series the first point of measurement was I mm inside the limb. In the 
exploratory Series I the diaphragm aperture was 0-3 mm and the measurement was 
madeatintervalsofi1mm. The first point of measurement was a variable distance 
from the limb, depending on the spectrum, the average distance being 1 mm. 


TABLE I 





Observations Measurements 





(3) | (4) (9) | (10) | (12) 





ons 9 W.V. 
Series I 9 sol. 
Exploratory 
1950 Oct. 5 i 9 W.V. 
9 sol. 


36 








36 





Series Ila 
1951 April 24 











Series IIb 
1951 May 25 














3h 46™ 


Series III 
1952 April 27 








3h 17m 












































Col. (1) Limb; (2) No. of sp. on each limb; (3) Exp. time in secs; (4) Slit width in mm; 
(5) Spectrum height in mm; (6) Time interval between first and last spectrum; (7) Seeing 
on a scale of 5; (8) Percentage sc. light; (9) No. of pts. per spectrum; (10) No. of lines 
per point; (11) Total no. of measured velocities. 
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the work; they were developed in Kodak Press Contrast Developer or Parkhurst 


two-solution developer 

In Serves | and II, on each day that a set of limb spectra was obtained, at least 
one “ sky brightness " plate was taken to find the ratio of the intensity of the light 
at a fixed distance outside the limb to that at a fixed distance inside. Each plate 
conmsted of one exposure at the centre of the Sun through a Zewms calibrated step- 
wedge and one without it to show any vanation in illumination along the slit; 
then narrow spectra, ; mm high, were taken of the light from points 40° arc radially 
inside and outside the limb. ‘The exposure time was the same for each of the four 
spectra and suitable densities on the plate were obtained by using N.P.L. calibrated 
filters. he plates were not developed until 24 hours after they had been taken, so 

op 


| 


SOLAR AXIS, 














Fic. 2.—Arrangement of observations in Series II. (The diagram is not drawn to scale.) 


that the initial fading of the photographic images did not introduce any spurious 
differences in density on the plate. Table I, col. (8) shows the ratio of the bright- 
ness of the sky spectrum as a percentage of the brightness of the limb spectrum, 
as found in this way. 

1.2. Measurements.—Spectra are marked perpendicular to the dispersion, at 
either end of the spectral region, by a needle attached to a Hilger measuring 
machine. This determines to approximately 0-01 mm the position on the spectrum 
and hence the position on the Sun at which measurement is to be made. The 
height of spectrum measured is defined by a limiting diaphragm of variable 
aperture in the eye-piece of the Evershed machine(Paper1). Apart from personal 
equation, error in measurement is chiefly produced by irregularities in the 
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distribution of silver grains. An investigation showed, as one would expect, that 
the consistency of measurement increases as the height of the limiting diaphragm 
increases, up to a certain limit. This limit is approached asymptotically and 
occurs at about a height of 1 mm Increase of diaphragm height decreases, 
however, the resolution on the solar image and in Serves I] and I 11 the spectra were 
measured with a diaphragm aperture of o-75 mm at intervals of O-75 mm so that 
the entire spectrum was measured between the first and last marked points. In 
these series the first point of measurement was 1 mm inside the limb. In the 
exploratory Series | the diaphragm aperture was 0 } mm and the measurement was 
made at intervalsof 1 mm. The first point of measurement was a variable distance 
from the limb, depending on the spectrum, the average distance being 1 mm. 


Taste ! 
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Col. (1) Limb; (2) No. of sp. on each limb; (3) Exp.time in secs; (4) Slit width in mm; 
(5) Spectrum height in mm; (6) Time interval between first and last spectrum; (7) Seeing 
on a scale of 5; (8) Percentage sc. light; (9) No. of pts. per spectrum; (10) No. of lines 
per point; (11) Total no. of measured velocities. 
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The relation between distance along the slit and position on the solar image is 
known when the position of the limb on the plate is known. The problem of 
finding the position of the limb proper was dealt with in two ways. In the first 
place, in Series I, microphotometer tracings were made perpendicular to the 
dispersion. ‘The intensities were found from the calibration curves of the “‘ sky 
brightness ” plates, which are described above. The true position of the limb is 
then found by comparison with Wanders’ distribution (23) for suitable values of 
the scattering parameter. In the second place, the position of the limb on the 
plate was estimated visually with a cross-wire ; this method becomes more accurate 
astheseeingimproves. For Series I, the visually estimated positions differed from 
the true ones on the average by approximately 0-03 mm and seldom by more than 
o-Imm. The first method takes a prohibitively long time and the average 
accuracy attained by the second method is comparable. A visual estimate of the 
position of the limb was made in Series II and III. 

Nineteen lines were finally selected and measured at each point on the plates of 
Series Il and III (see Table II). Ten of these are solar lines and they are bracketed 
by 9 water-vapour lines of comparable intensity which are used as standards of 
zero velocity. In this way, any errors caused by non-uniform illumination of the 
collimator are avoided. 


Tas_e II 
Water-vapour lines Solar lines 


Wave-length Rowland Wave-length Ele Rowland Km sec~' per micron shift 
in A intensity in A ement intensity (All other 

(Series I) measurements) 

5901468 
5919°055 
5919°643 
5924°272 
5932°093 
5946-006 
5947 066 
5957880 
5958-622 


5905 680 Fe I 0°0744 0'0746 
*5909°981 Fe I ae 0748 
5916-258 Fe I "0749 ‘O751 
5927°796 Fe I ‘0754 ‘0757 
5929 682 Fe I 0755 0758 
5930°190 Fe I 0756 0758 
5934665 Fe I 0758 0760 
5948548 Si I ‘0765 ‘0767 
5952°726 Fe I 0767 0769 
5956-705 Fe I ‘0769 ‘0771 


ae NWUhNU D 


(Wave-lengths taken from Rowland’s Revision and corrected according to I1.A.U. 
recommendations.) 


* Not measured in Series I. 


The spectra were measured at the various points along their height (see Table 1) 
by the Evershed method. _ Its principal advantages are the elimination of spurious 
shifts caused by the asymmetrical distortion of the profiles by the horizontal 
apparatus function of the spectroscope, and the use of an automatic recording 
device. The machine is fully described in Paper I and all measurements were 
made with it, except the direct measures of Series I made on the prototype 
machine. : 

General questions concerning the accuracy of the method are discussed in 
Paper I. A comparison between the velocities obtained by two measurers for 
direct and reversed measures made at 7 points on one of the spectra of Series IIb 
showed that the large R-D discrepancy of 0-7 km sec~? is inherent in the method, 
though it should disappear for water-vapour and solar lines of the same intensity. 
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Such a large value of the R-D difference is objectionable but the mean velocity 
agreed well with the mean of direct and reversed measures made with a cross-wire 
on the same spectrum. The root mean square error for a single line for these 
Evershed measures was + 0:23 km sec™}, which corresponds to a displacement of 
3 microns on the plate. ‘This may be compared with the values of I-54, 24 and 
I-5u obtained at Mt. Wilson 1906-8 (2), Ottawa 1911-13 (18) and Allegheny 1911 
(20), respectively. ‘The smaller dispersion of these earlier differential measures 
may, however, be a result of prepossessions (17) which are eliminated by the 
automatic recording device on the Evershed machine. 

A test was also made of the consistency with which a measurement is repeated, 
by measuring the 84 spectra of Series Ila twice. Table III shows the distri- 
bution of velocity differences which was obtained. All the velocity differences 
are small (which is satisfactory); the average arithmetic difference between the 
two measures of a single point is 1-420-:1kmsec"!. This is comparable with 
the root mean square error of measurement of a single velocity point, given by 
+ 0°23/1/10= +0-08 km sec“, since 10 solar lines were used to give the velocity 
at any point. The average algebraic or systematic difference between the two 
measures of a point is 0-24, which is negligible. 


TABLe III 
Limits of the No. of velocities 
difference interval in the interval 
Apu Av 
microns km sec™! 
4°0, 30 o'31, 0°23 
3°90, 2-0 0°23 O'ls 
2°0, 1°O O'r5s, 0°08 
‘0, °o'o 0-08, 0°00 
0:0, —I°0 0°00, —o-08 
—FO, —2°O —0'08, —O'r§ 
—2°0, —3°0 —O'I5, —0°23 
—3°0, —4°0 —0'23, —0°3I 
a a 0°38. 0-398 
—50, —60¢ —0°38, —0°45 


Altogether 480 measurements of velocity were made (see Table I). Before 
they can be used in the calculation of the velocity of solar rotation it is necessary 
to find the effect of any extraneous light superposed on the spectra. 


2. The effect of scattered light 

The spectrum of the sky remote from the Sun is an integrated spectrum of 
the whole disk with a mean velocity corresponding to that of the solar centre. 
If such a sky spectrum be superposed on the spectrum of the limb of the Sun, 
the resulting blended spectrum will have a mean velocity displacement less than 
that of the limb alone and the measured rotation will be distorted by an amount 
approximately linearly proportional to the ratio of the brightness of the sky 
spectrum tothelimb spectrum. This hypothesis, which ignores the mechanism of 
the scattering by which the sky spectrum is produced and which was apparently 
never tested, was put forward by DeLury (3, 4) to account for the varying values 
of solar rotation found between 1906 and 1915. 
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One quantitative test can, however, easily be made. The spectrum of the 
sky superposed on the limb must be that of the sky immediately adjacent to the 
limb, i.e. that of the solar aureole. If any distortion of the measured velocities 
is to occur, the velocity of the sky spectrum must be different from that of the 
limb. Accordingly, on 1952 September 18 a series of exposures of the West limb 
of the Sun was made with a rapidly rotating sector, having a 12° aperture, to cut 
down the flux from the disk to an amount approximately equal to that of the 
flux from the inner aureole. The height of the slit was limited to 12-5 mm and 
the image was guided so that the limb fell well within the portion of the slit 
shielded by the sector, while the aureole, at distances greater than a few seconds 
of arc from the limb, fell on the unshielded part and was fully exposed. Five 
of these exposures were made and a narrow spectrum of the centre, 3 mm in 
height, was taken before and after the series. 

















Fic. 3.—Velocity distribution in the aureole. 


Ordinates: velocities differenced from the centre velocity. 
Abscissae: distances from the limb along the E.-W. chord of the solar image. 


The spectra so obtained were measured on the Evershed machine with a 
diaphragm height of 0-75mm. Inside the limb the spectra were measured 
at distances of I mm, I-75 mm and 2°5 mm from the limb. Outside the limb, in 
the spectrum of the aureole, the velocity was again measured at 3 points separated 
by 0-75 mm intervals. (The first position at which the fully exposed spectrum 
filled the diaphragm aperture depended on the guiding of the limb.) The 15 
measured velocities in the aureole were divided into four groups and the mean 
velocity and position calculated for each group. For the centre spectra, a 
diaphragm aperture of 3mm was used. ‘The resulting velocities are shown by 
circles in Fig. 3; in each exposure the velocity at the centre was taken as the 
standard velocity. ‘The length of the vertical line through each circle shows the 
root mean square error of the velocity. 
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It will at once be seen that the basic postulate of DeLury’s hypothesis is not 
satisfied. The measured velocities in the aureole, instead of being zero, drop 
steadily from the velocity of the limb at the limb to reach zero only at a distance 
of several millimetres from the limb. Since at the places in the aureole where the 
velocity is measured the vertical apparatus function of the spectroscope accounts 
for less than one-thousandth of the observed brightness, the results shown in 
Fig. 3 are essentially due to the aureole alone. Hence if, as in DeLury’s hypo- 
thesis, the spectrum of the aureole exactly at the limb were superposed on the 
spectrum of the limb there would be no distortion of the measured velocity at all. 

The processes which occur are, however, more complex than pictured by 
DeLury. Scattering (as opposed to specular reflection) by the telescope mirrors 
and scattering by the Earth’s atmosphere redistribute the light from each point 
of the Sun over an area in the focal plane of the telescope. By far the greater 
part of the scattering takes place over small angles and only a minute portion, 
estimated at 10-* by van de Hulst (8), takes place over large angles to form the 
sky spectrum remote from the Sun. What is observed in the aureole near the 
solar image is presumably this small-angle, redistributed light. If this be the 
case, the main effect of scattering on line-of-sight velocities will be to smear out 
any differences between them and only to a completely negligible amount, as a 
consequence of second and higher order scattering in the Earth’s atmosphere, 
to superpose an integrated zero velocity spectrum on the limb spectra. 

T. D. Kinman has measured the surface brightness of the aureole as a function 
of its radial distance from the limb (Paper III). His measures are in good 
agreement with each other and with the surface brightness measures made at 
a single point (Table I) in the present investigation. From these measures 
Kinman has calculated the scattering function which must prevail if the aureole 
is wholly produced by small-angle scattering. If, as we have suggested above, 
this assumption is true and if there is a negligible amount of superposed zero 
velocity sky spectrum, then Kinman’s scattering function should predict the 
velocities which are observed in the aureole. Details of this calculation are 
given in the Appendix and the velocities predicted by the scattering function 
are shown by the continuous curve in Fig. 3. When it is recalled that Kinman’s 
scattering function refers only to average Oxford conditions, the agreement 
between the velocities predicted by the function and those actually observed 
may be considered not unsatisfactory. The difference between the observed 
and calculated velocities is in the opposite sense to the difference which would 
arise if there were any appreciable quantity of zero velocity light superposed on 
the small-angle redistributed light in the aureole. It would therefore appear 
that small-angle scattering is the dominant mechanism and we need only consider 
the effect of such scattering on the observed line-of-sight velocities. 

The order of magnitude of the effect can readily be calculated (see Appendix) 
from a simple example. If the equatorial region of the Sun has a uniform 
rotational velocity of, say, 2-07 km sec“, then after correction for the Earth’s motion, 
the observed line-of-sight velocity at the limb would be 2-07kmsec"!. The 
line-of-sight velocity observed (after correction for the Earth’s motion) as a 
result of small-angle scattering from all other regions of the Sun where the line- 
of-sight velocity is smaller would be 2-03km sec~!. Hence, at the limb, small- 
angle scattering produces an error of 0-04km sec-!. The error decreases as the 
point at which measurement is made moves further in from the limb and at a 
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distance of 10” arc from the limb (1 mm on the Oxford solar image) it would only 
be 0-03 km sec-4. Since the present observations refer to points 10” arc or further 
from the limb, the maximum systematic error, produced by scattering, will not 
be greater than 0-03 kmsec™!, that is, less than one-half the measurement error 
of 0-08 km sec~! for a single measured velocity. 

The nature and effects of scattered light are now clear. Small-angle scattering 
is the predominant mechanism and its only effect is a smearing out of velocity 
differences. Any systematic effects are well below the errors of measurement, 
so that measured velocities can be and are used in the subsequent calculations 
without further correction. 


3. Solar rotation and limb effect 


Reference to Table I shows that the 480 measured velocities are divided 
among the various series in the following way: Series I contributes 72 points, 
each day in Series II 84 points and Series III 240 points. In making a solution 
for the best values of solar rotation and limb effect for the mean epoch of the 
investigation it is advisable that approximately equal numbers of velocities 
should be used from each of the four days for which measurements were made. 
Accordingly, 4 representative spectra were selected from each limb of Series III 
so that 96 velocities only from the series are used. ‘The total number of velocities 
included in the final solution is, therefore, 336. 

The conversion of positions on the image to positions on the Sun and of 
measured line-of-sight velocities to velocities on the Sun are fully explained in 
Paper I. As shown in the last section, the effects of superposed light are 
negligible so the measured velocities only need correction for the relative motion 
of the Earth and Sun. Let this correction to the limb spectra be V,, and the 
correction for the standard centre spectrum be V, (V, includes a correction for 
a small line-of-sight component of solar rotation at points other than the centre 
of the image). Then the line-of-sight velocity on the Sun at any point is 
V+V,-V3,. 

For each point the equation of condition may be written 


(U+u) cosy, +vcos y,+wcos y3=(V+V,—V,)—A, (1) 


asin PaperI. U is the value of solar rotation, A the limb effect, and u, v, w are 
the components of any local velocity at the point, parallel to axes along the line 
of latitude, the line of longitude and the solar radius, which make angles +,, yo, y3 
respectively, with the line of sight. In a limited latitude zone, U the linear 
velocity can be represented by the Faye formula U =(a+6cos? B)cos B, where B 
is the heliographic latitude. Hence the complete equation of condition takes 
the form 


(a+bcos? B) cos Bcos y, +ucosy,+vcosy,+wcosy,;=(V+V,—V.)—A, (2) 


having six unknowns, a, d, u, v, w, A. 

Since u, v, w are the components of a local velocity there may be no relation 
between them from point to point, so the equations cannot be grouped together 
for a solution by the method of least squares. Furthermore, for the range of 
latitude covered by the observations | B|<6°, the coefficients of a and 6 are so 
nearly equal that the equations are ill-conditioned and a least-squares solution is 
inadvisable in any case (7). Consequently, a different method of solution has 
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to be found. The equations are therefore simplified and solved by an iterative 
method. 

Assume that the velocity is independent of latitude and u, v, w are zero. 
If U* be the equatorial velocity of rotation, the simplified equation of condition 
for any point becomes 


U*cos Bcos y;=(V + V,—V,)—A, (3) 


an equation only involving the two most important unknowns, where U* is 
constant and A is a function only of the fractional distance p/p, from the centre 
of the disk. 
Rewriting this in the form 
 V+t+V,-V, A 
U*= Ce: ’ (4) 


cosBcosy, cosBcosy, 





a first approximation to U* is easily found. ‘The numerical value of cos B cos y; 
depends mainly on the fractional distance from the centre of the disk, but the 
sign is opposite in the Western and Eastern halves of the disk. Hence, providing 
there are equal numbers of observations on the East and West limbs at each value 
of p/po, a condition roughly satisfied in practice, the first approximation to U*, 
say U,*, is given by 


V+V,-V, i 
a hie’ ae arate. ase (5) 


The resulting value of U,*=1-90 km sec" is used to give the first approximation 
A, to A, from 

A, =(V + V,—V.)— U,* cos B cosy, (6) 
for each point. For the 334 points lying in the interval 0-90 <p/py<0°995, the 
results so obtained are grouped together at twelve standard positions shown in 
Table IV (x=1—p/p, is the fractional distance from the limb). From these 
results, appropriately weighted, a least-squares solution for a quadratic expression 
for A, is calculated, giving 


A, = 0°236 — 2:640x + 2°07 1x2, 0°005 <x <O-10. (7) 


TABLE IV 


No. of A F No. of 
pipe * points . P/Po points 
©'908 0:092 8 -0-01 km sec™! O'951 0°049 42 +-o'12 km sec“ 
‘916 «084 8 —o'o9 ‘959 ‘O41 4I +0°17 
‘923-077 8 +o°12 ‘967. 033 42 +0°14 
8 
8 


4; 


°931 ‘069 +Or12 ‘972 026 40 +014 
-938 «©6062 +0°15 ‘982 «018 42 +021 
"943 °057 43 +0:°06 ‘988 -o12 44 +0'20 


Values from this analytical solution are then used to give a second approxi- 
mation to U* in each case 
V+V,-V, Ay 


a oT , 
cosBcosy, cos Bcos y; 


2 





(8) 


. 


The mean of the 334 values of U,* is 1-91 kmsec}= U,* say. 
So small a change as 0-01 km sec~! in the approximation to U* will produce 
a negligible change in the approximation to A, so the process has converged 


rapidly and satisfactorily. We shall adopt U*= U,*, A=A,. The values of 
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limb effect may be compared with those of M.G. Adam(x). Her interferometric 
results, necessarily obtained with very low resolution on the solar image, show 
both a larger value of A and a steeper rise to the limb. However, her more 
recent work with comparable high resolution on the solar image yields results 
not in disagreement with values found in the present investigation. 

Finally, it is necessary to test the assumption of the invariability of velocity 
with latitude. The individual values of U,* are grouped into eight zones of 
latitude and the mean found for each. Fig. 4 shows the result. The root 
mean square error of each mean is represented by a vertical line and the number 
in brackets in each case shows the number of observations used in calculating 
the mean. Clearly there is no systematic change in velocity and so no correction 
is needed. 





Br 
4) ao 6O SD % Gb (23) 








4 1 4 1 


o ° $ e ¢ p < 





Fic. 4.—The equatorial velocity; its independence of latitude. 


Abscissae: heliographic latitude. 


Residuals from the measured line-of-sight velocities were then calculated 
from the final values of U* and A. ‘The root mean square error of a single 
solar velocity from 334 observations is +0-I8kmsec"!. The root mean square 
error of U* is +0-o1kmsec"!, so that the first approximation U,* lies within 
the range of the root mean square error of the final value. One may take, therefore, 
for the mean epoch 1951-5, the equatorial velocity of solar rotation 

U*=1-91 + 0-01 km sec“, (9) 
and the limb effect in the range 0-01 <x <0-10 (where x is the fractional distance 
from the limb) 

A = 0:236 — 2-640x + 2:071x%, (10) 
as being the best possible values from the material. 


4. Evidence of local variations in velocity 

In the previous section the mean equatorial velocity U* and the limb effect 
A have been found from 334 line-of-sight velocities measured in Series I, II 
and III (cf. Table I). We now wish to see whether the individual measured 
values of the equatorial solar velocity, denoted by U, differ significantly from 
each other. For this purpose attention is confined to the homogeneous set of 
measures from Series II (168 values of U) and Series III (96 values of U). The 
70 measures of Series I, based on different lines measured over a smaller height 
of spectrum (cf. Section 1), are not included in the present discussion, though 
inspection shows that they support the conclusions which are derived for 
Series II and III alone. 
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4.1. Statistical discussion of errors.—The 264 values of the rotational velocity 
U obtained from Series II and III are each the mean velocity from 10 solar lines 
at a particular place in the equatorial photosphere. We shall call U a spectrum- 
velocity. ‘These 264 values have a dispersion about their mean with a root mean 
square error for a single spectrum-velocity of +0-17kmsec~!(+0-18 km sec 
is obtained for the 334 spectrum-velocities considered in the previous section). 
Can so large a dispersion arise solely from errors of measurements ? 

To answer this question we must first make a precise estimate of the errors 
of measurement. Let the measured velocity from a single line be wu, the line- 
velocity. ‘Then the spectrum-velocity, U, and the root mean square error, s, of 
the line-velocity at any point are given by 


n n 
Um=Sa, #2 ——-YU~a)V, (11a) and (11b) 
ny nm—I} 
where n=10, since each spectrum-velocity is the mean of 10 line-velocities u,. 
Now the operation of measurement is essentially one of taking samples of n 
measures at a time from a hypothetical, infinite population of line measures at 
that particular point in the photosphere. The value of the error of meas:trement 
will then be the value o of the root mean square error of a line-velocity for this 
infinite population. This cannot be calculated directly from equation (11 b) 
but an estimate of its size can be made if we suppose that the infinite population 
of all line-velocities has a Gaussian distribution about a mean (which is also 
unknown). 
The probability of obtaining a line-velocity between u and u+ du from the 
population is 
P(u) du =(4/27.c0)-1exp[—(u—)?/207] du. (12) 


Hence the probability of obtaining a sample of z line-velocities between u, and 
u,+du,,..., u, and u,+du,,..., u, and u,+du, respectively, is 


(1/27 .o)-"exp [— X(u, —)?/20] du,...du,...du,, (13) 
1 


since the probabilities of obtaining the individual velocities are independent. 
Written in terms of the easily calculable sample statistics U and s [equations 
(11 a) and (11 b)], this becomes 


(\/27.0)-" exp [ —n(U—)?/20?] exp [ —(n—1)s?/207] du,...du,...du,. (14) 


By a purely analytic process (10) the product of the differentials can be expressed 
in terms of the new variables U and s and n—2 other variables. After integration 
with respect to the latter, we find the probability of obtaining the sample is 


Cexp [—n(U —,)?/207] dU .s"-* exp [ —(n—1)s?/20*] ds, (15) 


where C is a constant. Hence the probability that the arithmetic mean of a 
sample of n line-velocities lies between U and U+dU is independent of the 
probability that the root mean square error of a single line-velocity, in a sample 
of n, lies between s ands+ds. Normalizing these two independent probabilities 
we find that the probability of the spectrum-velocity lying between U and U+ dU 
is 


P(U)dU=(V2a.0]ynyexp| —(U-w)9/= |av, (26) 
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and the probability that the root mean square error of a single line-velocity lies 
between s and s+ds is 


P,(s) ds = {2-921 [(n — 1)/2] }-* . [(n — 1) /o?]"-”, 5"? exp [ — (nm — 1)s?/20"] ds. (17) 

Integrating this over all s, we find that s the mean value of s is given by 
s=[2/(n—1)}1(n/2){T[(n— 1)/2]}*. 0, (18) 

so that for samples of n= 10, s=0-973¢ is a very good estimate of the unknown 


standarderroro. For moderate values of m equation (17) reduces to the Gaussian 
form 


P,(s) ds=[+/27 .0/4/2(n—1)}- exp | -«-a/ |«. (19) 


n—TI 


The problem now is to find a reliable value of s. For this purpose 106 
spectrum-velocities are selected from the 264 available observations so that 
they are fairly divided between the East and West limbs and the three days on 
which observations were made. The 106 values of s are then calculated according 
to formula (11b); the results are plotted as a frequency histogram in Fig. 5, 

















O10 O20 O30 O40 tans 


Fic. 5.—Distribution of sample standard errors. 


Ordinates: number of samples having a given error. 
Abscissae: root mean square error expressed in km sec. 


together with the theoretical distribution P,(s) (thin line, left-hand curve) from 
equation (17) for comparison. The mean value of s, s = 0-25 km sec (which agrees 
closely with the value 0-23 kmsec~! found for a slit position of Series II b and 
quoted in Paper I) gives for m = 10 (equation (18)) a value o=0-26km sec. The 
next problem is the reliability of the estimate of s, and hence of c, which is found 
from the dispersion of the values of s. ‘This is very easily calculated since the 
approximate Gaussian form of P,(s), given in equation (19) and shown as a thick 
line (right-hand curve) in Fig. 5, fits the observations as well as the exact form 
does. It gives for the standard error of o, simply o/\/2(n—1)=0-06kmsec"}, 
Thus we have as the best estimate of the root mean square error of measurement 
of a single line-velocity + (0-26 +0-06)km sec!. This estimate will be denoted 
by op. 

Consider now the distribution of the sample means. If the value of solar 
rotation were uniform and constant at all points in the equatorial region, the 
264 spectrum-velocities (plotted as a frequency histogram in Fig. 6) would 
themselves form a sampling distribution of the mean according to equation (16), 
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with a standard error o/\/10. Again o is the unknown root mean square error 
of a single line-velocity. The theoretical distribution P,(U) is shown as a smooth 
curve in Fig. 6. On this hypothesis o/\/10= +0-17km sec! or o= +0°54 
km sec!. Let this estimate of the value of o be denoted by oj. 
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Fic. 6.—Distribution of sample means. 
Ordinates: number of samples having a given mean. 
Abscissae: sample means expressed in km sec~. 

We have thus obtained two values for o, the standard error of a single line- 
velocity u. The value o,= + (0-26 + 0-06)km sec“! is derived from the sampling 
distribution of the standard error of a single line-velocity, calculated from the 
10 lines comprising a single spectrum-velocity. This estimate is independent 
of any variation of rotational velocity with time or position on the Sun. On the 
other hand, the estimate from the values of U gives o,= +0°54km sec“, but is 
valid only if the solar rotation is uniform and constant over the equatorial region, 
It differs from o, by nearly five times the root mean square error of the latter. 
The probability of the values o, and o, representing the same quantity in the 
infinite population is, therefore, less than one in a thousand, and we must reject 
the hypothesis of a uniform and constant rotational velocity throughout the 
equatorial region. That is to say, the large dispersion in the spectrum-velocities 
cannot be accounted for solely by errors of measurement and must be a con- 
sequence of actual variations in rotational velocity. 

4.2. Distribution of velocity residuals on the solar surface.—From this point 
onwards the investigation is purely exploratory, but it is interesting to see if there 
is any regularity discernible in the variation of velocity from point to point. 
Series III is obviously best suited to provide an answer to this question. In the 
solar rotation determination (Section 3) only 4 spectra from Series III are included 
from each limb, but in the present investigation the entire 10 spectra are included, 
giving a total of 120 measured velocities at each limb. 

The velocity residuals, U — U* cos B cos y, — A, for one slit position on the East 
limb are shown in Fig. 7. (The straight lines joining adjacent velocity points 
have no physical significance and are only intended to serve as a guide to the eye.) 
The figure shows no regular change in velocity over this portion of the solar 
surface. It may well be, however, that a regular variation of velocity is being 
obscured by the superposition of errors of measurement and randomly distributed 
local solar velocities. If a regular distribution exists, it can be revealed by taking 
the mean of a sufficient number of observations at any point. 





A. B. Hart 


£8 Ee eee 


oh. * 
t 


— Lins 
| 


t l 
00 





Fic. 7.—Typical distribution of velocity-residuals along an East-West chord of the solar image. 


Ordinates: residual velocity. 
Abscissae: heliographic longitude. 


In Series III the times of observation were such that 9 of the 10 exposures on 
each limb fall into 3 groups of 3, the time elapsing between the first and third 
exposure in each group being 40 seconds. ‘The time interval between the first 
and second groups is about 5 minutes and that between the second and third 
groups about 3} hours. (‘The seventh is the single spectrum remaining and it 
follows a few minutes after the second group.) The mean velocity at each point 
is calculated for each of the three 40-second groups. Figs. 8 and 9 show the 
results for the groups at the West and East limbs, respectively. Vertical lines 
show the root mean square error of the velocity (calculated from the agreement 
of the three individual velocities) at each point and, once again, adjacent velocity 
points are joined by straight guiding lines. 
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Fic. 8.—Distribution of mean residual velocities with longitude on West limb. 
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Ordinates: residual velocity. 
Abscissae: heliographic longitude. 


Attention may be drawn to the marked similarity of velocity variation 
between the three groups E,, E, and E, on the East limb. On the West limb 
there is little similarity between W, and W,, the groups separated by five minutes 
in time, though there is a similarity between W, and Ws, the groups separated 
by some three and a half hours. The steady variation in velocity at the East 
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limb over several degrees of longitude approximately equal to a distance of 
75000km along the solar circumference, together with the marked similarity 
in velocity pattern in these three groups, seem to point to the existence of a steady 
velocity field persisting over a period of some hours. The gradual dying-away 
of the velocity variations near the limb in all six cases presumably arises as a 
result of the averaging effect of measurement over a greater length of solar 
circumference. 
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Fic. 9.—Distribution of mean residual velocities with longitude on East limb. 


Ordinates: residual velocity. 
Abscissae: heliographic longitude. 


4.3. Possible systematic error.—At this point the objection may be made that, 
despite the precautions which have been taken, the velocity shifts are produced 
instrumentally. For example, it is conceivable, though not very probable, that 
a systematic velocity pattern exists in the standard spectrum (cf. Section 1) 
made from the out-of-focus image of the centre of the Sun; it is against this 
standard spectrum that the limb spectra are measured. Any such spurious 
velocity pattern is most easily found by taking the average of a large number of 
velocity residuals at each point on the slit, since for different areas on the Sun at 
different times the true residual velocities are randomly distributed about zero. 

Accordingly, the 264 velocities of Series II and III, included in the solar 
rotation determination and measured against the same standard spectrum, were 
divided into 15 groups, the members of each group corresponding to approximately 
the same point on the slit. The average residual velocity of each group was 
calculated in order to find the best available estimate of any correction for that 
region. Over the portion of the slit where a considerable number of independent 
solar velocities contribute to the mean of the group the correction does not exceed 
+0:03kmsec"4. This value should provide an upper limit to the size of any 
systematic error. Hence there is no evidence for supposing the correction differs 
significantly from zero and it is clear that the residual velocities cannot be attributed 
to instrumental causes. 

4.4. Evidence for velocity fields—Having confirmed the solar origin of the 
velocities, it only remains to find the probability of the similarity between W, 
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and W,, E,, E, and E, arising by chance. It can be decided by calculating the 
correlation coefficients r between the various groups. Correlations are calculated 
for the pairs W, W., W. W;, W; W, and for the corresponding groups on the 
East limb. 

In the case of E,, E,, E; the similarity in pattern is marked, but E,; appears 
to be shifted, relative to the other two, by about 1° in longitude. ‘This may 
be quite a spurious effect, since the longitude assigned to each point depends 
on the estimated position of the limb proper on the plate. The shift could arise 
from a difference of 0-I mm in estimating the position of the limb as between 
E, and E, on the one hand, and E, on the other. This difference lies within 
the possible error of éstimating the limb by visual means (cf. Section 1), and a 
change in seeing conditions, affecting the estimate of limb position, is not unlikely 
in a period of some three hours. The correlation coefficients are calculated, 
taking into account this shift. 

The average 7 of the 3 coefficients on the West limb and on the East limb are 
shown in Table V(a). The significant value of r, 74,,, for the relevant number 
of pairs of observations, at which there is a 10 per cent chance of the similarity 
between the two groups arising purely by chance, is obtained from tables (6). 
The value of 7,ign, given is the mean of the three values and is shown for com- 
parison. That the probability is only one in ten of the similarities being spurious 
in the case of the East limb may be seen from the table. 








TABLE V 
(a) (0) 

Group 7 Series Limb 7 
W,,W2,W; +0°32 Ila W. +0:26 
E,, E,;, E; +052 E. +0°34 

IIb W. +0°56 
Faign. +0°46 E. +0°47 
7 sign. +0:28 


Finally, we turn to Series II to see if any relation exists between the velocity 
variations along parallel chords of the disk slightly separated from one another. 
If the velocity field is steady over the period of some hours the 10 minutes elapsing 
between the first and last exposure on a limb (cf. Table I) is negligible and the six 
spectra may be regarded as a map of the region. Since only one exposure was 
made in each position, the observational evidence is very much weaker than in 
Series III and the results must be treated with reserve. 

Fig. 10 shows the velocity residuals along the West limb chords, marked 
in Fig. 2, for Series IIa. Similar variations are shown by the remaining obser- 
vations of Series II a and by Series II b on West and East limbs. There appear 
to be bands of velocity occurring approximately perpendicular to the solar 
equator. That this is not a completely spurious effect is checked by a method 
similar to that used for Series III. On any one limb, each of the six spectra is 
in turn used as the standard spectrum, and the value of the correlation of all 
remaining five velocities at each of the same longitude positions is calculated in 
this case. ‘The mean of the six coefficients on each limb is listed in Table V (6) 
in the column headed 7 and as before the value is shown at which there is a 10 per 
cent probability of the similarity being due to chance. ‘The mean observed 
correlation is well above the 10 per cent significance level, but this should perhaps 
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be regarded with some reserve since the single spectra being correlated have in 
them the full effect of errors of measurement and possible fluctuating solar 
velocities. 

An attempt was made to resolve the bands of constant velocity into their 
u, Vv, w components (see equation (1) in Section 3). The coefficients of u, v, w 
do not vary significantly, however, over the area of measurement, so that the 
equations are ill-conditioned and no solution for u, v, w could be found. 

Conclusion.—The general problem presented by measures of the equatorial 
solar velocity appears to be considerably clarified. The effect on the measured 
velocities, of light scattered in the Earth’s atmosphere and by the telescope, is 
shown to be only of the order of 0-03 kmsec"4. There remain variations in 
velocity of the order of + 0-17 km sec! which cannot be accounted for by errors 
of measurement. ‘These velocity differences are of the same order of magnitude 











Fic. 10.—Velocity residuals as a function of longitude and position on disk. 


Ordinates: residual velocity. 
Abscissae: heliographic longitude. 


as those originally found by H. H. Plaskett. The present observations go 
further, however, in showing that these differences may arise from velocity fields 
in the form of bands of constant velocity separated by distances of the order 
of 75000km. These bands possibly extend perpendicular to the equator over 
distances of not less than 60 000 km. 

If such velocity fields exist they might go far towards accounting for the 
results obtained by earlier workers, shown in Fig. 1. The dispersion of the most 
reliable values is + 0°05 km sec’! and this is the order of the variation in the mean 
values found on single days in the present investigation. Such large differences 
are to be expected for the mea: of a few observations made at a single place in a 
velocity field. 
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While the presence of velocity fields in the equatorial photosphere is suggested 
by the present investigation and accounts for the earlier discrepant values of the 
equatorial rotation, the existence of such velocity fields still requires demonstration 
from more extensive observations. Some of these observations are now in progress 
at Oxford. 
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APPENDIX 


In the calculation of the velocities to be expected in the aureole and inside 
the limb, the arrangement and extent of the scattering and solar zones are the same 
as in Section 1 of Paper III. B,, is the average brightness of the mth solar zone 
from the limb; U, is the average value of the scattering function in the pth 
scattering zone, the outer radius of the eighth being s= 120” arc, and beyond this 
the total light scattered is approximated by a small constant quantity, 0. Let 
0», be the area common to the pth scattering zone and the mth solar zone when 
the point X, to which scattering is taking place, is a distance x from the limb. 
Then the total intensity of light scattered to X is given by 


I(x) ai {2(Zo,,U,)B,} + Q, (A I) 


Paper III, Section 1, equation (2). 

Each solar zone scatters light of its characteristic velocity; for the mth solar 
zone let this be V,. It is calculated from the equatorial velocity of rotation, 
taken as 2-07km sec, by assuming that the centre showed zero line-of-sight 
velocity while the limb showed the full value. Hence the product of brightness 
and velocity for the area of the image included in the scattering zones is 


X(Zo,,U,)B, Ve: 
np 


For light scattered from outside the area included in the scattering zones the mean 
velocity is taken to be zero, since the velocity distribution is symmetrical about 
the solar axis. This approximation improves as X moves further from the 
limb along a solar radius perpendicular to the axis. 

Hence, for points outside the limb in the equatorial region, the observed 
velocity V*(x) due to scattered light is given by 


X(Zo,,,U,)B,V, 
np 
2@e,.0,)B, +O" 





V*(x)= (Az) 
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(The validity of this step is discussed later.) The calculations were performed 
for x=7"-5, 15” and 25” arc along a solar radius, corresponding to distances 
o-8 mm, I-6mm and 2-7 mm from the limb along the E.—W. line of the telescope. 
The results are given in Table AI (a), and yield the curve given in Fig. 3. 

For points inside the limb, however, the light of intensity J(x) and velocity 
V*(x) is superposed on that remaining from the true intensity J after scattering 
has taken place. This residual intensity is approximately 90 per cent and in 
this paper a factor 0-905J has been used. The value of J, to a sufficient accuracy, 
is found from an observed mean limb-darkening curve. If V is the true velocity 
at the point considered, the observed velocity is 
0-905] x V +1(x)V*(x) 

ooosl+ M(x)’ (A3) 
and the values of V(x) at points 0”-o and 10”:0 arc inside the limb are given in 
Table A1I(4). These lead to the results discussed in Section 2. 


TasLe Al 
(a) (6) 


Radial distance Calculated Radial distance Calculated Velocity without 
outside the limb velocity inside the limb velocity scattering 
y ih 1°25 km sec™? 0”'o 2°03 km sec 2°07 km sec? 
15”-0 0°93 10”-o 2°02 2°05 
25”°0 0°80 





V(x)= 


The only point remaining is to justify the assumption in equations (A 2) 
and (A 3), that a change in intensity produces a linearly proportional change in 
velocity. This was done by measuring the effect of blending a spectrum of one 


intensity and velocity with another of differing intensity and velocity. On a 
narrow (3 mm) spectrum from a point 20” arc inside the limb a spectrum of the 
centre of the disk was superposed directly with the same exposure time, the dark 
slide of the plate holder and the plate position remaining untouched between the 
exposures. The intensity of the centre spectrum was reduced to a known 
fraction of its original intensity by N.P.L. calibrated filters. Each plate also 
carried an unblended centre spectrum and two unblended limb spectra, taken 
at the same position as the spectra used for blending, to give the velocity difference 
between the unblended and superposed spectra. Five plates were taken at the 
West equatorial limb and four at the East equatorial limb. 
Taste A II 
Fraction of : Cohaiaet 
eunnndl tnalte Io/Iy, (Vp—Vp)/Vi—Ve)  =(Vi—V)/(Vi—Vo) 
per cent per cent per cent 
O'O1I7 2 os = $3 
0184 4 o'4 3°0 
-0382 7 6-0 2°3 
0597 II 6°8 $°3 


Column 1 of Table A II shows the fractions of the centre intensity transmitted 
by four combinations of filters. From a mean limb-darkening curve, the ratio 
of the intensity (at 6000 A) of the centre of the disk to that at 20” arc from the 
limb is 1-9 so that the values of J,/J;, given in the second column of the table 
represent the ratio of the superposed centre intensity to the limb intensity. The 
whole height of each spectrum was measured on the Evershed machine and the 
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third column shows the values of V;,— Vx expressed as a percentage of Vi,—V, 
together with the root mean square error of each measurement. V3, Vo, Vy 
are the measured velocities of the blend, centre and mean of the two limb velocities, 
respectively. A least-squares solution was made for a linear relation between 
these four pairs of values and the point J,/J,=0, (Vy—Vx)/((Vit—Vc)=0. The 
values obtained from the solution are shown in the last column of the table. 
For values up to J,/J;,=10 per cent the representation is obviously adequate. 
We may, therefore, safely assume that addition to a spectrum, of a small quantity 
of light of differing velocity results in a fractional change in velocity varying 
linearly as the fractional change in intensity. 
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I. SHowER METEORS 
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Summary 


Apart from a slight influence of the directional properties of the aerial 
system, the form of the radio-echo height distribution for a homogeneous 
velocity group of meteors constituting a stream with a well-defined radiant 
point is found to be a function only of the atmospheric scale height and the 
mass distribution of the incident meteors. The relation between these 
parameters and the width of the height distribution (expressed as the r.m.s. 
deviation from the mean) is derived, hence if one is known the other may be 
determined from the observed meteor heights. ‘The mean height corresponds 
to a value of the atmospheric pressure which is related to the physical properties 
of the meteors, their ionizing efficiency and velocity, and the radio-echo 
equipment parameters. The theory enables the observed dependence of the 
mean height on the velocity to be interpreted in terms of the variation of 
atmospheric pressure’ with height. 





1. Introduction.—The passage of a meteor through the atmosphere produces a 
column of ionization which may be detected by the scattering of radio waves from 
it. In the absence of distortion of the linear trail (by winds, etc.) the reflection is 
specular, i.e. the important part of the trail is a length of the order of 1/(RA) 
(R=range, A= wave-length) centred about the normal intersection of a radius 
vector from the observing station with the trail. We will call this the reflection 
point. ‘The characteristics of the observed radio echo are dependent primarily 
upon the electron line density in the trail and on the properties of the atmosphere in 
the vicinity of the reflection point (1, 2). 

Clegg and Davidson (3) have described an apparatus for determining the 
elevation and range, and hence the height of the reflection points. The purpose 
of this paper is to develop a theory of the dependence of the observed height 
distribution on the properties of the incident meteors and of the atmosphere in the 
meteor ionization region, and on the parameters of the radio equipment. The 
discussion is restricted here to a homogeneous velocity group of meteors constitut- 
ing a stream with a well-defined radiant point. 

2. Theory of meteor ionization.—In passing through the upper atmosphere, a 
meteor is subject to individual collisions with air molecules, a great proportion of 
which are trapped in its surface. The heat liberated causes the meteor to 
evaporate atoms which leave with thermal velocities relative to the meteor surface 
but carry the meteor velocity relative to the surrounding atmosphere. Subsequent 
collisions between meteor atoms and air molecules produce ionization and 
excitation, mainly of the meteor atoms (4). 
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Herlofson (5) and Whipple (6) have summarized the modern theory of the 
motion of a meteoric particle in a rarefied atmosphere, from which we obtain the 
following expression for the rate of evaporation : 


n= (wH)* mv cos x (P/Pmax)[1 — $(P/Pmax)]” 5 (1) 
Pmax is the atmospheric pressure at the point of maximum rate of evaporation and is 
given by 
Pmax = [(2/g)/(Av*A)]}m*® cos x, (2) 
where A = .Am-*, of is the initial cross-sectional area of the meteor in its line of 
flight in cm* and m is its initial mass in g. It should be noted that A, which 
Whipple calls the ‘‘ shape factor ’’, is independent of m, being a function only of 
the shape of the meteor and of its density. ‘The other quantities appearing in (1) 
and (2) are as follows: 


n (sec!) = number of meteor atoms evaporated per second, 
v(cm sec~*) = meteor velocity, 
(g)= mass of an individual meteor atom, 
l(erg g-!) =latent heat of evaporation of the meteor, 
x = zenith angle of the meteor radiant, 
p(dyne cm-*) = atmospheric pressure, 
H (cm) =atmospheric scale height, 
g (cm sec~*) = gravitational acceleration, 
A = heat transfer coefficient. 


The heat transfer coefficient (6) is defined as the ratio (heat energy liberated)/(air 
mass intercepted by the meteor x $v”). Herlofson takes A=1, i.e. he assumes 
that all impinging air molecules are trapped in the meteor surface, while Whipple 
suggests a value of about 0-7. 

In deriving (1) and (2) the deceleration of the meteor is assumed to be small, 
and approximations have been made which require v?>12/. The error at small 
velocities mainly affects equation (2). The zenith angle is assumed to be constant 
along the length of the trail, i.e. (go° — x) must be at least several degrees. 

The electron line density produced by the meteor is «= fn/v, where f is the 
probability that a single evaporated meteor atom will give rise to a free electron; 
B may be a function of v. 

We will restrict ourselves to the case of a meteor shower, comprising a homo- 
geneous velocity group of meteors with a well-defined radiant. Writing 


ky=B\(uH), Ry=2lg/(Av*A) 


we find that « electrons per cm are produced, at pressure p, by a meteor of mass 


_[ (ke2\22 p 3 
m=| (7) ~ato | ; (3) 


3. The meteor mass distribution.—Since pmax is a function of the meteor mass, 
the height distribution will depend upon the mass distribution of the incident 
meteors. 

From a statistical analysis of visual observations Watson (7, 8) and Williams (9) 
have shown that, over a wide range of visual magnitudes, the frequency increases by 
a constant ratio, a, per magnitude. This leads to the mass frequency law 


vy =bm~, (4) 
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where vydm is the number of meteors in the mass group m—(m-+dm), crossing 
unit area normal to the stream per second, 4 is a constant and 


s=1+2°5 logiya. (5) 
The above authors give values of a between 2 and 4, hence we may expect s to be 
in the vicinity of 2. 

4. Radio-echo height distribution.—In the main part of this discussion we shall 
assume that the sensitivity of the apparatus is such that the majority of the meteor 
trails detected have electron line densities less than 10! cm™ (short-duration 
echoes) and take the echo amplitude as directly proportional tox. The case where 
all echoes are from trails of line density greater than 1012 cm will be dealt with in 
Section 4.5 below.* 

On account of the specular reflection from the trails, all of the echoing points 
lie in a plane ABCD (Fig. 1) passing through the observing station, O, normal to 
the radiant direction. ‘The constant height (and pressure) contours are slightly 
curved lines, SS’, in this plane. 











A 


Fic. 1.—The echo plane ABCD is normal to the radiant direction and is at elevation x above the 
horizontal, m/2—8 is the angle in the echo plane, measured from the horizontal. 


1.1. Height distribution from echoes obtained in the sector 0 — (0 + d0).—Consider 
the point P, at height A, in the sector 6— (6+ d6) of the echo plane (Fig. 1). The 
minimum electron line density detectable at P is given by 


up = KR*? S-%(6), (6) 
subject to ~p<10!*, .S(@) represents the aerial polar diagram (amplitude) in the 


echo plane, R(cm)=OP, and K is a constant determined by the equipment 
parameters. It is given (1) by 
K=6:3 x 104 (9 G,,)-* [e/(PA®)]¥*, (7) 
where G,,, = aerial power gain (referred to an isotropic radiator) in the direction of 
the maximum of the beam in the echo plane, 
A(cm) = wave-length, 
P= peak trahsmitter power, 
€= minimum detectable echo power, 
7 is a constant which depends on the orientation of the incident electric 
field vector to the meteor trail. Its value lies between 1-0 for parallel 
polarization and 2-0 for transverse polarization. 


* For a discussion of the variation of echo amplitude and duration with electron line density, 
see Kaiser and Closs (1); Closs, Clegg and Kaiser (2); and Greenhow (10). 
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From (3), writing 
2 = (p/ks)[9up cos? x/(4ky)}**, (8) 

we find that the smallest detectable meteor passing through P has a mass 
ap 


tps ——— 
12k, cos x 


(22-8? + 2). (9) 


The number of detectable meteors passing per second through unit area at P is 
obtained by integrating the right-hand side of (4) with respect to m, between limits 
mp and infinity. ‘Thus the rate at which meteors will be observed in the sector 
6—(6+ dé), and within the height range h — (h + dh), is v, d0 dh, where 
oe 12k, cos x\*1 R I 
i Gab. oe (RIR,+cosOsinx)(@at=+apen = (7°) 
R, is the Earth’s radius. 

Except for «p in (8), z varies directly as the pressure. When z=1, let ap=«, 
p=p, R=R,, h=h,. Then 

== (P/P1) (%1/%p)"* =(0/P1) (Ri/R)*”. 
Now R cannot vary more rapidly than directly as h, so putting R,/R=h,/h we get 


r= (rH in?) 
Pi 2h, py 





where 

Pi= kz [9% cos* x/(4k,)]**. (11) 
Since meteor echoes are observed only over a few scale heights, with a mean 
height in the vicinity of 100 km, we can put z=p/p,. For the same reason we find 
that the only significant height-dependent term in (10) is that containing z and, 
although R is a function of height, negligible error is involved in putting R= R, 
throughout. The normalized height distribution for echoes observed in the 

sector 6 — (6+ d6) is thus 

3(s—1) 

sts) (12) 





b (“ cos x“ R, 


1 o~ * 
— (s—1) (R,/R,+ cos 8 sin yx) ach (13) 
subject to the restriction that y does not approach 90°. 
In most cases cos @ sin x > R,/R,and 


b (4k, cos x\*1 h, 
r (s—1) ( Qa, ) cos? @ sin y sie (14) 


4.2. Properties of the height distribution v,*.—Putting x, = —1In z we have 


3e%1 3(s—1) 
o Cae F eee 
= \3e8r24 7 ? (15) 


where, for constant scale height, 


Qa, 





x= (h—h,)/H. (16) 


Fig. 2 gives v,* as a function of x, for a range of values of s. ‘The most probable 
height occurs at x=0, i.e. at height h,, pressure p,, and is independent of s. 
It is easily shown that «,, the minimum line density detectable at range 
R,=h,/(cos @ sin x), is that produced by the meteor which has a maximum rate of 
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evaporation at pressure p,. The width of the height distribution is a function only 
of s and of the scale height. 


Due to the slight asymmetry of the height distribution, the mean value of x, 
is not zero. Denoting it by Z, we have 


| xv, * dx, 
€,=(h—-h,)/H= T we . (17) 
vy * dx, 


The r.m.s. deviation of x, from the mean value we will denote by 5x,. Thus 


‘¥ (x,—2,)?», *dx, 
(5x,)?= —- . (18) 
3 = v4 * dx, 


Both <, and 5x, are independent of 0, being functions only of s. They have been 
determined by numerical integration and are given in Fig. 3. 

The integral of v, * with respect to x, will be of interest for later work concerned 
with the total meteor echo rate. Toa good approximation it is given by 





vy* dx, = 2°3 (s—1-13)-™, (19) 
which is accurate to better than 1 per cent for s>}. 

4.3. Complete height distribution, polar diagram symmetrical about 6=0.—The 
complete height distribution, v, is obtained by integrating v, with respect to 6, 
remembering that v,*, through x, and f,, is a function of 6. To avoid confusion, 
we will let the values of h,, «1, p, for @=0° be denoted by h,,, %&,; Pm, and put 
x=(h—h,,)/H. ‘Thus they are related by the same equations (11) and (16). 

Anticipating that the normalized distributions will not be strongly dependent 
upon the exact form of the polar diagram and on the other 6-dependent terms 
multiplying v,* in (13) and (14), we will restrict ourselves to consideration of (14) 
and neglect the variation of h, with @. Thus we can put 


a, = a,, [cos®? §.S2(9)}-1 = a,, f-*(u), (20) 
where u=tan 06, f(u)=cos** 0S (0) and «,, is the minimum detectable line density 
at range h,,/cos x in the direction of the principal maximum of the aerial beam (in 
the echo plane). Hence vis proportional to 


c 


| © PN (u) v4 * de 


Now = —In(p/p,)= —In(p/pp»)—1n (Pp/P1) =*— $10 (p/). 
Therefore x, =x— $#lnf(u), 


and is e* %(e-1) 
“= |" Lame) 


Mean height 
The mean height is given by 


| Jxreu)n* du dx, 





#=(h—hp)/H= 


| | f*(u)vy* du dey 
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Fic. 2.—Normalized distribution in height, v,*, of echoes observed in the sector @—(0+-d0) of the 
echo plane. x,=(h—h,)/H. 
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Fic. 3.—Parameters of the distribution v,* as a function of s. 
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the integrals being taken for both u and x between — «© and + 00. Substituting 
(21) in (23) and remembering that 


2) «eo 
| %,¥,*dx, and | v, *dx, 
—- eo —-@ 


are independent of #, we have 


[few In fu) 
{ f2-(u) du 





E—Z,=§% 


R.m.s. height deviation from the mean 


Ther.m.s. deviation of x from £, denoted by 5x, is obtained from the elementary 
height distribution, v,*, through 


[Je — £)f**\(u)v,* du dx, 


a | | f-(u)v,* du dx, 





Using the previous expansion (21) for x,, we find 
(x— &)? =(x, — 4)? — 2(x, — %,)(F— %,) + (Z— %,)? + $(x, — Z) Inf (u)+§ ln? f(u), 
leading to 
[Pema int fu) de 
(5x)? = (5x1)? —(#—#,)? + §—z 
| f?®%\(u) du 





(25) 


The value of (24) and (25) is that they enable the mean height and the r.m.s. 
height deviation to be evaluated in terms of the corresponding quantities for the 
elementary distribution v,*, eliminating the necessity for evaluating v* explicitly. 

Simple physical considerations show that the main effect’of the polar diagram on 
the height distribution derives from the behaviour of f(u) in the vicinity of u=o. 
For this reason, the parameters of the height distribution have been determined for 
the following forms of f(u): 


(a) (1+c*u?)-! (b) exp(—c*u?) (c) cos*cu for —2/2<cu<n/2 
and zero otherwise. 


All vary as (1 —c®u?) for small u but approach zero at different rates as u—> 00. 
# and dx have also been evaluated numerically for S(@)=cos(7/2 sin 6), 
which corresponds approximately to the case of a horizontal half-wave dipole 
spaced <A/4 above a reflecting ground plane. Almost any actual polar 
diagram will lead to a function f(u) between the limits of (a)and(c). For example, 
the polar diagrams derived from (a), (b) and (c) for c?= 1-61 are given in Fig. 4(a), 
with cos (7/2 sin @) (which has the same form near 6 = 0) for comparison (curve(d)). 
The case of c?>1, when we can put 6=sin 6= tan @, is given in Fig. 4 () together 
with the polar diagram [sin (c/66)]/(c./6@) for comparison (curve (d)). 
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Fic. 4 (a).—Aerial polar diagram, S(@), in the echo plane. 
(a) f(u)=(1+1°61u*)™, (6) f(u)=exp(— 1-610), (c) f(u)=cos? (1-27u), (d) S(0)=cos (m/2 sin 8). 


LO; 











Fic. 4 (b).—Aerial polar diagram, S(0), in the echo plane for a beam sufficiently narrow that 
6 = sin 6 = tan 0. 


(a) flu)=(1+c*u*)-", (6) flu)=exp (—c*u*), (c) f(u)=cos* (cu), (d) S(0)=[sin (c+/68)}/(c/ 68). 
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Fic. 5.—Height distribution, v*, as a function of x for an aerial polar diagram (in the echo plane) 
represented by f(u)=(1+c*u*®)—.  x=(h—h»)/H. 
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Case(a) f(u)=(1+cu*?) 


The integrals which we require for the evaluation of Z and 5x become, using the 
substitutions cu=tanv, g=4s—6, 


ea 2 7/2 
| fu) du = = | cos! v dv, (26) 
“—@ Clo 
a on/2 
[ fu) Inf (u) du= 3 | cos‘ v ln (cos v) dv, (27) 
J oe J0 


= 8 22 
| f**M(u) In? f (u) du = . | cos? v ln? (cos v) dv. (28) 
J/—@ “0 


For integral values of g, (26) is easily evaluated, and for (27) and (28) we obtain the 
recurrence formulae: 
“1/2 m/2 rxf2 
q| cos*vln(cos wv) dv=(q—1) cos** v In (cos v) dv + | cos** vdv 
“0 “0 0 
7/2 
- | cos‘ v dv, (29) 


“0 


7/2 [2 
q { cos‘ vIn? (cos v) dv =(q—1) a cos** y In* (cos v) dv 
0 0 


7/2 n/2 
+2 | cos*? vIn (cos v) dv — 2{ cos? vIn (cos v) dv. (30) 
0 0 
These integrals have been evaluated numerically for g=0, I. 
From (22) the height frequency distribution, normalized to unity at x=0, is 
found to be 


is Zet 3-772 
v* (x) =e" Ls | , (31) 
which is plotted in Fig. 5 for various values of s. 


Case(b) f(u)=exp(—c*u*) 
Putting cu=v, we have 


[Mw du= = | e%-0¥'do, 
-—@ 0 


| a P2(u) In f (uw) du= — aa | ; eX) # do, 


3 _ feu) In? f(u) du= sy [ e%0* ao, 
I 
6(s—1)’ 
I 


(dx)? = (8x,)? + 186-1)" 


2 
0 


The function v* has not been determined for this case. 
Case(c) f(u)=cos*cu 


Putting cu=v, g=4(s—1I), we obtain the same equations as for Case (a), 
nameiy (26), (27), (28). 
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Equation (22) is integrable for integral values of 3s—7/2, giving 
, so— 7/8 $:4:6:.55% me 
ya(a) oc et-02 sixty —yy(r—y9) ES (r+ SO yee) | 5g 


where y = 2e°*? + 1, which is given in Fig. 6 for various values of s. 
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Fic. 6.—Height distribution, v*, as a function of x for an aerial polar diagram (in the echo plane) 
represented by f(u)=cos* cu. x=(h—h,»)/H. 
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Figs. 7(a) and 7 (5) show & and 5x as a function of s for the cases considered 
above. It will be noted that, as anticipated, the height distribution is remarkably 
insensitive to the exact form of f(u), and hence of the polar diagram. For the 
same reason, the complete height distribution for v, given by (13) will not be signi- 
ficantly different from those derived above. 

4.4. Complete height distribution, polar diagram symmetrical about §=06,.—In 
this section we will let h,,,, ~,,, etc., be the values of hy, «,, etc., at 0= 4p, i.e. as before, 
the quantities appropriate to the direction of the principal maximum of the 
aerial beam in the echo plane. Since the foregoing has shown that the height 
distribution is not critically dependent on the exact form of the function of 
multiplying v,* in (13) and (14), we will restrict the discussion to the case of an aerial 
beam sufficiently narrow that we may write «,—~z,, S~*(0’), cos @=cos 6), where 


6’=0—6,. Hence, from (14), 
b (4k,cosx\"1 hy, P 
<5 ( Qu, ) sin? 62” (35) 


where sin ¢=cos 6, sin x, i.e. ¢ is the elevation of the principal maximum of the 
beam intheecho plane. Putting tan 0’=u, f(u)—S(0’), v* is given by (22), i.e. is 
identical with the previous case, and the characteristics of the height distributions 
for the various f(u) are given by Figs. 5, 6 and 7. 

As the elevation ¢ tends to zero, v, may be found from (13), which, provided 0, 
does not approach 90°, leads to a result identical with (35) except that the factor 
h,/sin® ¢ is replaced by R,. 

4.5. Equipment sensitivity such that «,,> 10% cm~!.—We will deal briefly with 
this case for which only the long-duration echoes are observed. 

We proceed as before, except that instead of (6) we find «,,"* to be proportional 
to R®2,S-*(6) (1). Thus «,=«,,[cos*@S*(@)}-*. Equations (13) and (14) remain 





y= 
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Fics. 7 (a) and (b).—Parameters of the height distribution as a function of s for an aerial polar 
diagram (in the echo plane) represented by: 


(a) f(w)=(+cu?), (6) f(u)=exp (—c*u*), (c) f(u)=cos* cu, (d) S(0)=cos[(z/2) sin 6]. 
The broken curve in Fig. 7 (b) gives the r.m.s deviation, 5x,, for the elementary distribution v,*. 
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valid, as does the remainder of the analysis except that f4(u) replaces f(u) through- 
out. For f(u)=(1+c*u?)-! and f(u)=cos* cu both and 5x are more nearly 
equal to the values obtained for f(u)=exp(—c?u*), which remain unchanged. 
Fig. 8 gives & vs. s for the three cases ; the values of 5x are not given, since they do 
not differ significantly. 

For «,, in the transition region (i.e. echoes distributed more or less evenly 
between the two types) simple physical considerations show that the previous 
results are still satisfactory, with the values of # and 5x slightly nearer to Z, and 5x,. 


0.1 (c) 
x (b) 








(a) 


5 


at | J J 
1.0 1.5 2.0 2.5 
Fic 8.—x=(h— hm)/H as a function of s for an equipment of such low sensitivity that 
&m > 10!? cm-. 
(a) f(u)=(1 + c?u?)-1, (6) f(u)=exp (—c*u*), (c) f(u)= cos* cu. 








5. Diseussion.—The above theory makes possible the interpretation of the 
meteor height distribution observed on a radio-echo equipment during the 
activity of a major meteor stream. 

The width of the height distribution, expressed as the r.m.s. height deviation 
from the mean, depends on the atmospheric scale height, the mass distribution 
of the incident meteors, and only very slightly on the form of the aerial polar 
diagram. It will be shown in a later paper that the exponent s in the mass law may 
be obtained with fair precision from the echo amplitude and duration distributions ; 
thus the scale height may be determined. 

The mean height is close to the height h,, at which the pressure, p,,, is given by 

Pm = Re [9% COS? x/4R,]**. (36) 
It depends strongly on the meteor velocity which enters (36) through k, (which is 
proportional to v~*) and possibly through the probability of ionization, 8, con- 
tained in k,. It is also dependent upon «,,, the minimum electron line density 
detectable at range h,,,/sin ¢ in the direction of the principal maximum of the aerial 
beam in the echo plane; «,, may be evaluated from the equipment parameters or 
derived from the amplitudes and durations of the long-duration echoes. Thus 
Pm can be obtained in terms of the physical properties of the meteors (latent heat, 
density, velocity), the radiant zenith angle y, the heat transfer coefficient and the 
probability of ionization. The latter need not be known with great accuracy, since 
it enters (36) as the cube root. 

The zenith angle y has been treated throughout as constant. It is thus 
desirable that the height distribution be measured over a selected time interval 
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during the passage of the radiant across the sky, sufficiently short that the variation 
of x may be neglected. This is not a severe restriction, since p,, is proportional to 
cos**y, As an example, if x varies from 35° to 55°, the mean height changes by 
about 0-25 scale heights and if observations were extended over this range of x, 
the width of the height distribution for s~2, would not be more than some three 
per cent greater than the theoretical values derived above. ‘The assumption of a 
well-defined radiant point merely requires the spread in meteor radiants to be small 
compared with the aerial beam width. 

We have represented the polar diagram of the aerial beam in the echo plane by 
S(@). In the case of an equipment using separate transmitting and receiving 
aerials with polar diagrams S,(@), S,(@) and gains G, and G, respectively, we 
simply write S?(@)=S,(0) S,(0), G?= G, G,. 


University of Manchester, 
Jodrell Bank Experimental Station, 
Lower Withington, 
Cheshire: 
1953 July. 
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THEORY OF THE METEOR HEIGHT DISTRIBUTION OBTAINED 
FROM RADIO-ECHO OBSERVATIONS 


II. Sporapic METEORS 


T. R. Kaiser 
(Communicated by A. C. B. Lovell) 


(Received 1953 August 21) 


Summary 


The theory of the meteor height distribution is extended to the case of a 
selected velocity group of sporadic meteors. 

The width of the height distribution (r.m.s. deviation from the mean) is 
determined as a function of the atmospheric scale height and the mass distri- 
bution of the incident meteors. This latter is shown to be almost exactly 
inverse square. For a given velocity the mean height is shown to occur at a 
value of the atmospheric pressure which is related to the physical properties 
of the meteors and of the ionization process and to the radio-echo equipment 
parameters. 

The measurement of meteor heights by the radio-echo technique is 
suggested as a new and important method for determination of the atmo- 
spheric pressure and scale height in the meteor ionization region. 





1. Introduction.—In Part I (1) we have derived the height distribution of radio 
echoes from the trails formed by meteors belonging to a homogeneous velocity 
group, constituting a stream with a well-defined radiant point. In this paper 
the theory is extended to the case of a selected velocity group of meteors whose 
radiants are distributed over the celestial sphere. On account of the specular 
reflection property of the trails, only those meteors will be observed whose radiant 
points correspond to echo planes intersecting the aerial beam. Thus, referring 
to Fig. 1, the radiant points of meteors observed in the conical element of solid 
angle of half angle % must lie within the shaded strip of the celestial sphere of 
width 2y. 

We will make the assumption that the geocentric radiants are uniformly 
distributed over the celestial sphere. At any instant, the radiants for a given 
velocity group will in fact tend to be concentrated about the apex of the Earth’s 
way to an extent dependent upon the heliocentric radiant and velocity distri- 
butions. However, if the meteors are observed over an extended period of time 
from a station in northern latitudes with an aerial beam directed northwards, the 
instantaneous radiant distribution will sweep across the shaded region of Fig. 1 
and the effect will be similar to the case of a uniform distribution. In any case, the 
relatively small differences between the results obtained in the following for 
sporadic meteors and those given in Part I for shower meteors show that the 
height distribution is not very sensitive to the exact form of the radiant distribution. 
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On the above assumptions, the flux of meteors per second crossing unit 
area normal to their direction of motion, with velocities v to (v+ dv), masses m to 
(m+ dm) and with radiants within a solid angle dQ is vy dv dm dQ, where, as in 
Part I, we may write 

vy =5(v)m—. (1) 


The discussion in this paper is restricted to the case of a radio-echo equipment 
with sensitivity such that echoes from trails with electron line densities less than 
10% cm~! may be detected over most of the solid angle covered by the aerial beam. 
Thus the minimum detectable electron line density at a point P, distant R cm from 
the observer, is (2) 

ap = KR®?S-(¢4, 8), (2) 


where ¢, # are respectively the elevation and azimuth of P relative to the observer, 
the zero of # being chosen as the optimum azimuthal direction of the aerial beam. 











Fic. 1.—Meteors observed in the conical element of the aerial beam with semi-angle & have 
radiant points Q (coordinates (€, £)) within the shaded strip of the celestial sphere. 


S(¢, #) is the aerial polar diagram (amplitude) and is normalized to unity at (¢,,, 0), 
K is as defined in (7) of Part I, the gain G,, referring to the direction (¢,,, 0). 
The elevation, ¢,,, will be taken as the maximum of the distribution of echo 
elevations. It is not, in general, identical with the elevation of the aerial beam 
maximum (especially for a vertically directed beam), hence on the above 
definition the maximum value of S(¢, #) may be greater than unity. 

2. Meteors observed in a cone of solid angle dw at elevation ¢.—We will first 
consider those meteors, radiating from the element dé d{ at Q on the celestial 
sphere, which are observed in the narrow cone of half angle #% at elevation ¢ 


(Fig. 1). 





54 T. R. Kaiser, Theory of the meteor height Vol. 114 


The rate of echoes from meteors with velocities v to (v+dv) observed in a 
sector 6 to (@+ d@) of the echo plane corresponding to Q (see Fig. 1, Part I) and in 
the height range / to (h+ dh) is v, dv dh dé dé dl, where v, is given by equation (13) 
of PartI. Now the cone intercepts an angle 2(,;* — ¢?)" of the echo plane, hence, 
by putting cos y=cos¢ cos , where x is the zenith angle of Q, m/?= dw, in (13) 
of Part I and integrating with respect to é and ¢, the rate of echoes observed in 
an element of solid angle dw, having coordinates (¢,#), in the height range 
h to (h+dh), is vgdvdhdw, where, in the notation of Part I, 


— (v) (4k, cosd\*1 R, ae 
7. (s—1)\ ga, (K,/R,.+sin¢). 


3e*s 3(s—1) 
Vv i °*(x,)= Eas | > (4) 


a = KR,**S~*(4, 3), : (5 
v, = —In(p/p,) = (4—A,)/A, (6) 
1 = ke 9x, cos? $ cos? £/4k,)*%. (7) 
The height /,, at which the pressure is p,, is the most probable height at which 
meteors from (€, ) will be observed in the element dw. Thus A, is a function of & 
although, as in Part I, wherever it appears only as a multiplying factor (rather 
than an exponent) it may be treated as constant and taken as approximately 
equal to the mean observed height, A. R, is the range to height A, in the direction 
(¢, #), and except for the case of a narrow aerial beam with elevation [10° we 
can put sin ¢>R,/R,, R,=A,/sin ¢. 
Let hg, %, Xa, Py be the values of h,, «1, x,,p,at€=0. Thus, from (5), (6) and (7), 
Pi=p2cos** €, x,=x,+ §ln(cosé), x,=(h—h,)/H 





cos*! £y, * dé, (3) 


—a/2 


and y, is a function of x,. ‘To a sufficient approximation we can set «, =a». 

Following an analysis similar to that in Section 4.3 of Part I, we obtain the 
following expressions for x, and 5x2, the mean and the r.m.s. deviation from the 
mean respectively : 


-n/2 
| cos*€ In (cos ¢) dé 
2 0 





-n/2 ’ (8) 
3 | cos*é dé 
0 


4 4 cos? € In? (cos €) dé 
(5x2)? — (8x1)? = —(%_— %,)? + 9 rr) , (9) 
| cos*é dé 
0 


where g=s—1. The quantities Z, and dx,, which refer to the function v,*(x,), 
have been evaluated in connection with the shower meteor height distribution 
(see Part I), as have the integrals in (8) and (9), for a range of values of q. 

3. Complete sporadic meteor height distribution.—From (3), making the approxi- 
mation sin > R,/R,, we may write 


V9(X_) < (cos p/x_)*? sin-* dv_*(xp), (10) 


where v,* is normalized to unity atx,=0. The complete height distribution, v, 
is obtained by integrating (10) with respect to ¢ and #, remembering that a, 
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is a function of (¢,%). The aerial polar diagram function S(¢, #) can in most 
cases be adequately represented as a product S,(¢)S,(#) of the polar diagram 
functions in the principal planes of the beam. Thus S,(0)=S,(¢,,) =I. 

Let the values of hy, x, py, in the direction (¢,,,0) be A,,, %p» Py», and put 
x=(h—h,,)/H. Now a, is the minimum line density detectable at range 
h,/sin¢d in the direction (¢,#), thus, neglecting the small dependence of h, on 
(¢, #), we obtain, using (5), 


Hg = Op (Sin $,,/sin $)??S,-*($)S_-*(9). (11) 


From (10) and (11), with dw = cos ¢ d¢ dd, we obtain, for the complete height 
distribution, 


a)ac f |” CA(G)S (OP >gay) apd, (12) 


0 


where 





F114) = ( cos ¢ y ( sind ie S114). (13) 


COS¢,/ \SINd» 


The function F**-)(¢) is proportional to the number of echoes observed per unit 
elevation interval; thus ¢,, is the elevation at which F(¢) has its maximum 
value, which we have normalized to unity. It is clear from (13) that ¢,, is not in 
general the elevation of the aerial beam maximum, however, it becomes approxi- 
mately so for a narrow elevated beam. As an example, the theoretical distribution 
in elevation of echoes observed with a horizontal half-wave dipole situated A/4 


= 


7 ‘ 


, 


La , es 
—— Rate per unit elevation 





“90 


Fic. 2.—Theoretical distribution in elevation of sporadic meteors observed with a half- 
wave dipole aerial situated \/4 above ground. 





above a reflecting surface, for s=2, is given in Fig. 2. In this case ¢,,~40°, 
whereas the maximum of the aerial polar diagram is at 6=90°. The form of 
F(¢) will not, in general, depend strongly on s. 

The mean value of x (=) and the r.m.s. deviation from the mean can be 
evaluated in two steps. 
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vaM(xs) oc SP (D)rg%(x) 40, (x4) 


and let hg, %3, %3, Ps be the values of hy, %», xg, Pp, at F=0. From (5), (6), (7), 
a= ag5_ (8), pe=PgS.*%(8), xp=xg—§lnS,(9), x3=(h—hs)/H, 
hence, following Section 4.3 of Part I, 


: | S,2-0(9) In (8) a9 
&3—%,= 3 = me ’ 
J S.2¢-0(8) a6 





: | * §,4-2(8) In? S,(9) dd 
(5x3)? — (Sx_)? = — (Zs — Zp)? + 9 — . 





$2009) dd 


J—n 


Finally, we have . 
4a)ac{ F*%(4)09%( a) dp. 
Using (6), (7) and (11), we obtain 
Ps=PmG(¢), %3=x— §ln G(¢), 


G@)=(Se)(SeF)" sae): 


where 





sin dm 


hence, following Section 4.3 of Part I, 


fv" F(4)In G14) dp 


i—Z,= 





_ ’ 


), 


(5x)? — (8x5)? = —(%— Z5)? + 


/2 
F*-1($) dp 


,|, P@ Int) as 


9 rn/2 
9 [Fe -004) de 
Jo 

Equations (9), (16), (20), together with the computed value of 5x,, determine 
the r.m.s. height deviation, 5x (in units of the scale height), for a given aerial polar 
diagram S(¢,#). Similarly (8), (15) and (19), with the computed value of %,, 
enable the mean observed height to be related to the height h,,, at which the pressure 
is given, from (7), by 





(20) 


Pm = Ral 9%m COS" bn/4h,)"%, (21) 


om = K(Ah,,/sin b_)*”. (22) 

Since «,, appears in (21) as the cube root (and since f differs only slightly 
from h,,), it may be taken as the minimum detectable line density at range h/sin 4, 
in the direction (¢,,,0). The function F(¢) will in general be approximately 


symmetrical about ¢,,, which may therefore be taken as the mean observed 
elevation. 


where 
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3.1. Beamed aerial system.—As in Part I, we will assume various forms for the 
aerial polar diagram and evaluate ¢ and 5x as a function of s. We consider here 
the case of an elevated aerial beam with zero (or low) gain at ¢=0°, go°; thus 
¢,, will correspond approximately to the beam maximum and the controlling 
factor in the expressions for F(¢), G(¢) is S,(¢). We will anticipate that the 
values of & and 5x will not be strongly dependent on the exact form of these 
functions, and put F(¢)=G(¢) (as the beam width tends to zero, both F(¢) 
and G(¢) tend to S,(¢)). 

In order to obtain an estimate of the variation of the parameters of the height 
distribution with the form of the aerial polar diagram, we will consider several 
cases which between them should encompass the range of variation likely to be 
met with in practice (see Part I). They are: 

(a) F(¢)=(1 +66”), SA) =(1 +20)", 

(b) F(¢)=exp(—c.*¢"), S,(8) = exp (—c,"*), (23) 

(c) F(¢)=cos*(c,¢’), S,(9) =cos* (c,9), 
where ¢'=¢—¢,,, ¢, and c, are constants determining the beam width, and in 
(c) the functions are as specified for c,¢’ and c,# between — 7/2 and 7/2 and are 
zero otherwise. 

The integrals appearing in (8), (9), (15), (16), (19) and (20) for the above 
cases have been evaluated in connection with the shower meteor height 
distribution. The resulting values of the parameters % and 4x are given in 
Figs. 3(a) and 3(6) for a range of values of s. 

It is apparent from Fig. 3 that the lower the mean height in relation to 4,,, 
the greater is the r.m.s. deviation, the effect being most marked for the smaller 
values of s. The reason for this is that an aerial beam which falls slowly to zero 
at its extremities (e.g. Case (a) of (23)) tends to see relatively more large meteors 
which occur, in the main, below the maximum of the height distribution. This 
trend is also noticeable in the case of shower meteors discussed in Part I. 

3.2. Horizontal half-wave dipole, \/4 above a reflecting surface.—In this case 
we cannot express S(¢,#) exactly 4s a product S,(¢).S,(#); however, to a 
sufficiently good approximation we can take S,(¢)=S(¢,0), S,(#)=S(4,,, 9), 
where ¢,,~40°. Although when ¢=7/2, G(¢) now goes to infinity, the integrands 
in (19) and (20) all go to zero. 

The values of and 5x have been computed for s=2-0 and are given in 
Figs. 3(a) and 3(6). These points do not quite fit the general trend referred 
to in Section 3.1, the reason being that whereas ¢,, = 40°, the aerial beam maximum 
is at 6=90°. 

4. Height distribution for meteors observed within a selected elevation interval.— 
By selecting a narrow elevation interval about ¢,, (which in this case need not 
necessarily be the elevation at which F(¢) is a maximum) we can eliminate the 
final step in Section 3 which requires knowledge of S,(¢), F(¢) and G(¢). It is 
sufficient to restrict the range of ¢ to that for which G(¢) is substantially constant 
(say, to within 30 per cent), in which case the expressions on the right-hand sides 
of (19) and (20) are negligible and »* = v,*, = , 5x = 4x3, which can be evaluated 
from (15) and (16) with S,(#)=S(¢,,, 4). 

The parameters of the height distribution are given in Figs. 4(a) and 4(d) 
for the three forms of S,(#) defined in (23), as well as for a horizontal half-wave 
dipole, A/4 above ground (in this latter case for s=2-0 only). 
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Fics, 3 (a) and (b).—Parameters of the sporadic meteor height distribution. 


(a), (b), (c)}—Aerial polar diagram defined by (23) (a), (b), (c) respectively. 
O—Half-wave dipole aerial, 4/4 above ground. 
The broken curve in Fig. 3 (6) is the r.m.s. deviation, 8x,, for the elementary distribution v,*. 
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It will be noted that if the azimuths of the meteors were also known, the effect 
of the aerial polar diagram could be completely eliminated by selecting, in addition, 
a sufficiently narrow range of #, in which case = Z,, 5x=8x,. It would seem, 
however, that this would unnecessarily complicate the experimental apparatus. 


O5r 
(b)\ XO? 
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Fic. 4 (a) and (b).—Parameters of the height distribution of sporadic meteors observed 
in a selected elevation interval. 
(a), (b), (c)—Azimuthal aerial polar diagram defined as in (23) (a), (6), (c) respectively, 
©—Half-wave dipole aerial, 4/4 above ground (selected elevation ¢m= 40°). 
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5. Sporadic meteor mass distribution.—The total observed echo rate, N sec, 
can be related to the incident flux of meteors by integrating (3) with respect 
to v, h, ¢, #, over the complete range of variation of these quantities. 

If we have two equipments, A and B, with aerials of identical polar diagram, 
but with sensitivities defined by K, and K,, then using (3) and (5) with h,=h, 
the echo rates N, ,N, are related by 


N,/N, = (K,/Ka)*"(hip/hg)*2. 


The mean height will vary with equipment sensitivity to the same extent as the 
height h,,, hence from (21) and (22) 


h,~h,+4H1n(K,/K,), 
and since (h,—h,) will be small compared with h,~h,=h, 
N,/Np=(Ky/Ka)™ 11 + H(H/h) In (K,,/ Ky), 
where, from (7) of Part I, 
Kq/Ky=(€aP Ay?) /(€oPara®)- 


The second factor in the right-hand side of (24) will be close to unity, hence we 
may use an approximate value for H/h which we will take as 0-07 (h ~ 100km, 
H ~7km). 

McKinley (3) has measured the sporadic meteor rate as a function of trans- 
mitter power (in one case with a small change in wave-length). ‘The results, 
together with the values of s derived from them (using (24)), are given in Table I. 
The parameter, a, which defines the meteor magnitude distribution (see equation | 
(5), Part I) is also given. 

TaBLe | 
K,/Ka 
N, al No 
s 
a 


We will adopt the values 
$=2°00+0°02, a=2:51+0°04. (26) 


To a good approximation we could have ignored the second factor on the right- 
hand side of (24) and used N,/N,=(K,/K,)** leading to s= 1-99 + 0°02, which 
does not differ significantly from (26). 

McKinley has also measured the relative rates during the Quadrantid meteor 
shower, using the two equipments for which K,/K,=43. Owing to the large 
sporadic background on the high-power equipment he found that N,/N, varied 
during the period of observation. For the period og" 00™24>o00™E.S.T. 
(Canada) on 1951 January 3 he obtained N,/N,=25:2, while for the shorter 
period 12" 00™-195 00™E.S.T., when the radiant position was most favourable, 
he found N,,/N,=19°6. McKinley states that the hourly rate on the high-power 
equipment was twice the normal sporadic rate for the same hours during December 
and January. Depending upon the period of observation to which this applies, 
the true ratio for Quadrantid meteors becomes either 17:5 or 12°5, leading to 
$=1°76, a=2-02 or s=1-67, a=1°85 respectively. 

In the above we have assumed that all of the equipments used by McKinley 
were capable of detecting the short-duration decay-type echoes over most of the 
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solid angle covered by the aerial polar diagram. ‘That this must be so is indicated 
by the good agreement between the values of s deduced for the three different 
sensitivity ratios (Table I). Confirmation is provided by the values of «,, for the 
various equipments calculated from McKinley’s data (with ¢,, = 40°, G,, = 2-1 for 
the dipole aerial A/4 above ground, h~o5km, 7=1-0), namely 101, 2-5 x 10", 
40 x 10), 4-3 x 10% cm"!. 

6. Discussion.—As in the case of shower meteors, the width of the height 
distribution observed with a selected velocity group of sporadic meteors depends on 
the atmospheric scale height and the mass distribution of the incident meteors, 
but only slightly on the form of the aerial polar diagram. For the assumed mass 
distribution (inverse power law) it is independent of the equipment sensitivity. 
Since an accurate value for the exponent in the mass distribution has been obtained, 
the radio-echo height measurements provide a powerful technique for determina- 
tion of the atmospheric scale height in the meteor ionization region. The complete 
height distribution, v*, has not been evaluated explicitly ; however, it is clear that 
it will be of the same general form as for shower meteors (Part I), but slightly 
wider. 

The relation between the mean observed height, h, and the height, h,,, at 
which the pressure is p,, (equation (21)) has beenevaluated. Since k,ocv~* (Part I), 
Pm depends strongly on the selected velocity. It also depends on the probability 
of ionization, 8 (contained in k,) and upon«a,,. The form of the pressure-height 
variation may be obtained by integrating H-1 with respect to height, hence it is 
necessary to know f for one velocity only in order to obtain the pressure variation 
over a range of heights (limited by the range of meteor velocities, which is about 
20-60 kmsec~!). It follows that the probability of ionization as a function of 
velocity may be derived. 

The application of the radio-echo technique to the determination of the 
atmospheric pressure and scale height will be dealt with by S. Evans in a later 
paper.* 

From measurements by McKinley, using equipments with varying sensitivities 
corresponding to «,, between 10! and 4x10", the exponent in the sporadic 
meteor mass distribution is found to be s = 2-00 + 0°02, i.e. the absolute frequency 
increases by a=2-51+0-04 times per visual magnitude. If we assume that a 
meteor of zenithal magnitude +5 produces about 101” electrons per cm, then it is 
clear that the inverse-square mass law must hold accurately over the range of 
zenithal magnitudes +5 to +10. Theabove result is in good agreement with the 
most recent value given for telescopic meteors by Watson (4), namely a= 2:5, but is 
in marked disagreement with the value a= 2-08 derived by Williams (5). 

In contrast to the radio-echo measurements, the visual and telescopic observ- 
ations require large corrections, including subjective factors, in converting 
observed to absolute rates and values of a ranging from 2:0 to 4-0 have been quoted 
by various authors. 

The exponent in the mass law for the Quadrantid meteor shower is found to be 
1-67 or 1-78, depending upon the correction to the experimental results which 
should be made for the sporadic meteor background. It is clear that the exponent 
is considerably less than that for sporadic meteors, i.e. there is a relatively greater 
proportion of bright meteors in the shower, in agreement with visual and telescopic 
observations. 

* See p. 63. 
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An extension of the theory developed in these papers (I and IT), enabling the 
observed echo rate to be related to the absolute incident meteor flux, to the total 
rate of meteor ionization, etc., together with a more detailed investigation of the 
mass distribution of sporadic and shower meteors, will be presented elsewhere. 
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Summary 


A split beam radio-echo technique has been used to measure the heights 
of the reflecting points of meteor trails. The measured width of the 
distribution of heights has been used in conjunction with the theories 
previously developed by Kaiser to determine the atmospheric scale height. 
The observations of sporadic meteors have been divided into velocity groups 
and values of the scale height have been calculated for each. A value 
of the atmospheric pressure at the mean height of each group is calculated 
from the evaporation theory of meteors. For purposes of comparison these 
results are converted into atmospheric densities which are compared with 
those obtained from rocket flights and meteor photographs made in New 
Mexico and Massachusetts. The present results for scale heights and 
densities are in close agreement with those obtained in New Mexico. In 
view of this a latitude variation in the density profile, which has been suggested 
as an explanation of the difference between the New Mexico and Massachusetts 
results, now seems unlikely. 





1. Introduction.—This paper is concerned with the measurement of atmo- 
spheric pressure and scale height in the region 88 to 100km altitude, deter- 
mined from the observations of meteor heights with a radio-echo height-finding 
equipment. In two previous papers (1, 2) Kaiser has developed a theory of the 
distribution of heights of meteors observed on a radio-echo equipment. The first 
(referred to as Paper I) deals with shower meteors, and the second (Paper. II) 
with sporadic meteors. In these papers Kaiser, following Herlofson (3), shows 
that the pressure p, at the height of maximum ionization of a meteor of velocity v, 
producing a maximum of « electrons/cm path, is given by 

= 2gl cos** y (gHpua\"s 
AA ot Cz! () 
where 

H is the atmospheric scale height in centimetres, 

p. is the mean mass of a meteor atom in grams, 

B is the probability of ionization of an evaporated meteor atom, 

x is the zenith distance of the meteor trail, 

g is the gravitational acceleration in cm sec~?, 

lis the energy required in erg g~! to vaporize the meteoric material, 

A is the heat transfer coefficient (see Paper I), and 

A, the shape factor, is defined by A = %m-*?, where . is the cross-sectional 

area of the meteor in the line of flight, and m is the mass. 
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In principle this expression could be used to measure a series of values of p at 
different altitudes, but in the experimental measurement of meteor heights the 
reflecting point does not necessarily lie at the point of maximum ionization to 
which (1) applies. However, Kaiser’s theory shows that the mean height observed 
for a large number of meteors moving with the same velocity is close to the height 
at which the smallest detectable meteor produces its maximum ionization (referred 
to as the characteristic height, h,,). ‘The determination of h,, from the observed 
mean heights is discussed in Section 5. The smallest detectable value of the line 
density can be calculated from the parameters of the apparatus, and hence if ,, is 
determined for a series of values of v, the relationship between pressure and altitude 
can be obtained from (1), provided the other quantities are known. 

Kaiser also shows that the width of the height distribution is a function of the 
size distribution of the incident meteors and the atmospheric scale height, H. 
In Paper II Kaiser has derived a size distribution for sporadic meteors and this is 
used in Section 4 to determine the scale height from the measured height distri- 
bution. 

In the case of sporadic meteors, an integration over a range of zenith distances is 
taken in deriving the theoretical distribution of heights. The value of x to be used 
in calculating the pressure from (1) is shown to be the most probable elevation 
of the reflecting points. Since the equipment measures the individual elevations, 
it is sufficient to take the mean value of these. The remaining constants, /, A, A, 
» and B can be determined from data on photographic meteors given by Jacchia, 
and from correlations between visual and radio-echo observations. 

2. Apparatus.—The observations presented in this paper have been obtained 
with the radio-echo height-finding equipment described by Clegg and Davidson 
(4). The transmitter works on a frequency of 69 Mc/s and radiates pulses of 
6s duration and 50kW peak power from a single dipole which floodlights a 
major portion of the northern sky. A conventional receiver, switching rapidly 
between two aerials, is used. The ratio of the signals received in the two aerials 
enables the elevation and hence the height of the effective echoing point to be 
measured. The limitations and accuracy of the method have been discussed 
by Clegg and Davidson. Considerable improvement in accuracy has been attained 
in the present work by the use of a photographic recorder of the type described by 
Davies and Ellyett (5). The error in a single height measurement is + 1km. 
Some of the observations on major showers have been taken with a similar equip- 
ment using a slightly different aerial arrangement which can be set up at different 
azimuths, and working on a frequency of 36-6 Mc/s. The accuracy is of the same 
order. 

3. Observations.—Observations of the heights of sporadic meteors have been 
carried out using the 69 Mc/s equipment at Jodrell Bank during the period 
1952 January to July. Echoes for which it was possible to measure the velocity 
by the technique of Davies and Ellyett (5) comprised about 10 per cent of the total. 
Table I shows the number of meteors used in this analysis for each period of the 
day, foreach month. These were divided into the seven velocity groups shown in 
Table II. ‘The mean height, mean velocity and r.m.s. deviation of the heights are 
given here for each group. 

In Paper II it is shown that the scale height can be calculated from the width of a 
sporadic height distribution if the observations are selected to cover a range of 
elevations over which the gain of the aerials is roughly constant. Fig. 1 shows a 








No. I, 1954 atmosphere from radio-echo observations of meteors 65 


number distribution of the elevations of all the echoes given in Table II. The 
peak at 23° may well be due to the fact that at this elevation the ratio of the signals in 
the two aerials is unity and the chance of recognizing a velocity-type echo much 
greater. Accepting this explanation, echoes at elevations less than 20° or greater 
than 36° were rejected. These echoes are indicated by asterisks in Table II; the 
heights and r.m.s. deviations calculated after their exclusion are given in the lower 
section of the table. 


TABLE I 
Numbers of sporadic meteors used in this analysis 





Period of day oo” oom 06! oom— 12 oom 18" oom: 
062 oom 122 oom 182 oom 245 oom 





1952 
January 56 37 
February 28 21 
March 20 7 
April 32 17 
May 24 45 
June 13 dna 
July 10 7 





Totals 183 134 












































20 30 


Fic. 1.—Distribution of elevations. 


Ordinates: Numbers of meteors. 
Abscissae: Elevation in degrees. 


Table III gives similar data for four major showers. Further analysis and 
discussion of the astronomical significance of these results will be given in another 
paper. The data for the Perseid stream are used in Section 5 to estimate the 
behaviour of the aerial system. 
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TasLe II 
Heights of sporadic meteors 
Velocity group | 20-24} 25-29 30-34 35-39 49-44 45-49 50-59 
(km/sec) 
Mean velocity 22 27 32 37 42 47 53t 
(km/sec) 
Heights (km) 72 73 97 87 107 97 95 | IOI 103 97 100 | 108 103 
94 92 92| 84 107] 95 93 | 95 91] 99 102] 107 92 
93 65* 98 | 101 96 | 100 103 | 100 99] 100 go 99 IOI 
103 96 gI 88 100 83 105 96 105 97 99 99 107 
95 95 97| 95 890] 92 94] 100 103 | 93 101 | 100 94 
88 99 95 | 97 100] 98 97] 98 102] 89 99 | 103 113 
83 99 100*|} 87 100 87 94 89 99 | 102 106 99 97 
84 92 95 98 100 gI 100 92 105 | 100 99 | 102 I05 
go 88 84 94 89 99 9°] 100 102 | 106 g2 93 109 
78 96 95 | 102 103 94 105 97 103 | 101 89 98 100 
95 93 88| 82 91 | 94 94] 97 100] 96 97 | 101 108 
go 88 103 98 100 | 76* 92 | 101 96*| git 100] 118 92 
96 95 100] 110 92] 87 106] 98 113 89 76 | 103 92 
86 88 98 | 94 97| 91 77] 95 92] 100 93 | 130 108 
97* 82 87 87 86 79 94] 102 82 83 105 QI 103 
86 88 95 904 77 | 105 90] 100 99] 101 109 QI 100 
93 103 99 89 88 9° 97 | 100 95*| 100 103 | 109 86 
| 97 96 94] 102 95 | 85* 81 | 97 109 | 104 97 104 
| 89 95 88] 89 91] 96 99*| 100 107 100 104 
84 90 93| 95 89] 96 89] 106 92 81 94 
| go go 86} 89 81 | 90 94] 98 93 go 
83 92 87] 96 102] 116 88] ror 110 
88 g2 81] 101 85 95 92] 101 98 
88 9° 75| 90 91 | 96 97] 106 85 
| 97 84} 93 86| 92" 95 | 70 94 
| 9° 93 78 88} 113 90 Q2 103 
gt 84 | 103 79 | 116 99] 81 99 
89 82] 84 101 | 95 91] 94 93 
85 86] 83* 95 | 88 86] 96 84 
85 84] 79 97| 98 91 | 100 
100 IOI | 103 90 
97 93 | 91 100* 
103 95 | 100 94 
93 88 | 98 94 
1o1r 96 94 100 
86 99 94 
92 90 
96 93 
gI 
Number 24 60 71 77 59 35 41 
Mean height 89:2 90°7 93°2 94°4 97°6 97°3 100°7 
R.m.s. deviation 6-63 7°05 7°47 7°22 7°22 7°35 8°45 
Selected elevations 
Number 23 58 70 72 57 35 41 
Mean height 89°0 90's 93°3 94°7 97°6 97°3 100°7 
R.m.s. deviation | 6-58 6-12 7°44 7°13 7°36 7°35 8-45 


























+ There is some difficulty in recognizing high-velocity meteors on the 69 Mc/s height-finding 
equipment, causing the numbers to be asymmetrically distributed in this group. There are insufficient 


numbers in higher velocity groups as yet. 
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The velocities are given to the nearest km/sec from the results of radio-echo 
measurements made at Jodrell Bank (6, 7). 

4. Determination of scale height.—If the mass distribution of the meteors 
incident upon the Earth follows a law of the form dN=m~ dm, where dN is the 
number of meteors in the mass group m to m + dm, it has been established by Kaiser 
in Paper II that for sporadic meteors the value of s is 2-00+0-02. Referring to 
Fig. 4(5), Paper II, it can be seen that if 5/ is the r.m.s. deviation in a height 
distribution of sporadic meteors of uniform velocity and selected elevations, then 
for this value of s 

5h =(1-05 + 0°02) H, (2) 
where H is the atmospheric scale height. The probable error quoted is due 
partly to the uncertainty in the value of s and partly to the uncertainty in the aerial 
polar diagram. The latter error is small, but it will be shown later that a curve 
lying somewhere between Case (a) and Case (6) of Fig. 4(), Paper II, is most 


appropriate to the equipment used here. 


TaBLe III 





Shower 


Equipment 
used 


Number 
of echoes 


R.m.s. 
deviation 


(km) 





Daytime 
Arietids 
Perseids 


Geminids 


Geminids 


1952 
June 8-9 


1952 
Aug. 7-10 


1951 
Dec. 14 


1952 
Dec. 11-12 


69 Mc/s 


36 Mc/s 


69 Mc/s 


36 Mc/s 


$5 


Quadrantids 1953 
Jan. 2-3 36 Mc/s 


























The weighted mean value of the r.m.s. deviation of all the velocity groups in 
Table Il is 7-17km. Thus (2) gives, as a first approximation, a value of 6-80 km 
for the scale height, over the range of heights covered. 

In determining a more accurate value of the scale height, the effect of the finite 
range of velocities included in each group must be considered. Apart from v, 
the only quantity in equation (1) which may be velocity-dependent is f, the 
probability of ionization of a meteor atom. Adopting a law of the form 

B=av", where a is a constant, (3) 
then 
p=const./v?+"8, 
Hence 
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In estimating m, consider the two extreme velocity groups in Table II. 
Substituting the approximate value of H derived above we find (2 +/3)= 1-99. 
As a first approximation a value of n=o has been adopted, i.e. the probability of 
ionization is a constant and p varies inversely as v*. 

If 54 is the true r.m.s. deviation, let suffix m refer to measured values, 
suffix v to velocity variations, and suffix e to errors of measurement. Then 
dh? = d5h,,2—dh,2—5h,? since the distributions are approximately Gaussian. 
dh, is of the order of 1 km. 


TABLE IV 





Measured | Corrected 

M 
slits r.m.s. r.m.s. Number of 
(km/sec) deviation deviation observations 
(km) (km) 





6°58 6°43 23 
6-12 6-00 58 
7°44 7°32 : 7o 
7°13 7°02 72 
7°36 7°25 . 57 
7°35 7°24 ‘ 35 
8-45 8-30 , 41 


























The corrected r.m.s. deviations for each group are shownin TableIV. Values 
of the scale height have been calculated for each, using equation (2), and these are 
plotted against altitude in Fig. 2. The weighted mean value is 6-70 + 0-3 km at 
94°7 kmaltitude. The errors in the individual values of 54 used are of the order of 
(2N)-"? 8h, where N is the number of observations in a single distribution. This 
error is much larger than that due to the uncertainty in the value of s and the aerial 
polar diagram, but can clearly be reduced by taking more observations. The 
magnitude of the error is indicated by the horizontal lines in Fig. 2. 

The full line shows the linear relation fitted to these data by least squares and 
represents the values which will be adopted in this analysis. ‘The broken line 
shows the smoothed values adopted by the Rocket Panel (8). 

It is unlikely that the selection of a particular type of echo would have a serious 
effect on the height distribution. In order to check this, the r.m.s. deviation of 
92 sporadic velocity-type echoes of all velocities was found to be(9°5+1)km. The 
r.m.s. deviation of 106 sporadic echoes of all types obtained during the same period 
was found to be (g-8+1)km. If this small difference is found to be real as more 
data become available, it will be necessary to apply a small correction to the figures 
given in Table IV. In the meantime it is assumed that the velocity selection has 
no effect. 

5. Determination of characteristic heights.—As discussed in Section 1, the 
atmospheric pressure can be calculated at the characteristic height, h,,, which has 
to be determined from the measured mean heights. The theoretical deviation of 
the measured mean height from f,, is shown in Fig. 4(a), Paper II. 

For values of s of the order of 2 it can be seen that there is a difference of 0-12H in 
this deviation for the curves relating to the two extreme types of aerial polar 
diagrams considered. The polar diagram of the height-finding aerials cannot 
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easily be represented by a single curve, but the behaviour of the system can be 
derived in the following way. 

An appropriate value of s for the Perseid meteor stream is obtained first from 
the r.m.s. deviation of the height distribution given in Table III. This is 12 km 
or, using the value of scale height determined in Section 3, 1-65 H. Fig. 7 (6), 
Paper I, then shows the value of s to be 1°55 + 0°05. 

Unfortunately sufficient numbers of sporadic meteors of the same velocity 
as the Perseids are not available, but we can extrapolate from the highest velocity 
group available. It was shown in Section 4 that the pressure is inversely propor- 
tional to v*. 
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Fic. 2.—Atmospheric scale heights. 


The full line has been fitted by least squares to the present results. The broken line 
represents the smoothed values adopted by the Rocket Panel. 


Ordinates: Altitude in km. 
Abscissae: Scale height in km. 


I dinp 
We have ; Raha at 


Hence hy —h,= 2H 1n0,/v9. (5) 
Thus we should expect the mean height of sporadic meteors of velocity 60 km/sec 
to be 102-5 km. ‘The mean height of the Perseids is 100-3 km. The difference is 
2:2 km or 0-3H, and this we attribute to the difference in the value of s. A compari- 
son of Fig. 7 (a), Paper I with Fig. 4(a), Paper II, shows that the behaviour of the 
height-finding system is best represented by a line lying about half-way between 
curves (a) and (6). Thus for sporadic meteors of selected elevations the observed 
mean height is (0-15 + 0-05)H above the characteristic height. 
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The values of h,, for each velocity group are plotted in Fig. 3. 

The above argument presupposes that the behaviour of the 69 Mc/s equipment 
used for the sporadic observations and the 36-6 Mc/s equipment used for the Perseid 
observations is identical. ‘The polar diagram of the 36 Mc/s equipment in an echo 
plane of elevation 45° will be sufficiently near to the azimuthal polar diagram of the 
69 Mc/s equipment at constant elevation 27°-5. The low-frequency equipment is 
intrinsically of greater sensitivity but was operated at lower gain. In order to 
trigger the photographic recorder and provide a satisfactory height measurement 
the received echo must be of the order of five times noise. By measuring the peak 
power of the transmitter and the sensitivity of the receiver and using the formula 
of Lovell and Clegg (9) the minimum detectable line density at the mean elevation 
and mean height was found to be 2-3 x 101! electrons/cm for both equipments. 
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Fic. 3.—Characteristic heights of sporadic velocity groups. 


Ordinates: Altitude in km. 
Abscissae: Mean geocentric velocity in km/sec. 





6. Determination of constants.—Before equation (I) can be used to determine 
the absolute pressure/height relationship a knowledge of the constant 


Oo- <4(5) (6) 


is required. ‘This can be determined from photographic meteor data and radio- 
visual correlations. 
Jacchia (10) has given values of 


logio Ky, = — 11-75 + 0°02, 


A 
where K,,, = — and y is the drag coefficient, and log,, (yA)*® ty = — 9°43, where 75 is 


yl 
defined by 7 =7,v and7 isthe luminousefficiency. Thus(6) can be written 


a T 
logig O = logy 0°5 — logiy Ky, — $ logo (yA)>79 + $logio ef . (7) 
B 
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The visual magnitude of an object at 100 km emitting light of intensity J erg/sec 
uniformly in all directions is given by 


M = 24:6 — 2:5 logy I. (8) 


dmv? 
Fora meteor l=r =< (9) 


_ Bdm 
and a= ote ee (10) 


Tr 


Hence M = 25°35 —2'5 lo of) —10lo 19 U—2°5 logyy%. (11) 
810 B 810 


Greenhow and Hawkins (11) have shown that the zenithal magnitude of a 
meteor can be related directly to the electron line density and from the observations 
of Millman and McKinley on Perseid meteors a value of M=3 for «=7-3 x 10 
electrons/cm is quoted. An error of a factor 1-5 in « would produce an error of 
only 0-05 in logy). 

Hence from equation (11) logy, (*#) = — 31°05 for Perseid meteors. Sub- 


B 


stituting in equation (7) we find 


logio O = 4°24. 
This value will be adopted for all velocity groups in this analysis since the value of 8 
was shown to be independent of velocity in Section 4. 
7. Absolute pressure versus height.—Rewriting equation (1) we have 


cos? y 
v 


p=289 —Z~ (9Ha)". (12) 
Using an inverse-square law from the Earth’s centre, the value of g is 951 cm sec~? at 
95 km altitude. 
The value of « was shown to be 2-3 x 10" electrons/cm in Section 5. 
Referring to Fig. 1, the mean elevation is y=27°-5. Hence the pressure at the 
characteristic height is given by 


logy) P= 11°42 + $ log, H— 2 logy, v. 


The pressures caiculated from this equation are plotted in Fig. 4. The smooth 
curve was fitted to these points for the minimum algebraic sum of the residuals in 
log,)p. The form of the curve was obtained, by integration, from the scale height 
curve in Fig. 2. 

8. Comparison with other results.—In Fig. 2 the present measurements of scale 
heights are compared with the New Mexico rocket results. It can be seen that, ata 
given height, the difference in the scale heights is only of the order of the standard 
deviations. Thus there is no significant difference in the results of the two methods 
of measurement. 

In order to calculate the air temperature it is necessary to assume a value for the 
mean molecular weight of air molecules. At these altitudes partial dissociation 
of molecular oxygen may be expected, but to a first approximation we can use the 
molecular weight measured at ground level, which is 29 g/mole. Using H=6-70 
km at 95 km altitude the air temperature is 220deg. K. Because of the narrow 
range of heights available the temperature gradient cannot be estimated with any 
accuracy but it is between 2 deg. and 6 deg./km, increasing upwards. 
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Photographic meteor results fundamentally measure air densities, and thus for 
purposes of comparison the Jodrell Bank results have been converted into densities 


eit p=plg. (13) 
It will be seen that with a direct measure of H it is not necessary to use a value of the 
mean molecular weight in order to derive the gas density. 
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Fic. 4.—Atmospheric pressures. 
The circles represent absolute values of the pressure calculated from velocity groups. The form 


of the smooth curve which is adopted here was obtained by integration of the scale height curve 
in Fig. 2. 








Ordinates: Altitude in km. 
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Fic. 5.—Atmospheric densities. 
R, density profile from rocket and photographic meteor observations over New Mexico, latitude 32°N. 
M, profile from photographic meteor observations over Massachusetts, latitude 42° N. 
J, profile from radio-echo meteor observations at Jodrell Bank, latitude 53° N. 
Ordinates: Altitude in km. 
Abscissae: Log, density (g cnt), 
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Photographic meteor data obtained by Jacchia in New Mexico are now in 
agreement with the rocket data, but meteor data from Cambridge, Massachusetts, 
reduced in an identical manner, show a systematic difference. In Fig. 5, R is the 
most recent profile from rocket flights over New Mexico (12), M is the profile 
over Massachusetts from photographic meteor data, and J is the profile from the 
Jodrell Bank observations. 

It was pointed out by the Rocket Panel that the difference between the New 
Mexico and Massachusetts results may be due to a latitude effect or to the differing 
location with respect to upper atmospheric circulation patterns. The rocket site 
in New Mexico has a latitude of 32°, Cambridge, Massachusetts, has a latitude of 
42°, and Jodrell Bank, 53°. Thus, in view of the present results, a monotonic 
variation with latitude is unlikely. Results from other sites will clearly be awaited 
with great interest. Further observations to establish the seasonal effect at Jodrell 
Bank are now being carried out. 
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AN ANALYSIS OF BURSTS OF SOLAR RADIO EMISSION 
AND THEIR ASSOCIATION WITH SOLAR AND 
TERRESTRIAL PHENOMENA 


R. D. Davies 
(Received 1953 December 24)* 


Summary 


Solar radio emission records on seven frequencies in the range 60 to 
10 000 Mc/s were studied over a period of 18 months. The analysis, which 
includes a number of histograms, shows that many of the properties of bursts 
change with frequency. A study of the time delay between bursts on different 
frequencies revealed that 3000 Mc/s bursts often occur first; a simultaneous 
up and down movement of ionized material from the level of zero refractive 
index at that frequency is postulated. It is found that bursts were more 
frequent than flares, fadeouts and crochets and that they almost always 
accompanied these effects. The commencement of bursts appears to be 
simultaneous with that of flares and crochets but precedes that of fadeouts 
by about two minutes. In addition, there is a rough correlation between the 
intensity of bursts, flares, fadeouts and crochets. 





1. Introduction.—In the course of routine observations of radio emission from 
the Sun taken at the Radiophysics Laboratory, Sydney, a considerable accumula- 
tion of records has become available covering a large frequency range. Although 
a day-by-day summary of these records is published in the Quarterly Bulletin 
of Solar Activity, there is no opportunity to present in that publication a detailed 
analysis of outstanding occurrences, in particular the bursts of radio emission 
characteristic of the records. Such an analysis seems worth while in order to find 
what generalizations emerge from the data and how the recorded bursts are 
related to other solar and terrestrial phenomena. 

Previous work on this subject includes detailed studies of (a) many bursts at 
one frequency (1, 2, 3), (4) occasional observations of simultaneous bursts on several 
frequencies (4, 5), and (c) spectrum analyser observations at the lower frequencies 
(6, 7, 8). However, no systematic investigation of bursts recorded over a wide 
frequency range has been described. 

In the present paper a detailed analysis is made of bursts recorded in the period 
1950 January to 1951 June. The records were taken at seven frequencies in the 
range 60 to 10000 Mc/s. The analysis includes a detailed study of the correlation 
of bursts with flares, sunspots, ionospheric fadeouts and magnetic crochets. In 
most cases no attempt has been made to draw physical conclusions from the 
analysis. However, certain results which appear to be of immediate interest are 
discussed briefly in the final section. 

The main analysis is given in Section 4. In treating a long series of largely 
disconnected properties, no attempt has been made to maintain continuity or to 
discuss each item fully. The essential results are presented in the accompanying 
diagrams and tables and only the main conclusions are stated in the text. 


* Received in original form 1953 August 24. 
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2. The data.—The radio data, most of which were tabulated by Davies and 
Hindman (9g), were obtained mainly from records at 200*, 600, I 200, 3 000 and 
9 400 Mc/s; these were extended to 62 and 98 Mc/s in a later part of the analysis. 
The following burst properties were measured as a preliminary to the analysis. 

(i) Time of start—the time when the radio noise rose to 5 per cent of the 
maximum burst intensity. 

(ii) Time of upswing—the time at which the major part of the burst (lasting 
more than 10 per cent of the total life-time) rose to 20 per cent of the peak intensity. 

(iii) Maximum (peak) burst intensity. 

Flare and sunspot information was obtained from the Quarterly Bulletin of the 
I.A.U. and the Monthly Bulletin of the Tokyo Observatory. 

Fadeouts, observed as a fall in intensity of the received galactic radiation, were 
derived from 18 Mc/s galactic noise records obtained at this laboratory. Observed 
fadeout intensities measured in decibels were corrected for obliquity of ray path 
through the ionosphere by multiplying by the secant of the Sun’s zenith distance. 

Magnetic solar flare effect (crochet) information was obtained from the magneto- 
grams of Toolangi Magnetic Observatory. The quarterly reports in the Journal 
of Geophysical Research and the 1950 Bulletin of 1.A.'T.M.E. were also used. 

Ina few cases complete burst, flare, fadeout and crochet information is available. 
Two typical examples are given in Fig. 1. In other cases the data are incomplete 
but are used in the statistics. 

3. Qualitative properties of bursts. (a) Definitions.—The occurrence called a 
burst in this paper is any clear-cut solar radio emission rising above the daily 
level. In view of the close relationship of bursts with flares (discussed in 
Section 4, Part B), many may be identified with the “ outbursts”’ studied by other 
workers (6, 10, 11). 

(b) Burst outline.—A characteristic of some bursts is their “‘jaggedness ”— 
short-duration rises and falls in intensity. The degree of jaggedness decreases 
with increasing frequency, as can be seen in Fig. 1. In an average burst at 200, 
600 or I 200 Mc/s about two peaks were observed each minute. However, bursts 
at 3000 and 9 400 Mc/s consisted of a smooth rise and fall. This decrease of the 
amplitude of fluctuations with increasing frequency was also noted by Covington 
and Medd (12). 

Burst traces at the lower frequencies may be considered to consist of fluctua- 
tions superposed on a smooth minimum level. ‘The smoothed minima are similar 
in shape to the burst outlines at the higher frequencies. 

Burst activity is more profuse at 200 Mc/s than at the higher frequencies. 
Over some periods as many as 25 separate bursts, most of which were not 
coincident with activity on the higher frequencies, were recorded each hour. 
Therefore the statistics applied to 200 Mc/s in Part A of the analysis (Section 4) 
relate only to bursts which coincided with those on the higher frequencies or with 
flares. In sections where all bursts (irrespective of their coincidence and flare 
associations) are analysed, the 200 Mc/s data have been excluded altogether. 
This applies to Sections 4 A (a) and 4 B(a). 

A number of bursts of long duration were examined to determine the extent to 
which the shape was preserved on different frequencies. ‘The smoothed outlines 
were similar on all frequencies and also the similarity was more pronounced the 
closer the frequencies. The fine detail, as seen on high-speed records, does not 


* The 200 Mc/s records came from the Commonwealth Observatory, Mt. Stromlo, 
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Fic. 1 (a).—Burst, fadeout, crochet and flare observations on 1950 July 21. 
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Fic. 1 (6).—Burst, fadeout, crochet and flare observations on 1951 June 22. 
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correspond on different frequencies. These effects are illustrated in Fig. 2, 
which shows a burst occurring on 600 and 1 200 Mc/s. 
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Fic. 2.—(a) Outline of burst on 1950 February 17, at 600 and 1 200 Mc/s. 
(6) A comparison of a section of the above burst at 600 and 1 200 Me|s. 

4. The analysis.—The analysis consists of a systematic study of the properties 
of bursts and their correlation with flares, sunspots, fadeouts and crochets. The 
analysis has two parts: 

(A) Relations between bursts on different frequencies. 

(B) Relations between bursts and sunspots, flares, fadeouts and crochets. 

Part A: Relations between bursts on different frequencies. (a) Rate of occurrence.— 
The number of bursts occurring per hour of observation (and also the average 
interval between bursts) on the frequencies 600, I 200, 3000 and 9400 Mc/s is given 
in Table I. It is evident that bursts become progressively less frequent with 
increasing frequency. The value of one burst at 3000 Mc/s for each 17 hours of 
observation is in fair agreement with Covington’s value of 23 hours estimated at 
2800 Mc/s for the period 1947 to 1950 (1). 


TABLE I 


The average number of bursts per hour observed at 600, 1200, 3000 and 9400 Mc/s in the 
period 1950 January to 1951 June 


Average interval 
Frequency Average number of b b 
(Mc/s) bursts per hour oe ae 
(hours) 
600 o'10 10 
1200 0:08 12 
3000 0:06 17 
9400 0°04 26 


(b) Lifetime of bursts.—Fig. 3 shows the distribution of the lifetimes (defined as 
the time between 4, peak intensity points) of bursts recorded on the five fre- 
quencies 200, 600, 1200, 3000, 9400 Mc/s. It is seen that the proportion of bursts 
of short duration (say less than four minutes) is greatest at the lower frequencies. 
On the other hand, an examination of individual bursts showed no clear relation- 
ship between lifetime and frequency. 
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(c) Time of decay of bursts.—The distribution of averaged “‘decay times” 
(determined by fitting sections of the tail of each burst to an exponential form 
e~") for 40 bursts is plotted in Fig. 4. The mean decay timeat the low frequencies 
is shown to be less than at the high frequencies. However, it was observed that 
the smoothed minima of the low frequencies had decay times of the same order ‘as 
the high frequencies. At no frequency did the majority of bursts show an 
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Fic. 3.—The distribution of lifetime of bursts at 200, 600, 1 200, 3 000 and 9 400 Mc/s. 
exponential decay over the whole length of the “tail”. The manner in which the 
decay time changed over the tail of a burst did not vary significantly with frequency. 

(d) Intensity of bursts—Histograms showing the distribution of peak burst 
intensities at various frequencies are shown in Fig. 5 (a). ‘The modal values are 
plotted as a function of frequency in Fig. 5(5). The mean spectrum obtained by 
joining these points (full line) is see 1 to show a minimum at 1200 Mc/s in contrast 
to the smooth downward trend of the intensity of the quiet Sun (dotted line). 
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Fic. 4.—The distribution of decay times of bursts at 600, 1 200, 3 000 and 9 400 Mc/s. 
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However, the intensity of individual bursts showed no unique trend with 
frequency. The spectra of 44 bursts were distributed as follows: 
8 showed a downward trend with increasing frequency, 
7 an upward trend, 
7 a maximum at an intermediate frequency, 
14 a minimum at an intermediate frequency, and 
8 fell to a minimum at the lower frequencies and showed a maximum 
at the high frequencies. 


(e) Time delays between bursts on different frequencies. (i) Time of start.—In order 
to discover the time relationship between bursts, delay times were calculated 
relative to 1200 Mc/s, the middle frequency. The distribution of the delay of 
times of start is shown in the dotted histograms in Fig. 6. They are symmetrical 
about coincidence on all frequencies except 9400 Mc/s, which shows a lag. This 
unexpected 9400 Mc/s lag was further tested by comparison with 3000 Mc/s. 
The times of start and maximum intensity of 9400 Mc/s were found to be delayed 
4 minute behind 3000 Mc/s, 
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(ii) Time of upswing.—Most bursts at 200 and 600 Mc/s which preceded those 
at 1200 Mc/s in the dotted histograms of Fig. 6 showed preliminary activity of low 
intensity before the sudden upswing characteristic of the lower frequencies. The 
results of a comparison of times of upswing (as defined in Section 2) with 1 200 Mc/s 
are given in the full-line histograms in Fig.6. Although a large number of bursts 
are coincident within the error of measurement there is a tendency for 200 and 
600 Mc/s to lag behind 1 200 Mc/s. 
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Fic. 5.—(a) The distribution of the intensities of 44 bursts on 200, 600, 1200, 3000 and 9400 Mc/s. 
(b) The full curve gives the modal intensities derived from Fig. 5 (a) plotted against 
frequency. The dotted line gives the intensity spectrum of the quiet Sun. 

The sequence of times of upswing on different frequencies was determined for 
35 bursts. They were distributed as follows: 

(a) 7 coincident on all frequencies. 

(5) 4 lowest frequencies occurred first. 

(c) 15 highest frequencies occurred first. 

(d) 9 anintermediate* frequency occurred first. 

(e) ©Oan intermediate* frequency occurred last. 
1200 and 3000 Mc/s bursts occurred first most often. Double bursts with two 
clear-cut upswings separated by several minutes sometimes occurred at the lower 
frequencies. In some cases, all of group (4), the first burst of the pair appeared 
to fit the sequence best; in other cases, all of group (c), the second. The latter 
effect was observed by Reber (5) in one occasion. 

Part B: Relations between bursts and sunspots, flares, fadeouts and crochets. 
(a) Sunspot statistics and the rate of occurrence of bursts—To make possible a 
comparison of the radio data with sunspot data, an estimate was made of the 
average number of bursts per hour foreach month. In Fig. 7 are displayed three- 
month running means of 600, 1200, 3000 and 9400 Mc/s burst indices along with 
monthly average sunspot areas and numbers.} It is concluded from the figure that 


* A frequency other than highest or lowest recorded. 
+ Obtained from Washington Naval Observatory circulars. 
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(i) time plots of the burst index are similar at different frequencies; and 

(ii) the burst index correlates well with sunspot area and number ; the number 
of bursts per hour at any frequency may be used as an index of solar activity. 

(b) Coincidence of bursts, flares, fadeouts and crochets. (i) Bursts and flares.— 
The distribution of 261 coincidences between bursts and flares is given in Table II. 
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Fic. 6.—The dotted histograms give the distributions of times of start on 200, 600, 3000 and 
9400 Me|s relative to 1200 Mc/s. The full-line histograms are for the times of upswing. 
Three flares only were unaccompanied by excess radio noise at one or more 
frequencies in the range 60 to 10000 Mc/s and in two of these cases the radio 
records were incomplete. Thus at least 99 per cent of all flares are accompanied 
by excess radio emission in some part of the frequency range. However, the 
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possibility that this high correlation could be due in part to chance coincidences 
requires special examination, especially when activity was recorded only on one of 
the low frequencies where activity may persist over long periods. As shown in 
Table II, 31 out of the 222 flares were accompanied by extended activity on 
200 Mc/s. In these cases it was difficult to recognize a discrete burst and it is 
possible that in some of them a significant coincidence did not exist. 

These results extend those of Dodson, Hedeman and Owren (13), who observed 
that 78 per cent of all flares were associated with energy excess at 200 Mc/s. 
Flares are more often coincident with low- than with high-frequency bursts. 


9400 Mc/s 





Mc/s 




















Ton Feb Max Ape Moy he yy Sept Nor Dec Feb Mor Ape. May June 


1951 


Fic. 7.—A comparison of the three-month running means (number of bursts per hour) with 
the monthly average sunspot areas and numbers for the period January 1950 to Fune 1951. 


TABLE II 
The distribution amongst the classes (a), (b) and (c) of flares and bursts on each of the frequencies 
200, 600, 1200, 3000 and 9400 Mc/s 
Class (a) Class (5) Class (c) 
Flares with Flares without Burstswithout Totals 
bursts bursts flares 


200 Mc/s number 180 36 29 245 
percentage 74 14 12 100 
600 Mc/s number 109 61 34 204 
percentage 53 30 17 100 
1200 Mc/s number 94 79 28 201 
percentage 47 39 14 100 
3000 Mc/s number 21 50 4 75 
percentage 28 67 5 100 
9400 Mc/s number 40 44 3 87 
percentage 47 50 3 100 
All frequencies number 211* 11t 39 261 
: percentage 81 4 15 100 

* These included 31 flares accompanied by extended activity on 200 Mc/s. 

+ 8 of these 11 have associated burst activity on 98 and 62 Mc/s. 


However, a burst on a high frequency is more likely to be accompanied by a flare 

than one ona low frequency. Thirty-nine bursts with a wide range of intensities 

were unaccompanied by flares of importance 1~ (the minimum value reported) or 

greater. In order to test whether these bursts were accompanied by less intense 
6* 
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flare-like activity, 200 Mc/s records were searched at the times of flares of 
importance 0.* Of 20 flares 0, 7 were associated with distinct bursts and 7 
with prolonged intense activity. Thus, less pronounced forms of flare-like 
activity can produce excess radio noise. 

(ii) Fadeouts, bursts and flares.—Sixty-eight galactic fadeouts (sudden absorp- 
tion by the ionosphere of galactic radiation) were reported in the period 1950 July 
to 1951 June. ‘They are compared with Tokyo flare records and burst data in 
Table III. Of 68 fadeouts, only one was not coincident with a burst, and then the 
fadeout intensity was small. Fadeouts are more often associated with bursts than 
flares; this may be due in part to the fact that flares are more easily missed, and in 
part to effects intrinsic to these phenomena. 





TasLe III 
The distribution of 68 fadeouts among the 9 classes of flare and burst conditions 
Bursts 
Occurrence Nor-occurrence No record 
Occurrence 25 
Flares < Non-occurrence II 
No record 31 I 


In order to test the association of a fadeout with each burst, the galactic radio 
records were searched at the time of 20 intense bursts. In 15 cases fadeouts 
(larger than average) were detected and in 5 no fadeouts were visible on the records. 
It is concluded that there are more bursts than fadeouts and that all fadeouts are 
probably accompanied by a burst on some frequency. 

(iii) Crochets, flares and bursts—The coincidence between flares, bursts and 
the crochets visible on the magnetograms of the Toolangi Magnetic Observatory, 
near Melbourne, are given in Table IV. The feature of this table is that more 
crochets are associated with bursts (67) than with flares (35). This is due to the 
smaller optical observing time and to the fact that bursts occur without an observed 
flare. 

TaBLe IV 
The distribution of 92 crochets amongst various classes of flare and burst conditions 
Bursts 


A. 





~ 


Occurrence Non-occurrence No record ‘Total 


{ Occurrence 29 ee 6 35 
Fasee 4 Non-occurrence II I 2 14 
| No record 27 I 15 43 
| Total 67 2 23 92 


(c) Time delays between bursts, flares, fadeouts and crochets. (i) Bursts 
and flares.—The differences between the times of start of flares (I.A.U. data) 
and the times of upswing of the corresponding bursts on 200, 600, 1200, 3000 
and 9400 Mc/s were determined and the distributions plotted in Fig. 8. The 
histograms show no marked asymmetry at any frequency. ‘This agrees with 
the results of Dodson, Hedeman and Owren (13) who found that the starting 
times of most types of bursts seem to cluster about the rise in H« brightness. 


* Phenomena less intense than flares of importance 1~, recorded in the Mount Stromlo 
spectrohelioscope log. 
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(ii) Bursts and fadeouts.—To determine the time relationship between fadeouts 
and bursts the times of starts of the galactic fadeouts were compared with the 
burst upswing on the frequency which occurred first. The distribution of time 
differences is shown in Fig. 9 along with the fadeout-flare comparison. The burst- 
fadeout histogram (Fig. 9 (a)) shows a smaller spread than the burst-flare histogram 
(Fig. 8) or the flare-fadeout histogram (Fig. 9(6)). The median delay of fadeouts 
after bursts is I-5 minutes. 
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Fic. 8.—The distribution of the times of start of flares and the upswing of the corresponding bursts 
on 200, 600, 1200, 3000 and 9400 Mc/s. 


Fig. 10 shows an example of a burst preceding a fadeout. Hey, Parsons and 
Phillips (3), however, found that on the average the commencement of radio- 
communication fadeouts, as distinct from galactic fadeouts, preceded by 13 
minutes bursts at 4-1 metres in the period 1946-47. 
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(iii) Bursts and crochets.—The distribution of the differences between the time 
of upswing of the burst (on the frequency which occurred first) and the time of 


start of the corresponding crochet is given in Fig. 11. The median delay of the 
crochets after the bursts is about } minute. 
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Fic. 9.—(a) The distribution of time differences between the start of a fadeout and the upswing of 
the burst which occurred first. 
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(6) The distribution of time differences between the commencement of flares and fadeouts. 
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Fic, 10.—An 18 Mc/s galactic record showing a fadeout preceded by a burst. The dotted 
curve is a record from an aerial directed away from the Sun (the burst is consequently smaller and 
the beginning of the fadeout is visible). 


(d) Intensity relations between bursts, flares, fadeouts and crochets. (i) Burst 
intensity and flare importance.—The distribution of burst intensities among the 
4 flare classes (no flare, and flares of importance 1, 2 and 3) is shown in Fig. 12. 
Flares of the various classes are seen not to be restricted to any burst intensity 
range; nevertheless, small flares tend, on the average, to give smaller bursts. 
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Fic. 11.—The distribution of time differences between the start of the crochet and the upswing on 








the frequency which occurred first in the related burst. 
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Fic. 12.—The distribution of flares amongst the various ranges of burst intensity at 200, 600, 
1200, 3000 and 9400 Mc/s. 
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(ii) Fadeout and burst intensities.—The corrected* value of fadeout intensity is 
compared with the peak intensity of the largest simultaneous burst observed in the 
range 600 to 9400 Mc/s in Fig. 13. Although a close correlation does not exist, the 
intense fadeouts tend to be associated with large bursts. 
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Fic. 13.—A plot of corrected galactic fadeout intensity against the intensity of the largest burst 
which occurred in the range 600 to 9400 Mc/s. 
, 


Correlation coefficients were calculated between corrected fadeout intensity 
and 

(a) intensity of largest burst in the range 600 to 9400 Mc/s; 

(5) 200 Mc/s intensity ; 

(c) flare importance (evaluated as follows, 1- = 0-67, I =1, I” = 1°33, etc.); 

(d) Ha linewidth of flare. 
It can be seen from the results in Table V that no index of solar activity is more 
closely related to fadeout intensity than is flare importance. 


TABLE V 
Correlation coefficients between galactic fadeout intensity and burst and flare indices 


: P ' ’ i : Standard 
Quantity correlated with fadeout intensity Correlation coefficient dusletiem 
600 to 9 400 Mc/s burst intensity o'51 +o°12 
200 Mc/s burst intensity 0°40 +0'16 
Flare importance 0°56 +0°14 
Ha linewidth of flare O14 +0°25 


(iii) Crochet and burst intensities.—No clear correlation was observed between 
the crochet H-element intensity and either (a) maximum burst intensity in the 
range 600 to 9400 Mc/s, or (b) intensity of 200 Mc/s bursts (shown in Fig. 14 (a)). 
McNish (14) has postulated that the size of a crochet is roughly proportional to the 
simultaneous daily magnetic variation. We should then expect that the quotient, 


* See Section 2. 
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H-element intensity divided by the daily variation, might show a correlation with 
burst intensity. Fig. 14(b) shows that high quotient values tend to occur only 
with large bursts, though a detailed correlation is not present. 

A note on the correlations between bursts and the geomagnetic phenomena, 
magnetic storms and world-wide crochets, is given in the Appendix. 
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Fic. 14.—(a) H-element intensity of crochet plotted against 200 Mc/s burst intensity. 
(6) Quotient (H-element intensity/daily variation) plotted against 200 Mc/s burst 
intensity. 
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5. Discussion. (a) Change of characteristics with frequency.—The outstanding 
feature is the occurrence of the jaggedness at only the low frequencies and the 
consequent smaller decay times observed at those frequencies. The jaggedness 
decreases in intensity with increasing frequency while the smoothed (minimum) 
burst increases in intensity. At the extremes of frequency one type of burst 
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activity can occur without the other. ‘Thus it is suggested that there may be two 
separate components of bursts, one which shows rapid fluctuations and prepon- 
derates at the low frequencies, and one which is smooth and is characteristic of the 
high frequencies (although it may occur at low frequencies also). The latter has a 
spectrum similar to that of the quiet Sun and may possibly be thermal in origin. 
It is possible that plasma oscillations give rise to the jagged burst component.  __ 

(b) Delay times.—Time delays of the order of minutes between bursts on 
different frequencies have been observed by a number of workers (15, 12, 6, 4, 5» 
16). These time delays are thought to be due to the outward movement of ionized 
material through successive ‘‘ plasma” frequency levels. Such a movement, 
followed by an inward movement, has been detected by radio interferometric 
methods (11) and by visual methods (17). A number of bursts studied in the 
present analysis showed a frequency sequence which corresponded to a movement 
outwards and subsequently inwards similar to that detected by Payne-Scott and 
Little (11). Another group of bursts suggested a movement of exciting material 
towards the Sun. This is consistent with Ellison’s observations (17) of surge 
prominence material. 

The statistical analysis showed a tendency for all frequencies to lag behind 
3000 Mc/s. Many individual bursts showed a similar effect. This time sequence 
would suggest a simultaneous movement up and down from the level of zero 
refractive index at the frequency which occurred first. At 600, 1200, 3000 and 
9400 Mc/s these levels were calculated by Smerd (18) to be 10 000, 5.000, 3000 and 
500 km for the quiet Sun. These values would need to be increased by a factor of 
approximately five to allow for the effect of coronal condensations, supposed by 
Waldmeier and Miller (19) to occur above sunspots (from the vicinity of which the 
ionized material is ejected). The revised values approximate to the heights of the 
top of flares, at which level surge prominences are known to originate (Ellison (17)). 
Dodson and McMath (20) give supporting evidence for the suggestion that the 
radiating material moves up and down simultaneously from some levels in the 
solar atmosphere. On 1951 May 8 they observed visually a surge prominence 
associated with a limb flare, the top of which was at 10000km; ionized material 
rose to 50 000 km at 700 km/sec and then to 200 000 km at 300 km/sec; in addition 
there was a simultaneous movement of material downwards from the 5000km 
level. 

(c) Association between bursts, flares and terrestrial flare accompaniments.— 
During the period of radio-frequency observations between 1950 January and 
1951 June there were 400 bursts, 300 flares, 92 crochets and 68 galactic fadeouts 
(the latter cover one year only). The difference in the frequency of occurrence of 
these phenomena is due to the following effects : 

(a) Method of observation—e.g. overcast conditions limit flare observations 
and magnetic storminess prevents the detection of crochets. 

(5) Intrinsic differences in the frequency of occurrence—even under favourable 
observing conditions more bursts can be detected than flares, fadeouts or crochets. 
Nearly all flares and their terrestrial accompaniments are coincident with bursts. 
(The proportion of flares, fadeouts and crochets unaccompanied by bursts is 
I per cent, I per cent and 2 per cent respectively). Hence bursts are a sensitive 
indicator of solar activity and may be used to confirm doubtful fadeouts and 
crochets. Eighty-one per cent of bursts are simultaneous with flares of importance 
I~ or greater. Some of the remaining bursts are accompanied by less intense 
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forms of flare activity while others accompany phenomena (fadeouts and crochets) 
normally associated with flares. It is thus concluded that bursts are caused not 
only by flares of importance 1~ or greater but also by less intense flare-like activity, 
perhaps in the bright floccular regions near sunspots. 
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APPENDIX 


Magnetic storms and world-wide crochets.—Normally it is difficult to establish 
the association between a flare or burst and the sudden commencement of a 
magnetic storm a day later. However, if these phenomena occur in a quiet 
solar period, the relationship may be detected. On four occasions tabulated in 
Table VI, a sudden commencement was preceded by a burst although no flare was 
observed. This substantiates the conclusion above that weak flare-like activity 
(which can produce the sudden commencements) as well as larger flares may 
cause bursts. 


TaBLe VI 
Sudden commencement of magnetic storms accompanied by bursts but not by flares 


x Time delay 
Sudden commencement of a Burst (U.T.) Ste Tee Flare records 
magnetic storm (U.T.) and S.C ILA.U. Tokyo 


1950 
Sept. 154 10 ;gm Sept. 147 05% 05m 298 13 N.F.* N.F.* 


5% ogh 22m Feb. 49 072 24m 204 53m N.Rt N.R-+ 
164 062 52m April 174052 21m 25h 3;m N.F.* N.R.F 
2549 ogh 27m June 249072 15m 22h ;2m N.F.* N.R.t 


*N.F. No flare observed although the Sun was watched, 
+ N.R. No recording with spectrohelioscope. 


It is known to geophysical workers that a crochet observed in the sunlit 
hemisphere of the Earth may sometimes be accompanied by a simultaneous effect 
in the dark hemisphere. Such world-wide crochets are rare and their association 
with flares is not well established. Sixteen world-wide crochets observed at the 
Toolangi Magnetic Observatory were associated with observations at daylight 
stations; a flare was simultaneous with only one of these. Eleven world-wide 
crochets not accompanied by a reported flare were associated with bursts and are 
listed in Table VII. 
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TaBLe VII 
A list of world-wide crochets accompanied by bursts but not by flares 

Time (U.T.) 1950 Time (U.T.) 1950 

Feb. = 174. 23" o8™ July = 214 oo 46m 
March 254 04" 53™ Aug. 24 22h o4m 
May 52 052 som Aug. 39 orb 20m 
May 204 20h ;5m Aug. 94 oo! 04m 
May 204 20! 36m Nov. 15% 06% 45™ 


July —_ 197 oo 33™ 
These phenomena thus appear to be genuine and may be supposed to have a 
flare-like origin judging from their association with bursts. 


Division of Radiophysics, Jodrell Bank Experimental Station, 
Commonwealth Scientific and Industrial Lower Withington, Macclesfield, 
Research Organization, Cheshire. 


Sydney, Australia: 
1953 December 21. 
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THE HELIUM STAR HD 168476 
A. D. Thackeray 


(Received 1953 October 27) 


Summary 


Wave-lengths of 206 absorption lines have been measured in HD 168476 
between 3700 and 6700 A. No lines of hydrogen or oxygen are found, while 
lines of He1, Cir and Nel are prominent. A list is given of the various 
elements represented. A large range of excitation and ionization potentials is 
exhibited. Titt, with the same I.P. as H and O1, is definitely present. 
Except for H and O, a fair match with the spectrum of the B5 supergiant 
67 Oph can be made, but equivalent widths in HD 168476 are about 
twice those in the supergiant; this corresponds roughly to increasing 
the number of atoms above the photosphere by a factor 10. The observations 
are consistent with a great deficiency of hydrogen causing increased trans- 
parency and the spectral characteristics of a supergiant. 

The velocity —165 km/s appears to be constant. In view of the large 
positive velocity of HD 160641, a similar H-weak star, these stars must be 
regarded as high-velocity objects rather than as exhibiting peculiar atmo- 
spheric motions. Proper motion data suggest an absolute magnitude of —2 
or brighter; this is supported by evidence from a line in HD 168476 which is 
probably due to interstellar calcium. 





1. Introduction.—There are now three high-temperature stars known in whose 
spectra no lines of hydrogen, in emission or absorption, can be found. 
These stars, in order of discovery, are: 


HD m, Sp. Reference 
124448 101 B2_ OD. Popper, P.A.S.P., 54, 160, 1942; 59, 320, 1947. 
160641 98 O - W. Bidelman, Ap. 7., 116, 227, 1952. 
168476 8 9 B_ A.D. Thackeray and A. J. Wesselink, Obs., 72, 248, 1952. 


The present paper is concerned with the spectrum of the brightest of the three, 
HD 168476, as recorded with the Cassegrain 2-prism spectrograph of the Radcliffe 
74-inch reflector. Seventeen plates with the //3-7 camera (49 A/mm at Hy) and 
five plates with the //6 camera (29 A/mm at Hy) are available. The latter plates 
have been measured thoroughly for identification of absorption lines, while the 
lower dispersion plates have been used primarily for measures of velocity. One 
two-hour exposure with the //6 camera, taken in excellent conditions, is heavily 
overexposed in the ordinary photographic region, but permits a study of the 
interesting region 3700-3900 A at the relatively high dispersion of about 16 A/mm. 
Another plate on 103 aF emulsion permits a study of the region 5000-6700 A. 

2. Identification and representation of elements.—A total of 206 absorption lines 
have been measured between 3700 and 6700 A. Of these, 39 lines, in addition to 
being recorded only once (chiefly because only one plate was available for the 
relevant region), are marked as of questionable reality. But many of them can be 
identified with some confidence as weaker members of multiplets known to be 
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present. In all, some 170 lines have been identified in HD 168476 and no strong 
lines remain unidentified. 

The list of measured wave-lengths and identifications is not published here as it 
is hoped that eventually the star may be studied with still higher dispersion. Buta 
list has been deposited with the Society, and duplicated copies will be issued on 
request to the author. 

Identifications were made with the help of Mrs Moore-Sitterly’s Revised 
Multiplet Table in the usual manner. The representation of the various elements 
is as follows: 


H No Balmer lines are recorded, not even on spectra widened more than 
usual which would bring out broad shallow features. It can be stated 
that if Hy had a normal profile with equivalent width 0-15 A it should 
have been measured on at least one plate; the limit might be put even 
lower at about 0-10 A. The upper limit for H« is much larger owing to 
the lower dispersion in the red. 

H, Although it is excessively unlikely that hydrogen could exist in 
molecular form at the temperature of HD 168476, this question was 
investigated by searching for lines of H, of intensity 10 in the list by 
Gale, Monk and Lee. Only one coincidence occurred with a faint 
unidentified line in HD 168476 and this must be attributed to chance. 

Hel All accessible lines in the R.M.T. except the forbidden appear. The 
diffuse singlets are very strong relative to the triplets. The diffuse series 
21P°~n'D appears strongly in the ultra-violet and probably extends 
beyond the laboratory limit of 3756(m=14). Details of identifications of 
this line and the next two with extrapolated wave-lengths are as follows: 


A measured Int. Suggested identification (with blends) 

3756°80 Bd.NN He 66 56-10, Fe 11 154 55°56, (Ca 11 8 55-6, ?Cr 11 20 54°59) 

3745°92 4 He p 45°8, V 11 15 45°81, ?Fe15 45°56 

373098 = 4 He p 37°3 Ca 11 3 36°90, ?Fe 15 37-13 
The line 3756 appears as a remarkably diffuse shallow feature, attributable to 
blending to a large extent. However, it seems possible that the diffuseness is 
partly connected with the fact that the upper state 14 1D lies only 0-02 cm“ from 
14°D; the upper state may thus be broadened by perturbation. The corre- 
sponding line in the diffuse triplet series is at 3487 A and lies outside the scope 
of the present investigation. 


He Not found. 

C1 Not found, but no good lines are accessible. 

Ci Well represented by strong lines. Components of the following 
R.M.T. multiplet numbers are present: 2, 4, 5, 6, (13, 15), 16, ?17, 
(18, 22), 28, 33, 36, ?37, (39, 40), 245. Numbers in brackets, here and 
in the following, refer to possible contributions to blends. 


C ill Probably absent. 4650 is suspected as an extremely weak feature on 
one plate. 

NI Present; multiplets 10, ?16, 21. The blend at 6483 is moderately 
strong; this feature has been found strong in the hydrogen-weak star 
v Sgr. 

Nu Present; multiplets 5, 6, (8), 12, 15, (17), ?19, (48), ?(58, 59). 


Nu Not found. 
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O1, u, 01 Not found. A strong broad feature is measured at 4075-76 but this is 


Vu 
Cr il 


Mn 11 


attributed entirely to the C m blend of multiplet 36. This is supported 
by the measured wave-length and by the fact that none of the other 
expected O 11 lines (e.g. 4649, 4414, etc.) are found.* It would be of 
interest to search for the strong O 1 triplet at 7774 A, since the level of 
ionization in HD 168476 is probably near the I.P. of O 1. 

Present; multiplets 1, 3, 5, 6, ?9, ?19. Although the lines do not 
appear strongly, on account of the low dispersion in the red, they must 
be enormously strong relative to other B stars. Several Ne lines are 
assigned the same visual intensity (3) as He 5876; the latter line appears 
strongly in the Michigan tracing of 8 Ori, while the Ne lines appear on 
the margin of visibility. These lines of Ne I are also strongly recorded 
by Greenstein in v Sgr, the first star in which they were astrophysically 
observed. 

Not found. 

Presence in stellar or interstellar form is masked by C 11. 

Not found, although 3838 might contribute weakly to a blend. 

Present; multiplets 4, ?5, ?9, 10, ?18, (19), 20, (21), ?25, 26. 

Present; multiplets (1), ?45, ?46,57. The evidence rests mainly on 
the line at 6201 A. 

The line 3905 is very doubtfully present as a contributor to a blend. 
Present; multiplets 1, 2, 3, (4), 5. 

May be represented by multiplet 2, 4552 being blended with N 11, 4567 
appearing on the margin of visibility. 

Possibly represented weakly by members of multiplet 6 but 5253 is not 
found. Higher dispersion in the green is required to settle the 
question. 

Not found. 

Well represented by numerous faint lines; multiplets 6, ?8, 9, 11, 14, 
19, (26), ?29, 736, 40, 243, 44, 45, 49, 59- 

A few lines of multiplet 2 appear to be present. 

Not found. 

Not found. 

The resonance line 4226 is not found. 

H and K strong; multiplets 3 and 8 also found in the ultra-violet but 
blended. 

The few accessible lines have E.P. 30 volts and are probably absent. 
Possibly represented weakly by the lines 4246 and 4314. 

Present; multiplets 13, 19, 20, 34, 41, (82), ?104, ?105. ‘The presence 
of this ion of a relatively rare element is of special interest in view of the 
fact that its ionization potential is practically identical with that of H 
and of O1, which are not found. 

Probably represented by the multiplets ?10 (15). 

4558 possibly appears weakly on one plate, and two lines of multiplet 
20 may contribute to blends. 

Not found; the best lines lie beyond the ultra-violet limit. 


* It appears probable that the same argument applies to the feature at 4075 in HD 124448 
which Popper attributes entirely to O 1. None of the other expected O 1 lines are measured by 
Dr Wesselink in his Radcliffe spectra of HD 124448. 
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Fel May be represented by the strongest lines in multiplets 4, 5; there are 
no other satisfactory identifications for lines at 3719-69 (int. ?4n), and 
3860-23 (?2n); other lines may contribute to blends. 

Fe Well represented by fairly weak lines; multiplets 27, 28, 37, 38, ?40, 
42, 48, 49, 154, ?172, 173, PIgO. 

Fe 11 4419 is a weak line, and 4396 may contribute to a blend. Some lines 
of multiplet 118 (E.P. = 20-5) also seem to be present. 

Nill 4067 and 3849 may be weakly present. 

Other elements; not found. 


In general there appears to be fair correlation between the estimated intensities 
in HD 168476 and Struve’s* estimates in the supergiant 67 Oph (B5 HD; B8s 
Vict.; B5 Ib MK) except for H and On. In Fig. 1 the visually estimated 
intensities in these two stars are plotted against each other. Note that the groups 
of lines with widely different excitation potentials have a fairly random scatter about 
a straight line. Note also that the estimated intensities in HD 168476 run 
systematically higher than for 67 Oph despite the much lower dispersion used for 
the former star; (Struve used in part Mt. Wilson Coudé spectra). 
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Fic. 1.—Comparison of visual intensities in HD 168476 and 67 Oph. 


@ Feu,Can E.P. 0-3 
Meu, Siu E.P. 8-10 

4 Nu,Su,Cu E.P. 14-19 
Het E.P. 20-21 


Arrows to the left indicate lines not observed in 67 Oph. The Balmer lines would appear 
outside the figure below the right-hand corner. 


3. Measures of equivalent widths.—A few selected lines have been measured 
for equivalent width, the plates having been calibrated with a tube sensitometer. 
These measures must be regarded as highly provisional and it is hoped to replace 
them eventually with more extensive ones of higher dispersion. ‘The accom- 
panying table gives the results, which depend in most cases on one plate only. 


* O. Struve, Ap. F., 74, 230, 1931. 
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TaBLe I 
Provisional equivalent widths in HD 168476 (angstroms) 

Line E.W. Line E.W. 
Si tr 3853 0°33 He 1 3964 0°74 
Sim 3856 0°59 N Il 3995 0°30 
Sir 3862 0°48 He I 4120 0°54 
He 1 3927 1°2 Sill 4130 056 
Ca II 3933 0°76 Hel 4143 21 
Ca 11 3968 062 Cir 4267 0:92 


The helium lines and C 11 4267 are stronger than in any star in Williams’ list.* 
The lines with widths greater than 1 A have wings, perhaps due to Stark effect. 
All lines have sharp cores. 

The ultra-violet Si 11 triplet lines have theoretical relative intensities 1: 9: 5 
and since the strongest, 3856, is practically equal to Cat 3968, they can be 
combined with H and K to give a rough idea of the shape of the curve of growth in 
the range 0-3 to 0-8 A, with a range in Nf of relative intensity 18:1. The mean 
slope (A log W/A : Alog Nf) is about 0-3. 

4. Excitation and ionization.—The lower excitation potentials of lines appearing 
strongly in HD 168476 range from 0-o (H and K) to 21-1 (He 1) and 24-3 (C 1). 
The ionization potentials of well-established ions range from 11-8 (Ca I1) to 29°5 
(N 11); but outside this range there is doubtful evidence for Fe 1 (7-86 volts), 

?Mg1(7-61) and Si1 (8-11), and at the high level, Fe 111 (30-5) and Si 111 (33-3). 

The excitation and ionization potentials represented are illustrated in Fig. 2, 
the joining lines or dashes giving a rough indication of the intensities of the 
corresponding lines. There are two points requiring special mention concerning 
this figure. First, the wide range in excitation and ionization, which would 
probably be still better indicated were higher dispersion available, has a parallel 
in the wide range exhibited by the star v Sgr and emphasized by Greenstein and 
Adams.t Secondly, the range in E.P. easily includes those of the Balmer lines 
and aiso includes near the upper limit those of strong lines of O 1 which have not 
been found. Moreover, the ionization potentials of H and O 1 are both practically 
identical with that of the relatively rare ion Ti 1 which is observed. As with 
v Sgr, but to an exaggerated degree, there appears to be no escape from the 
conclusion that hydrogen really is extremely deficient in the atmosphere of 
HD 168476, being far less abundant than helium, probably less than carbon and 
neon, and that oxygen is less abundant than carbon. It is unfortunate that the 
present observations do not cover the ultra-violet series limit of helium where 
Popper has observed a depression in the spectrum of the similar star HD 124448, 
but there can be little doubt that the same effect could be observed in HD 168476. 

5. Consequences of hydrogen deficiency.—The cosmic abundance of hydrogen 
has so many important consequences that one would expect to see many marked 
peculiarities in a hydrogen-weak atmosphere. The effects have been considered, 
for instance, by H. N. Russell § when he originally put forward the evidence for the 
great abundance of hydrogen in his classical papers. 

The first effect is due to the fact that in the atmosphere of a normal late B type 
star hydrogen provides essentially all the electrons and is thus the most important 

* E. G. Williams, Ap. F., 83, 279, 1936. 
+ J. L. Greenstein and W. S. Adams, Ap. 7., 106, 339, 1947; see esp. p. 349. 


t D. Popper, P.A.S.P., §9, 320, 1947. 
§H. N. Russell, Ap. 7., 78, 239, 1934 (esp. last two paragraphs p. 297). 
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source of opacity. If the hydrogen were removed, all lines due to other elements 
would be strengthened, due to the increased transparency. But the lowering 
of the electron pressure would increase the state of ionization and produce effects 
normally attributed to a supergiant atmosphere. 
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Fic. 2.—Representation of excition and ionization potentials in HD 168476. 


Strong lines are represented by continuous lines joining the I.P. and lower E.P. 
Intermediate lines are represented by dashes and weak lines by dots. 


In HD 168476 we observe such effects. All lines in such supergiants as 67 Oph 
and £ Ori, with the exception of H and O 11, appear to be strengthened. Increased 
transparency of the atmosphere compared even with a supergiant is definitely 
indicated. 

Secondly, there is the question of the hydrogen convection zone. If hydrogen 
should be less abundant than say carbon, then presumably helium would replace 
hydrogen as the effective agent in setting up a convection zone. At the temper- 
ature of HD 168476 this would presumably lie chiefly below the photosphere, 
just as does the hydrogen convection zone inthe Sun. However, Eddington* has 
shown how the solar convection zone can extend itself into much higher layers 
than the limits given by Unsdéld’s original calculations. It seems reasonable to 
suppose that the helium convection zone in HD 168476 extends up into the 
observable portion of the atmosphere. This, as well as the transparency of the 


atmosphere, may contribute to the large range in excitation and ionization in the 
observed .spectrum. 


* A. S. Eddington, M.N., 101, 177, 1941; 102, 154, 1942. 
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Thirdly, a deficiency of hydrogen might have a marked effect upon the damping 
portion of the curve of growth. Collisional damping, which may vary from line 
to line, is dependent in a normal atmosphere on the density of hydrogen and of 
electrons. It is hoped eventually to construct a curve of growth for HD 168476 
from plates of high dispersion. But what evidence there is at present available 
suggests that the damping portion of the curve of growth is not reached by lines as 
strong as H and K (see Section 3). 

6. The transparency of the atmosphere.—We can make a rough estimate of the 
increased transparency by means of the equivalent widths in TableI. HD 168476 
is matched much better by 67 Oph than by £ Ori, but unfortunately no equivalent 
widths are available for 67 Oph. Instead, the following indirect approach is 
adopted. Struve’s visual estimates in 8 Ori were calibrated by measuring the 
equivalent widths on the Michigan Atlas tracings of thisstar. It was then assumed 
that Struve’s visual estimates for 67 Oph followed the same calibration. The 
widths derived in this way for 67 Oph (W,,) were then compared with the widths 
directly measured in HD 168476(W*). We thus find: 


Line 3933 3964 3995 4120 4130 4143 4267 
log W*/W,, +°40 +°3I +°39 +°25 +°34 (+°62) +-28 


Excluding the strongly winged line 4143 we find as a mean value 


log W*/ We =-+ 0°33 + 0°02. 


Thus the E.W.’s in HD 168476 are about twice those of corresponding lines in 
67 Oph. This ratio applies to the portion of the curve of growth covered by the 


Siu and H and K lines, whose slope is about 0-3. Hence it would need a reduction 
by a factor 10 in the number of atoms above the photosphere of HD 168476 to 
bring the equivalent widths in this star to equality with those in the supergiant 
67 Oph. ‘This could be taken as a measure of the relative transparencies of the 
atmospheres of the two stars, provided the parameters of the curves of growth 
(turbulence, temperature, etc.), are effectively the same. 

7. The radial velocity.—As announced in the discovery note, the velocity of 
HD 168476 is high. From 21 plates we find a velocity relative to the Sun of 


— 165-0 + 0-8 (m.e.). 


Some of the earlier plates suggested possible variations in velocity over intervals 
of some months, but from the present evidence the velocity should probably be 
regarded as constant. The range in velocity for 4 “b” plates (29 A/mm) is 
7:7km/s, for 17 “‘c” plates (49 A/mm) 18-7km/s. The differences between 
individual plates may arise in part from the influence of faint blending lines, 
especially if the lines vary slightly in intensity—a circumstance for which there is 
some evidence. If the differences are real they cannot be reconciled with any 
period in excess of a day. There is no-evidence for systematic residuals for 
different lines which could not be attributed to blending. In particular there is 
no evidence for duplicity of the star. 

It is of particular interest to note that the other two hydrogen-poor stars, 
HD 124448 and 160641, also have high velocities of —65 and about + 100km/s 
respectively. The positive velocity of HD 160641 is important as indicating that 
in these three stars we are not dealing with peculiar expanding atmospheres. 


a 
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They must be regarded as high-velocity objects. It has already been pointed out 
that the proper motions of HD 168476 and 124448 are both small and apparently 
normal for B type stars of their apparent magnitude. 

8. The luminosity and colour.—It is exceedingly difficult to draw any secure 
conclusions regarding the luminosity of HD 168476. Popper has suggested an 
absolute magnitude of about — 2:5 to — 4 for HD 124448, partly from spectroscopic 
criteria and partly from the radial velocity which could be attributed to the effect 
of galactic rotation. However, since these stars are apparently high-velocity 
objects, the latter argument is quite inconclusive; and it seems dangerous to 
apply spectroscopic criteria to such peculiar spectra. Hydrogen deficiency will, 
as mentioned above, tend to produce the characteristics of a supergiant. 

However, it does not appear that the absolute magnitude can be fainter than — 1. 
A tangential velocity of 100 km/s at 1400 parsecs would yield an annual proper 
motion of 0”-015 of arc, which should be perceptible. The proper motion data 
suggest absolute magnitudes of the order of — 2 or brighter. 

Another line of approach lies in the intensity of interstellar lines, but here the 
high galactic latitude (— 20°) of HD 168476 is a drawback. A line of intensity 2 is 
recorded at 3970°51 (allowing for a velocity of —165 km/s). No satisfactory 
identification is possible, the contributions of C 11 3970-20 and S 11 3970-69 being 
probably negligible. It is more reasonable to attribute the line to interstellar 
Cai witha velocity of -10 km/s. Unfortunately the corresponding K component 
is lost entirely in the strong stellar line He 3935-9. ‘Thus the normal circumstance 
that interstellar K can be measured more accurately than interstellar H is reversed. 

The equivalent width of the ‘‘ 3970 line”’ has been measured on one plate as 
0:20 A. Assuming it to be entirely due to interstellar Ca 1, the width of inter- 
stellar K should be between 0-30 and 0-35 A. Miinch* gives a distance of 550 
parsecs for the star HD 84567 at b= +18°, with equivalent width 0-14 for the K line. 

It is probably safe to conclude that if the “ 3970 line” in HD 168476 is correctly 
identified with interstellar calcium, then the star must lie at a distance of the order 
of a kiloparsec, or more; this gives a distance of 340 parsecs from the galactic 
plane, which would take it outside most of the gaseous substratum. An absolute 
magnitude M = —1 or brighter is indicated. 

Dr Wesselink has very kindly undertaken at my request photoelectric 
measures of the star. He has derived a colour index on the Cape 1953 S system of 
—o-36. Using the linear relationship in Cape Mimeogram No. 3, this corresponds 
to B— V= —0-05 on the Johnson system. The colour matches that of a normal 
late B type star. 

9. Acknowledgments.—I am much indebted to Dr A. J. Wesselink, Dr M. W. 
Feast and Mr A. Morrisby for obtaining six p.ates of HD 168476 on closely 
adjacent nights in 1953 for the purpose of studying possible short-period variations 
in radial velocity ; also to Dr Wesselink for undertaking the photoelectric measures 
of colours of the star. The author is responsible for the other plates and all 
spectroscopic measures and reductions, 

Radcliffe Observatory, 

Pretoria: 


1953 October 22. 
* G, Miinch, P.A.S.P., 64, 312, 1952. 








PHOTOMETRIC EFFECTS OF REFLECTION IN CLOSE BINARY 
SYSTEMS 
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Summary 


The aim of the present investigation has been to study the mutual 
illumination of two luminous spheres, reflecting light in accordance with 
Lambert’s law, by taking into account the finite angular size of the light 
source. 

After an expository section containing an outline of our problem 
(Section 1), exact expressions will be set up which govern the amount of 
light of one sphere reflected by the other for any geometrical configuration 
(Section 2). In Section 3 we shall specify the flux of radiation incident 
upon the fully illuminated as well as penumbral zones of the reflecting 
component; while in the subsequent Section 4 the exact results of Section 3 
will be simplified to the order of accuracy to which two fluid components 
of a close binary system can be regarded as spheres. 

Section 5 will then contain a derivation of the law governing the 
variation of the amount of reflected light with the phase; while the concluding 
Section 6 will be devoted to a brief analysis of this law, as well as to an 
outline of its generalization to the case of distorted binary systems, or to 
systems exhibiting eccentric orbits. 





1. Introduction.—It is inevitable in close binary systems that a part of the 
radiation of either component, sent out in every direction, will be intercepted by 
the surface of its mate where it will be absorbed and re-emitted (or scattered). 
This fact is bound to cause the hemisphere of each star exposed to incident light 
to be brighter than that turned away ; and as the components revolve, the amount of 
light* ‘‘ reflected” in the direction of the line of sight should vary with the phase. 
This effect has, in turn, long been known to give rise to by far the most intricate 
part of the theory of eclipsing variables. Even for a plane exposed to parallel light 
the law of angular distribution of re-emitted radiation is not yet known exactly (x) ; 
and in close binary systems the difficulty is greatly increased by the large fractional 
size of their components. In point of fact, the problem of mutual illumination of 
two finite spheres alone becomes geometrically so complicated that no other 
aspect of the interpretation of light curves of eclipsing binary systems has con- 
fronted us with greater difficulties. 

In the simplest case of the incident light forming a parallel beam (which 
would correspond to an infinite distance of the illuminating source) and re-emitted 
in accordance with Lambert’s law, the solution of our problem is, of course, well 
known and had been investigated in connection with planetary phases many years 
ago(2). A first attempt at taking account of the convergence of the incident beam 
(which must be considered if the light source is placed at a finite distance from the 


* Under this term we shall hereafter understand “ integrated radiation in all wave-lengths ’’. 
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reflecting star) was due to Takeda (3), who regarded the secondary component as a 
non-dimensional light-point illuminating the primary which (as seen from the light 
source) exposed to the incident light an apparent disk of finite angular size. 
Takeda carried through an approximate solution of this problem to quantities of 
the fourth order in the fractional radius of the primary (reflecting) component 
(i.e. to the order of accuracy to which a distortion of this star due to its own axial 
rotation as well as to the tidal action of its companion can be ignored). In doing 
so Takeda committed, however, a mistake of extending the cone of incident light 
beyond the horizon of visibility of the illuminating light source at the primary 
star—an oversight which vitiated his fourth-order term and which was eventually 
corrected by Sen (4). 

A reconnaissance of the effects of a finite angular size of the illuminating compo- 
nent on the amount of reflected light was first undertaken by Russell (5), by an 
approximate method which was essentially incapable of generalization and which 
merely permitted a surmise that the contribution of the penumbral zone (where 
the apparent disk of the illuminating component would appear partly set for an 
observer on the primary star) is likely to amount to less than one per cent of the 
total light reflected. Two years later a discussion of the same problem was 
resumed by Matukuma (6). Like Russell, Matukuma did not aim at a derivation 
of the phase-law of the reflected light, but limited himself to ascertaining the 
proportion of the total incident light to that reflected in all directions (expressed by 
what Matukuma called the “criterion formula”). However, in concluding his 
paper already quoted, Matukuma mentioned that... “‘recently the writer had 
taken up this problem” (i.e. of the derivation of phase-law) ‘‘and has almost 
completed the necessary calculations. But some minute point, which is obscure 
in a physical sense, makes him postpone publishing his results”. Unfortunately, 
shortly after writing these words, Matukuma died—apparently without being able 
to clear up his doubts; and an examination of his “ criterion formula ”’ reveals that 
its derivation was in error*, invalidating his result. 

These reasons led the present writer to open up the problem de novo, and to 
work out a detailed theory of the mutual illumination of two finite spheres reflecting 
in accordance with Lambert’s law. Such a theory should, in particular, reveal a 
variation of the amount of reflected light with the phase, taking account of the 
finite angular size of the illuminating component. _ Its application to astrophysical 
cases will be limited, on one hand, by the extent to which the actual reflecting 
properties of stellar outer layers depart from Lambert’s law; and, on the other, by 
the assumed spherical shape of both components. Approximate solutions of the 
problem of radiative transfer of reflected light in stellar atmospheres (1) reveal 
indeed the existence of departures from Lambert’s law which should be taken 
into account in more refined work. Such refinements could, however, be 
carried through only at the expense of detailed numerical integrations. On 
the other hand, the assumption that the components possess spherical shape, 
implying the neglect of first-order rotational and tidal distortion, will restrict the 
correctness of our work to quantities of the fourth order in the fractional radii 
of the components. A general theory developed below could, in principle, be 
extended to distorted stars without essential modification—at the expense of 


* The integral for N, in Matukuma’s equation (13)’ has been evaluated incorrectly—its true 
value being (27,/37p) sin? rather than 2(r,/7p) sin as given by Matukuma—and this slip 
avalanched through the rest of his paper. 
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nothing else but additional algebraic work. At present we propose, however, to 
limit ourselves to a consideration of both components as spheres, and to postpone 
complications arising from their distortion to subsequent communications. 

2. Geometry of the problem.—In order to fix our ideas, let us adopt a rectangular 
frame of reference XYZ, with origin O, at the centre of gravity of a star to which we 
shall hereafter refer as the primary component, and oriented so that the X-axis 
is constantly in the direction of the secondary (illuminating) component, while the 
Z-axis is perpendicular to the orbital plane. Let, furthermore, 7, 7, ¢ be the 
spherical polar coordinates in this system, related with the rectangular ones by 
means of the relations 

x=r cos ¢ sin 7, 
y=r sin ¢ sin 7, (1) 
=F COS 0, 


such that 7 denotes the polar distance of an arbitrary point P on the surface of the 
primary (see Fig. 1); and 4, its longitude measured from the “ principal meridian” 
in the XY-plane. 


Fic. 1. 


As the components revolve, this system of axes will evidently rotate with respect 
toa fixed line of sight. Let the direction of the latter be described by the direction 
cosines 


1, =cos ¢& sin 4, 
My = sin ys sin 2, (2) 
Ny = COS 1, 
where ys denotes the true anomaly of the secondary component in its relative orbit 
(measured from the moment of inferior conjunction), and i the inclination of the 
orbital plane to the celestial sphere. The angle « between this line and the X-axis 
joining the centres of both components (whose direction cosines are, therefore, 
I, 0, 0) will evidently be given by 
cos «=/,=cos y¢ sin 1. ~ (3) 
Moreover, the amount # of the secondary’s light “reflected” from the primary 
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component in the direction of the line of sight should then depend on the true 
anomaly ; only through e, and be generally expressible as 


¢=| F cos y ds, (4) 
c 


where / denotes the intensity distribution of reflected light; y, the angle of 
foreshortening ; and ds,, the corresponding surface element of the reflecting star ; 
while the domain C of integration is to be extended over the crescent visible to a 
distant observer at any particular phase. 

Consistent with the foregoing definitions, the angle y of foreshortening can be 
solved for from the spherical triangle NPE (Fig. 1) in the form 


cos y=sin 7 cos (e—¢), (5) 
ds, =r, sin 7 dn dd, (6) 
where r, denotes the (constant) radius of the primary component. 


As to the limits of integration in both variables let us, for the sake of brevity, 
introduce the integral operator K(x) defined by 


cos! « pa—sin— (x sec ¢) 
| * ride cos y ds,=K(x){...}, (7) 


sin~" (2 sec $) 


while 


and consider the illumination of one finite sphere by another. The fully illumin- 
ated portion of the primary’s surface (at which the whole apparent disk of the 
secondary component of radius r, will be visible to an observer situated at P) will 
evidently be subtended by the inner tangent cone to both spheres, within which 


m <0 <7—-My (8) 
where 

‘ mn4+Tr 

sin 7, = ~~ :, (9) 





- R=O,0, denoting the separation of centres of the two components. 

For 7=7, (or 7—7,) the lower limb of the apparent disk of the secondary 
component will touch the horizon of the observer at P; and for 7 outside the 
inequality (8) this disk will gradually sink below the horizon of the observer until, 
for 7 = 7 (or 7 — 2) such that 


sin 72= a ’ (10) 





it disappears completely. Therefore, as long as 
1 <€<a—M, (11) 


the fully illuminated zone will eventually be delimited by an application of the 
integral operator K(sin 7,); while the penumbral zone (illuminated by only a part 
of the apparent disk of the secondary) should be delimited by an application of the 
operator K(sin »,)— K(sin 7,). 

The foregoing statements are, to be sure, exact as only as long the phase angle e 
continues to lie between the limits of the foregoing inequality (11); for when 


— S€ Hy (12) 
the fully illuminated zone becomes completely exposed to the external observer, 
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and the lower limit « — 90° of integration with respect to 7 in the operator K(sin 7,) 
should be replaced by —cos~*(sin 7,).* When, on the other hand, 


m7, << 7+, (13) 


the whole fully illuminated zone becomes invisible and the integral K(sin ,){... } 
then vanishes. Similarly, when 


— 2 S€ Sy (14) 
the limb of the penumbral zone becomes a circle, and the lower limit «—9g0° 
of integration with respect to 7 in the operator K(sin 7.) should be replaced by 
—cos(siny.). When, ultimately, 


7—Ng SE <7+ He (15) 


and r,>r,, the whole penumbral zone will be lost out of sight of our external 
observer, and the integral K(sin 7) { . . . } will consequently vanish. Should, on 
the other hand, the illuminating component be the larger of the two the penumbral 
zone will not disappear for e constrained by (15) but will continue to surround the 
primary component with an illuminated ring—obtained again by replacing the 
lower limit « — 90° with respect to 7 in K(sin y2) by — cos~*(sin n,) and remembering 
that, for 7; <7rg, sin n, becomes a negative quantity. 

It may be added that the foregoing equations, as they stand, express the fraction 
of light of the secondary component reflected by the primary star. The fraction 
of the primary’s light reflected from the secondary at any phase can likewise be 
deduced from the foregoing equations, provided only that the values of r, and r, 
are interchanged, and that the phase angle « is consistently replaced by 7—«. The 
total amount of light reflected by a close binary system will then be equal to the 
sum S(r,, f2; €)+L(rq,7,; ™—€); and it is this sum which must be taken out 
of the observed light changes before the light curve of the system is analysed 
for photometric elements in the usual manner. 

3. The incident flux.—In order to be able to evaluate the integral (4) as a 
function of «, we must specify the intensity 4 (¢, 7) of radiation reflected at any 
point P of the illuminated portion of the primary star. This intensity can indeed 
be established provided that we can ascertain (i) the flux of radiation of the 
secondary component incident at P ; and (ii) its transfer (by absorption-re-emission, 
or scattering) outwards consistent with the requirement that the heat-albedo of the 
stellar surface must necessarily be unity. This second part of our problem, first 
considered by Milne (7), has recently become the subject of renewed investigation 
(x); but its proper solutions are so far only approximately known. In what 
follows we shall, therefore, by-pass this physical difficulty by assuming that the 
reflected radiation is uniformly distributed (which is always true to a first approxi- 
mation), and that its intensity varies in accordance with Lambert’s cosine law 


JF ($,) = S cosa, (16) 


where 7S denotes the flux, per unit area normal to it, originating in the secondary 
component; and « is the angle of incidence (cf. Fig. 2). In the present section 
our problem will, therefore, be limited to expressing both S and « in terms of our 
angular variables ¢, 7. 


* This case is, however, largely academic; for if | «| < ||, the secondary component is 
going to eclipse the primary star; and if r,>0, a part of the reflected light will be lost by eclipse 
—a situation whose analysis is outside the scope of the present paper. 
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In order to do so, consider a new rectangular system of axes, with origin 
at P (see Fig. 2), whose Z-axis coincides with O,P=p, and the XZ-plane coincides 
with that of the triangle O,PQ (i.e. with the plane of the diagram); the Y-axis 
being perpendicular to it. The direction cosines of an arbitrary radius-vector 
from O, through P then are 

COs a, SiN a, 0; 


while the direction cosines of PQ’=p’ in the same system become 
cos 8, sin 8 cos w, sin B sin w, 


where f denotes the angle O,PQ between the vectors p and p’, and w= 7 QP’Q’ 
(see Fig. 2). 





Fic. 2. 


If so, the angle { between p’ and the surface normal at P will evidently be 
specified by 
cos (=cos « cos 8+sin « sin B cos w, (17) 
where, from AO,PQ, 
sin B= > sin 6’ and cos B= meee (18) 
while (from the same triangle) it follows that 


2_ 7,2 9'2 os0—r, 
os (a f— 3 ot. (19) 
2p le Pp 





Therefore, 


p’ =p? +1,2—2rep cos 8, (20) 
while (from AO,O,P) 


p? = R? + r;? aad 2r,R cos Y> (21) 
where y’= 7 0,0,P, defined by 
cos y’ =sin 7 cos ¢; (22) 
and, ultimately, 
Reosy’ —rz 
p oe 


cos «= 


(23) 
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The last three equations (21)-(23) define « as a function of ¢ and 7; while 
the flux ~S incident at P can be ascertained as follows. Let the distribution of 
brightness 1(@’) over the apparent disk of the secondary component vary in 
accordance with the well-known cosine law of limb-darkening as 


1(6’) = I(o){1 —u+u cos 6}, (24) 


where J(0) denotes the intensity at the apparent centre of the illuminating star 
and u, its coefficient of darkening. The intensity .4(¢, 7) of light reflected at P 
can then be expressed as 


I($, 0) =1(0){(1 —u)J, + ud 9}, (25) 
# wo [EE ae 


where 
2» (26) 


the surface element ds, of a (spherical) secondary component being, in turn, 
expressible as 


ds,=r,? sin 0 d0 dw=r,(p'/p) dp’ dw; (27) 
while cos { and p’ continue to be defined by (17) and (20). 

As long as the whole apparent disk of the secondary component remains 
visible to the observer at P, the limits of integration in w extend from 0 to 27. 
If so, however, the integration with respect to w clearly annihilates the second 
part of the right-hand side of (26) arising from the second term on the right-hand 
side of the equation (17) for cos ¢, leaving us with 


2" cosa cosB p?—r,? 
2p 


Pi=P—T2 
neem, | 


and cos « as well as cos 8 have already been defined in terms of our variables of 
integration. Integrating, we thus establish that 





ie p”? 
J, = ie dp’ dw, (28) 


where 


cos « " (p? — r,”)? 
7 aoe 


<P ap = (30) 
—.——_ dp’ = cosa; 0 
2p Pi Pp 4 p? ’ 3 
while, similarly, 





® (2* cosa cos B (p?—r,2—p’")2 rap’ , , 
Jo= | | = { 1 — dp’ dw 
ale 2r ep mp 


2 
= ye Cos a. (31) 
Inserting (30) and (31) in (25) we find, eventually, that 


L, cos x 


JI(¢, 0) = ’ (32) 
where 

Ly = mrg4l(0){I — $u} (33) 
tepresents the apparent luminosity of the secondary component. At any point 
of the fully illuminated zone the flux incident upon P proves to be proportional 
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to the total light of the secondary component, and inversely proportional to the 
square of the distance of its centre from P; so that, by (32) and (21)-(23), 


ae L, Ru - ") 
JI(%, 0)= ry (R? +7,?— 27,Ru)>?’ (34) 





where we have abbreviated 
#=cos y’ =sin 7 cos ¢. (35) 
In the penumbral zone the situation becomes complicated by virtue of the fact 
that the secondary component now undergoes eclipse by the horizon of the 
observer (which acts as a straight occulting edge), and only the light of the 


uneclipsed portion will continue to illuminate the surface element at P. In order 
to evaluate the intensity of light reflected from any point of the penumbral zone, 


let us write again 
I'($, 0) =1(0){(1—u) Jy + uJ '}, (36) 
where now 


‘08 (r/p) ‘(sin A csc 0) 3 
nd,’ 2) (.. (p’)-? cos { cos" 6’ sin @.d0 dw, (37) 
A 0 


the limits of integration being now adjusted so as to extend over the visible 
portion of the secondary’s disk visible to the observer at P. The (plane) angle A, 
which will hereafter play the role of a ‘‘ geometrical depth of the eclipse”’ is, 
in turn, defined by the triangle HPO, (cf. Fig. 2) as 

A=a—cos-\(p/rg) cos «, (38) 
and normalized so as to vary between +cos~!(r,/p) as the apparent disk of the 
secondary component continues to set behind the horizon of the observer. 
We shall, moreover, reckon A so that its lower limit corresponds to the beginning 
of the eclipse (i.e. to the observer situated at the boundary between the 
fully illuminated and penumbral zones). 

The integration of (37) with respect to w offers again but little difficulty ; 

for inserting cos { from (17) we readily establish that 

J, =M,,cos«+N, sin a, (39) 
where 


cosX7/) Cos i 


cos“! (=) cos" 6’ sin 6 dé (40) 


7M, = ar? | 


72 
A Pp 


P cos (rs/P) sin B . 9 2 \\L/2 ie 
aN, =2r, f <n 6 — sin? A)*? cos" 6’ dé, (41) 


respectively. Unlike the case of flux in the fully illuminated zone, however, 
the second integration of the foregoing expression with respect to @ can no. 
longer be performed in a closed form without appeal to certain simplifications 
which we shall presently proceed to discuss. 

4. Solution of the equations.—In the two foregoing sections we have formulated 
rigorously the problem of mutual illumination of two spheres reflecting in 
accordance with Lambert’s law. The expressions for the flux and limits of 
integration could, in fact, be used to study the reflection of light from two spheres 
in actual contact, or from a sphere illuminated by an infinite wall—in both of 
which cases there would be no fully illuminated zone, but only a penumbral zone 
(extending, in the latter case, over the entire surface of the reflecting sphere). 
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Such extreme cases would, however, be of little or no astrophysical interest, 
since two stars brought so close together could not possibly retain spherical 
shape, and their mutual distortion would render the geometry of our preceding 
sections inexact. As was mentioned in the introductory paragraphs of the paper, 
it is our intention to carry through the solution of our problem only to the order 
of accuracy to which both components of our binary system can be regarded as 
spherical (i.e. to quantities of the fourth order in their fractional radii). 
In doing so, we surmise, in particular, that the amount of light reflected from 
the penumbral zone will represent a small quantity whose square and higher 
powers become negligible within the scheme of our approximation. 

If, however, only terms of the lowest order are to be retained, we may 
evidently set 


sin B=(r,/p) sin 0, | 
cos BI, 
(42) 


cos 6’=cos 6, 

p' =p, | 

in virtue of which the flux of radiation incident on the penumbral zone may be 
approximated by 


2 2 i 
M,= = r cos~? e ;) sin 6 cos 6 dé 
A 


r* {; A+sinA cos *} 


212 7 


2 pa/2 
= 25 | (sin? @—sin® A)? sin 0 cos 6 dé 
A 

3 

= 5 cos? A 

3p 

for uniformly bright disks, while 
2 px/2 
M,=—> | cos“? ee ;) sin 0 cos? 6 dé 
mp* J, sin @ 


r,* :{ (x+ cos* *)\ 
= +,41-sinA 
3p” 2 
™ i" (sin? 6 —sin® A)¥* sin 8 cos? 6 dé 
r° 
= 33 cos‘ A (46) 
for disks completely darkened at the limb. 


Inserting the foregoing equations in (39) we are then in a position to assert 
that, within the scheme of our approximation, 


2/1 oe 
fea COS 


2ar,* 3 X si 
“ a+ 373 cos* A sin a, (47) 


’ r : cos? X 
Jim 2{2—sin a (24 - ) feos «+ 3 gpa cost A sin a, (48) 
where, to the order of accuracy we are working, 


7 7 
— 3 SA<5- (49) 
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If A= —g0° (corresponding to the boundary between full-light and penumbral 
zone), the foregoing expressions reduce to 


, 7 r ° - 
J, (- 3) = Kees, (50) 
2 
J,’ (- *) = 34008 a, (51) 


in agreement with the results established previously by equations (30) and (31) 
for J, and J, in the fully illuminated zone; whereas if A= + 90° 


1 (Q)—4() = 0 


as they should, since the eclipse of the secondary component has now become 
complete. 

In order to simplify further the foregoing expressions for J, , it should be 
remembered that within the penumbral zone « will be in the neighbourhood 
of go°, thus rendering cos« a small quantity of the order of r./p. If so, 
equation (38) defining A can be replaced by 


A=sin~ (R/r,) cos « (53) 


and, in terms of the order of (r,/R)®, we may approximate sina=1, while 
cos «=r,/R. In consequence, equations (47) and (48) may be rewritten as 





, v® fsin*x+x(1—x*)¥? I 2r,° 
, tefe 
J» = Fa{g(t— aXe +a—2)+ 40-24), (55 


and expressed in terms of the single variable x=sin A which, consistent with (53), 
can be approximated by 
x=sin A=(R/rq) cos «, (56) 
and is by (49) constrained to vary within —I <x <I. 
Moreover, an appeal to equation (23) discloses that, to the order of accuracy 
we are working, 
R R ry 
aoc dalt Saat %, (57) 
2 
where, as before, »=cos y’=cos¢siny. If, furthermore, L, stands again for 
the luminosity of the secondary component—as defined by equation (33) above— 
the equation (36) governing the intensity of light reflected by the penumbral 
zone will, with the aid of the foregoing expressions for J’, », ultimately assume 
the form 


I— 
FG 1) = 9 9G, + FG, (58) 
where, as before, u denotes he piconets s picetori of limb-darkening, and 
IY (dy 0) = ZTE, (6x(sin-t 24+. xV Em 30) +4(2—")94}, (59) 


while 





I$, 0) = - 247 rea (1 — x)°(7x* + 14% + 3), (60) 
respectively. As x can, by means of (56) and (57), be readily expressed as a 
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function of u=cos¢sin7, the requisite explicit forms of %%°(¢, »)—and, 
therefore, of #£’(¢,»)—as functions of our angular variables are thereby 
established. 

5. The phase-law.—A possession of the foregoing results puts us at last in a 
position to ascertain the fraction of light of one component of a close binary 
system from the other at any phase; for this task calls merely for an application 
of the K-operator, as defined by equation (7), to the expressions for 4(¢, ) 
established in the preceding two sections. 

In the fully illuminated zone, K=K (sin ,), while the intensity .4(¢, 7) of 
reflected radiation is known to be given by equation (34). Moreover, an 
expansion of the right-hand side of this latter equation in ascending powers of 
r,/R in terms of the Legendre polynomials P;(u) of ~=sin ¢ cos n, permits us 
to replace the closed form of (34) by 


F (bs 1) = 2,4 Py(u) +22 Pyle) +3 Plu) + -.-$, (6r) 
a7R R R 


correctly to the order of accuracy adopted in Section 4. Hence, the corresponding 
amount of light of the secondary component reflected from the fully illuminated 
portion of the primary star in the direction of the line of sight should be given by 


Lilo 3 )=Kicin n){ Plu) +2 Pylu) +3 Fe Pu) +» SR. (62) 
Performing the actual integration in a closed form we find that 
K (sin ,){P,(#)} = $717{,—s sin® 71} cos ¢ + §r,°F cos? n,, (63) 
K (sin 9,){ P2(4)} = 37,°®, P, (cos €) 
+ 47,?®,(1 + 3 sin? 7) cos* 7, cos € 
+ 47,?F(4 +3 cos* n,) sin 7, 
K (sin ;){P3(#)} =71°@s sin® n, cos? 7, cos € 
+7,°F sin? y, cos* n, — 47,°F°, 
where we have abbreviated 
®, =cos™ (sin 7, csc €), 
®,=cos~! (—sec 7, cos €), 
®, = — cos“ (tan 7, cot e) 
and 
F =(sin? e — sin? 7,)"*. 
Since, however, 
®,= - —sin y, csc e—} sin® y, cos e— ..., 


®,=7—¢€+} sin® y, cot e+ 4 sin* y; (cot? «+3) cote+..., 


7 e ° 
®,=- > +sin 7, cot e+} sin* 7, (cot? «+3) cote+..., 


F=Ssin e—} sin? n, csc e—} sint n, cscke— ... 








= Pact ae 2 a Ce 
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correctly to quantities of the fourth order in sin 7,, it follows that, within the 
scheme of our approximation, 


2 (r71\2 (7 —€) cose + sine 
#(%, rid=laf=(3) { “ 


3 
+ (2) (3 cos* e+ 2 cos «—I) 





(74) 


) 4 sinecos*e 7,?r,” sin € 


for any phase angle « within 
sin? 22? <¢<n—sin we ° (75) 
In the immediate neighbourhood of the inferior conjunction, when 
. a t7, rn+r. , 
—sint+7— <e <sin1——, (76) 
the K-operator (7) must be modified as explained in Section 2, in order to account 
for the fact that the fully illuminated zone should become completely exposed 
to the external observer (though it is bound to be partially masked by eclipses) ; 
and the result of such integration reveals that, during this phase, 
2/2 mtr 
F,(1;, 723 €)= laf 3 + is" + -. } 608 6 (77) 
correctly to the same order of accuracy.* If, ultimately, 
‘ mn+T . +r 
7—sin-} . * <e<a+sin? ++ (78) 


and the fully illuminated zone becomes completely invisible to our observer, 
the foregoing expression for #,(r,, 72; €) is indeed seen to vanish for «=z, and 
to become an ignorable quantity of the fifth order in r/R for all other ¢’s within 


the limits of the inequality (78). 
The fraction #,(r,,72.;«) of the secondary’s light reflected from the 


penumbral zone of the primary star may, in turn, be es into 
3(1 aa ov 4 a #; o> 

3-u 3-u 
where u denotes, as before, the coefficient of limb-darkening of the illuminating 
star, and 


(79) 


LAr, 723 6)= 


£9 ={K (sin 42)—K (sin 1)}{ 4 (¢, 7)}- (80) 
Now the expression (59) for #¥ can be expanded in ascending powers of x to yield 


I°($, 9) = as += + wa (8) 


within an error of less than one per cent for the extreme values attainable by x; 
and equation (57) permits us to convert the foregoing expansion into a series of 
the form 


FXG, 1)= HAE CMY, (82) 


* A part of this light will, of course, be lost by the eclipse. 
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nV a)" + 
ad Fe 12 
r 
+2(2) PIs 
3 \"e 
r,\2 
*) ee 
ie) 


where 





On the other hand, by virtue of (60), 


L, 


D = v2) 3 a 
J weogee £4242 Ga 


. 

2 
_ "he ¥ p/ Ru 
~ aR ois o © (=) 





ee 


Now as the reader can easily verify, 


2r,* ’ 
K(a){u"} = K(o){u"} — oe gly sine+..., 
by virtue of which 


{K(sin 7) — K(sin ;)}{¥""°(¢1n)} 
L u,D 
=2 —_ oe zo . (*) (sin"*+! y, —sin"*+? ,) sin e 


= 2 (FB) ty sine SEO fate) (ony, 
R nwo M+1 lo lo 
where all CY’s for j >4 are identically zero. Inserting the appropriate coefficients 
C,, as given by the foregoing equations (83) and (85), in (87) we find, moreover, 
that the whole summation on the right-hand side of (87) reduces to 7! in the case of 
a uniformly bright secondary component (u = 0), and to 0-3 in the case of complete 
darkening (u=1). In consequence, the requisite amount of light reflected by the 


8 
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penumbral zone of the primary component may, to the order of accuracy we have 
been working, be approximated by 


2ry"s 


2 
He" (1% ; = {zal L,sine+..., (88) 


T1l2 


6 
FL? (11,72 5 = = {eh L,sine+. (89) 


so that the amount of light, sia by the penumbral zone, incident from a 
secondary component arbitrarily darkened at the limb becomes 


12 { mr2)2 (5 +(7—- , 

L3(t5%2 3 = — : = {a3 2 ee SM} Lysine... (90) 
for any phase angle ¢ within the limits of the inequality (75). The total amount of 
light P(r, 72 ; €) reflected from the fully illuminated as well as the penumbral 
zones of the primary component is naturally represented by the sum 

P(r 1. 5 \=Ly (ty 72 5 )+Lo(M 12 5 €) (91) 
where the quantities #7, and, are (within the scheme of our approximation) given 
by equations (74) and (90); which means that the geometrical solution of our prob- 
lem which we set out to obtain in Section 4 is now complete. 

6. Discussion of the results.—The phase law (91) with its constituents as given 
by equations (74) or (77) and (90) represents the final outcome of our investigation. 
The first term on the right-hand side of (74)—of the order of (r,/R)? — is, of course, 
well known from planetary photometry. The next term, of the order of (r,/R)°, 
was first deduced by Takeda (3), and that of the order of (r,/R)* by Sen (4). The 
terms of the order of (7,7./R*)? in both (74) and (90) are new, and their sum repre- 
sents the first-order photometric effect of finite angular size of the secondary 
(illuminating) component. These terms are approximately (r,/R)* times 
smaller than the first-order reflection effect. If we consider r,/R=0-3—a realistic 
value in close binary systems—the penumbral zone, in particular, turns out to 
contribute approximately 10 per cent of the total reflected light. Russell’s 
conjecture that . . . “the penumbral effect on the reflected light may be neglected 
unless the precision of observation surpasses one part in a thousand”’ (5) proves, 
therefore, to be ill-founded ; in actual fact the penumbral zone is relatively much 
more important. 

The ensemble of these terms renders, to be sure, our phase-law exact only as 
far as quantities of fifth and higher orders in fractional radii of both components 
are negligible. The geometrical analysis of Sections 2 and 3 is, however, rigorous 
and free from any approximation; and the general formulae given there could be 
developed to any order of accuracy—at the expense only of rather troublesome 
algebraic work. If, however, such an extension were to be physically meaningful, 
it would be necessary to take account of the departures of both components from 
spherical form, due to their axial rotation and mutual tidal action. In more 
precise terms, according to the equilibrium theory of tides, the radius of a centrally 
condensed primary component rotating with the Keplerian angular velocity 
should, to the first order in small quantities, be given by 


r1(4,n) =o {x -- =(1 = ws) (“z)’ P, (cos) + ms mS (ey P,(cos¢dsin »} A 
(92) 
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where m,/m, denotes the ratio of masses of the two components, and 74 the radius 
of a sphere having the same volume as the distorted ellipsoid (8). 

The foregoing expression for r,(¢,7) should be introduced in equation (21) 
and its sequel; and the fact that it is no longer constant will affect also the limits of 
the illuminated crescent visible at any phase. The corresponding distortion of 
the secondary component will, in turn, affect not only the geometry of our problem, 
through equation (20), but also the flux of the incident radiation; for the apparent 
surface brightness of a distorted star should—apart from limb-darkening—vary 
also in proportion to local gravity (g). All these phenomena are bound to render 
a determination of the proper amount of reflected light much more involved, but 
will have to be considered in any attempt to establish the specific form of the terms 
of the order of (r/R)* or higher in the corresponding phase-law. This, however, 
constitutes a task which must necessarily be postponed for future investigation. 

In order to facilitate application of the results established in the preceding 
section to practical cases, let us expand the phase-law (91) in ascending powers 
of cose in a series of the form 


P(11,%2; «-)=L, & C,, cos "e, (93) 
n=0 


where, within the scheme of our approximation, 
1 2 (n\_1(n\ R (nr? 
com se(8) -3(8) —- sGt) * 
it tet etter 
cn EG) #3(9) + 


_ 2 (nt, 3 (nm), Fn 
c= 5e(R) +8(R) +(e 


C3 = oO, 


eee) oe i) 
c= (2) -<(%) * 


where we have abbreviated 





12 5+(7—5)u ¥ 
we ccgee 7? (95) 


Alternatively, by expanding our phase-law (91) in a Fourier cosine series we get 


k=I- 


L(r1,%2; €)= Ly X C,,' cosne, (96) 


n=0 


_ 8 fn\? | 1 (n\8 2 (n\4 2k (ryr2\2 
3e\R) + ro\R) *3e2\R) ~ PVR) T 
I/n\*  t/n\? 

3\R TAR en 
8 Ay 3s 

C2 a7m (3) * 16 (3 


ve? 
C3 =0, 


C'= 16 /(r,\2 
4 67572\R 
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The phase-law (91) or its foregoing expansions (93) or (96) represents the 
amount of light of one (the secondary) component reflected at any phase from its 
mate. ‘The corresponding amount of light of the primary component reflected 
from the secondary is, however, bound to be given by the same equations provided 
that the radii r, and r, are interchanged, the phase angle « is replaced by 7—e, 
and L, is replaced by L,. The total amount of light added by reflection to the 
intrinsic light of a close binary system is, of course, represented by the sum of 
L(r1,%23 €-)+L(re,7,; m—€); and it is this sum which should be removed from 
the observed light changes of a close eclipsing system before its light curve is 
analysed for the geometrical elements. The explicit form of this sum is so 
easily written down with the aid of the expressions already given that it need not 
be repeated here. The reader may, however, notice that the difference of 180° 
between the phase of the primary and secondary component will change the sign 
of all odd functions of « in (91), but not of even functions. As a result, the terms 
arising from reflection of light of one star from another and varying as odd 
functions of the phase, i.e. odd powers of cose in (93), or cosines of odd multiples 
of « on the right-hand side of (96), will tend to neutralize each other in the combined 
light of a close binary system, while those varying as even functions of ¢ will 
always reinforce. In the particular case of a system whose components are equal 
in size and brightness all odd terms arising from reflection should cancel, and the 
surviving even terms should closely simulate the photometric effects of ellipticity 
produced by tidal distortion. 

In concluding the present investigation, one last point deserves explicit 
notice. All expressions derived above for the amount of reflected light and its 
variation with the phase hold good irrespective of whether the relative orbit 
of the two components is circular or eccentric. If it is circular, the phase-angle 
e is simply proportional to the time and the radius-vector R is constant. Should, 
on the other hand, this orbit happen to be an ellipse, it is easy to show (g) that 
equation (3) should be rewritten as 


cose =sin(v+w) sini, (98) 


where wv denotes the true anomaly of the secondary (illuminating) component 
in its relative orbit around the primary star, and w stands for the longitude of 
periastron. Furthermore, the separation R of centres of the components will 
vary so that 


Sa = 1 2(r +ecosv), (99) 
where h denotes the semi-latus rectum of the relative orbit, and e its eccentricity. 
In consequence, the reflection effect in eccentric binary systems will no longer be 
symmetrical with respect to conjunctions; and its observed asymmetry may 
indicate to us an approximate position of the periastron.* 

Ultimately, it should be recalled that the reflection effect—like the ellipticity — 
is bound to give rise to continuous variation of light of close binary systems, 
irrespective of whether or not the system happens to be an eclipsing variable ; 
and this variation should disappear only under certain special circumstances. 
If, for instance, the relative orbit of the two stars is circular and the line of sight 
happens to be perpendicular to the orbital plane (i.e. 7=0°), then ¢=g0° all 


* For an illustrative example of such a case, the reader is referred to a study of the reflection 
effect in the system of é Aurigae (e=0-41) by K. Walter, Zs. f. Ap., 14, 62, 1937. 
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the time, and the reflected light continues to contribute a constant amount to the 
intrinsic light of the system. Should, however, the relative orbit happen to 
be an ellipse, the amount of this additive light will vary (approximately) as R-*, 
being greatest at the time of periastron passage, and smallest at apastron—a 
fact which, if observed, should indicate by inspection an approximate orientation 
of the apsidal line. 


Department of Astronomy, 
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Addendum.—After this paper had been read before the Society, T. G. Cowling 
pointed out to the writer that, in order to obtain a consistent phase-law complete 
to quantities of the order of (r,7,/R*)*, we should adjoin to (74) the effects of 
“secondary reflection ’’—i.e. the amount of the primary’s light reflected from the 
secondary component back to the primary star. This can indeed be done when 
we ensure that if L, denotes, in equation (74), the intrinsic luminosity of the 
secondary component, its actual luminosity increased by reflection at “full” 
phase (« =0) will be equal to 

2r," 
L.+L, {ote + 


In consequence, to the phase-law as represented by equation (74) should be 
added the term 


2 2 ies +sin 
Le Cate: o= aah {« €) cos € + Si ‘ho.+ 


7 


27y%2\2{I cose coste coste 
= 57a} 3 peer deer ale teh Ly (98) 
which represents the phase effect of secondary reflection. 

Within the scheme of our approximation, the foregoing equation (98) alone 
would represent the light reflected by an eclipsing system whose secondary 
component was completely dark; for the latter would still be bound to intercept 
a part of the radiation of the primary star, and to illuminate it back by reflection 
in very much the same way as lunar ashen light illuminates our Earth. Tertiary 
reflection (arising from reflection of the secondary’s light on the primary, back 
to the secondary, and returned to the primary again) would clearly constitute 


a quantity of the order of (r,?7,/R*)?L,, and therefore negligible within the scope 
of our accuracy. 
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Summary 


Relative line strengths for the 3s—3p and 3p-3d transitions in OII 
and Neti and the 4s—4p transitions in S11 and A have been calculated 
by quantum mechanics, using the intermediate coupling approximation. 
Landé g-values have been calculated and compared with experimental 
values. The parameters in the calculations have been determined empirically 
from observational spectroscopic data. Absolute strengths have been 
calculated by the Coulomb approximation. The strengths of the arrays 
3p-ns are discussed. Applications are made to the curve of growth of y Pegasi. 





1. Introduction.—The theoretical analysis of stellar spectra requires a 
knowledge of the strengths of a large number of spectral lines. In some cases 
(e.g. iron) experimental values are available for many lines, but for most of the 
lighter atoms and for most ionized atoms experimental techniques do not permit 


of the measurement of either relative or absolute line strengths. Recourse must 
then be had to theoretical calculations. General formulae have been obtained 
for a large number of possible transitions when it is permissible to assume that 
the atom or ion obeys Russell-Saunders (LS) coupling (i.e. that spin-orbit 
interaction and configuration interaction are negligible) (1). For many atoms of 
astrophysical interest the assumption of LS coupling is insufficiently accurate 
and appreciable departures from LS coupling are observed for some arrays of 
atoms in the first row of the Periodic Table. In astrophysical work such 
departures are major contributors to the large scatter of points around a curve of 
growth. Dr L. H. Aller pointed out to the writer that the 3p—3d array in O11 
appears to show substantial departures from LS coupling, and for heavier atoms 
it is well known that many strong lines in astrophysical sources have a theoretical 
strength of zero in LS coupling. 

It is evidently desirable to attempt to carry out more accurate calculations, 
including the effects of spin-orbit interaction but neglecting configuration 
interaction. Atoms for which such calculations provide a fair representation of 
the physical state are said to obey intermediate coupling (IC). IC theory is 
now standard and is described by Condon and Shortley (2). Very few applica- 
tions have been attempted, largely because of the considerable labour involved 
and the need to consider individually each array of each atom. Most of the 
work has to be done numerically, and numerical results for one array are of little 
help when studying some other array. Shortley (3) described much of the 
general theory, worked out in detail the array 3s—3p in Ne1 and compared his 
results with dispersion measurements of Ladenburg and Levy. His results 
were not wholly satisfactory, the calculated intensities were sensitive to the 
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assumed values of the parameters involved. Gottschalk (4) studied the 
3d’4s—3d’4p array in Fe1 and obtained results which represented a considerable 
advance on LS coupling calculations. His results gave much better agreement 
with experiment for many line strengths and for the Landé g-values. For a 
large number of intercombination lines he obtained finite strengths, and many 
of these agreed with experimental data. 

It is concluded from this work that while, in some cases at least, reliable 
values of the line strengths are difficult to attain, in all cases the introduction of 
IC gives a great improvement over LS calculations. (In spite of the defect of 
sensitivity mentioned above, Shortley’s IC results are much better than the 
LS values for all but one or two faint lines.) Considerable confidence is felt 
in the value of extending this type of calculation to other ions even if no experi- 
mental data are available for comparison. Indeed, the need for acquiring data 
for ions which cannot be studied experimentally is the main justification for 
the work described in this paper. 

2. Method of calculation.—Line strengths in intermediate coupling can be 
calculated from the strengths in LS coupling if we know the transformation 
matrices. These matrices are obtained by diagonalizing the complete energy 
matrix (including both electrostatic and spin interactions). Let @,, @, be the 
LS wave functions of the initial and final configurations. , and @, are repre- 
sented by column matrices. Let the transformations to intermediate coupling 
be written 


P=AP, $y =Boe,;. 
Then the matrices of the square roots of the line strengths are related by the 


(S¥?),¢=B(S"”), <A, 


equation 


where A is the transpose of A. 

If H is the energy matrix of the initial configuration in LS coupling, A is 
determined by making AHA diagonal. Similarly B is determined for the 
final configuration. The diagonal elements of AHA give the IC energies of 
the initial configuration. 

The following arrays were chosen for numerical calculations of line strengths : 


Ou 2p? 3s—2p? 3p 
2p*?P)3p-2p*(*P)3d 
Ne Il 2p* 3s—2p* 3p 
2p*(*P)3p-2p*(*P)3d 
Su 3p? 48-3" 4p 
Al 3p* 48-3p* 4p 
The first step in the calculation is the determination of the matrices of spin- 
orbit interaction. Condon and Shortley (2) describe the method in detail and 
give the matrices for the p?s configuration. We have followed their method 
through for the p*p configuration.* Throughout this paper we use a prime 
(2D’) to denote a spectral term based on a !D parent and a double prime (?S”) 
to denote a term based on a1S parent. Terms based on a °P parent are written 
without a prime. In the p? p configuration ¢,,’ refers to the outer electron, 


* The writer will gladly supply a copy of these matrices to any interested investigators. 
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One approximation made by Gottschalk (4) was to neglect the interactions 
between spectral terms based on different parents. In the s-p transitions we 
shall retain these interactions and investigate their effect on the line strengths. 
In the p-d arrays, which are more complicated, we shall neglect interparental 
interactions. Examination of the s—p results shows that interparental transitions 
are weak. ‘Their neglect for p—d is likely to be a good approximation and gives 
a valuable reduction in the computational labour involved. Johnson (§) gave a 
method which is immediately applicable if we take the LS values of the parent 
term (*P) as playing the part of the first electron in his formulae. Using Johnson’s 
method and Condon and Shortley’s determination of absolute term values (6) 
we obtain the matrix for the p?d configuration.* 

The matrices for p*s, p*p and p*d are obtained from those for p?s, p?p and 
p*d respectively by changing the sign of ¢,, and leaving the signs of ¢,,’ and C, 
unchanged. 

The complete energy matrix is obtained by adding to each diagonal term of 
the spin-orbit matrix the appropriate electrostatic energy (the centre of gravity 
of the multiplet). 

The next problem is to fix the numerical values of the spin-orbit parameters 
(¢) and the centres of gravity of the multiplets. This is, in many cases, the most 
difficult and critical step in the calculations. In the p?s and p*s configurations 
the departures from LS coupling are small. It is sufficient to work to the second- 
order terms in ¢,. We have five parameters to adjust, ¢,, and four centres of 
gravity. In principle, all the parameters could be calculated from wave functions. 
Generally, however, wave functions are not available, and an empirical approach 
must be adopted. We choose the parameters to give the best possible fit to the 
spectroscopic energy levels. The observational data are taken from the recent 
publications of C. E. Moore(7). The observed energies are equated to theoretical 
expressions written down by second-order perturbation theory. The set of 
eight equations is solved by a least squares process or a method of trial and error. 
Having determined the five parameters, we compute from them the theoretical 
energies and compare with the eight observed energies to check the adequacy of 
fit. The adopted parameters for the p?s and p‘s configurations are given in 
Table I. E(*P) denotes the centre of gravity (electrostatic energy) of the *P 
term. With the exception of Am 3p4(1D)4s?D’ the fit obtained is excellent. 
The At ?D’ doublet is predicted to have a separation of 12 cm~, but is observed 


to be 222cm™. No reasonable choice of parameters would remove this 
discrepancy. 


TABLE I 
Adopted parameters for p*s and p‘s configurations 


Ion Ou Nell Si All 
Configuration 2p” 35 2p* 3s 3p? 4s 3p* 4s 
E(*P) 185 402°0 219 452 IIO O19 134795 
E(?P) 189 o10°1 224291 113 289 138 587 
E(?D’) 206 971°0 246 391 121 515 148 689 
E(??S”) 226 849 276 664 136 8007 167 242 
Ce» 190°5 615 524 I 030 





* The writer will gladly supply a copy of these matrices to any interested investigator. 
+ Estimated from the S 111 spectrum. 
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In the p?p and p*p configurations we have two spin-orbit parameters and 
ten centres of gravity of multiplets, to be chosen to give as good a fit as possible 
for 21 energy levels. The values of the spin-orbit parameters for the outer 
electrons were rather difficult to set. An examination of the energy levels arising 
from the 'D parent indicated that the parameter ¢,,, for O11 and Ne should be 
estimated theoretically, wave functions being available from unpublished work 
of D. R. Bates for O11 and of the writer for Nett. Numerical integrations gave 
{;, =23cm™! for Ou and 46cm“ for Nell, using the standard formula 


eh? ° 1d 
Lan gett I, (~ 7 ar) PAM) Ar: 


For the purposes of calculation values of 20cm™! and 40cm™! were adopted for 
O1 and New, taking into consideration the observed energy levels. Values of 
£4, =80cm~ and 100cm~! were adopted for Sm and Al respectively. These 
latter values, based on the observations, are quite uncertain ; they are appreciably 
smaller than the core parameters (,,, and no appreciable error is likely to result 
from their adoption. The remaining parameters were estimated by a process 
of trial and error. In some cases second-order perturbation theory was used. 
This was insufficiently accurate for many of the configurations. For the latter 
a method of successive partial diagonalization was used. Using a preliminary 
set of parameters the energy matrix was set up numerically. A series of unitary 
matrices P, Q, R, ... was constructed by inspection in such a way as to reduce 
or eliminate one by one the large off-diagonal elements in the energy matrix when 


the series of transformations PHP, QP HPQ, RQPHPOR, ... Was per- 
formed. Eventually a matrix [...R QP] is found which reduces H to a new 
matrix with small off-diagonal elements. The calculation of the energies and 
the final transformation matrix is completed by second-order perturbation theory. 
This method is especially suited to our work because we have to begin with a 
matrix H computed on the basis of assumed parameters. These parameters 
must be varied to give the best possible fit of the calculated and observed energies. 
When one of the parameters is altered, we compute a new matrix H and imme- 
diately apply the transformation [...R QP] obtained previously, finishing the 
calculation by second-order perturbation theory. This method makes the trial 
of new parameters a relatively simple operation. It would be possible to apply 
a least squares technique in the final stages of the adjustment of parameters ; 
but this refinement was not felt to be justified in the present series of calculations. 

The adopted parameters for the p?p and p*p configurations are given in 
Table II. On the whole excellent fits were obtained, but considerable discre- 
pancies were encountered in the ?P and ?P’ levels in Neu, Su and An. No 
satisfactory explanation of these discrepancies has been found. 

The p*d and p‘d configurations in Ou and Neti respectively were studied 
by the above method. Only terms based on the *P parent were considered. 
The adopted parameters are listed in Table III. The spin-orbit parameter (;, 
could not be obtained from the observed energy levels. {,, was computed from 
a wave function for Neu. The value obtained was 0-3cm™; accordingly it was 
sufficient to put (,, zero forO1and Net. The observed and calculated energies 
agree remarkably well for Ou. For Net the agreement is rather poor. 
Unfortunately, the calculated energy levels are very sensitive to the adopted 
parameters. In choosing trial values for the parameters special care was taken to 
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ensure that the levels with any given J came out in the correct order of increas- 
ing energy, as judged by the observed energies. This appears to be essential if 
transition probabilities are to be computed. It is possible that further adjust- 
ments of the parameters might improve the Neti energies, but in view of the 
general sensitivity of the calculations and the labour involved, further work did 
not seem justified at the present time. 

The matrices which transform LS coupling to intermediate coupling were 
computed, using second-order perturbation theory, combined, where necessary, 
with the approximate matrices [...R QP] already obtained for many of the 
configurations. 

TasBLe II 
Adopted parameters for p* p and p* p configurations 
(units : cm-*) 
Ton Ou Nell Si All 
Configuration 2p* 3p 2p* 3p 3p" 4p 3p" 4p 

E(?S) 203 945 252 852 125 520 160 845 
E(*P) 214214 254 139 133 379 159 900 
E(?D) 211 635 251 206 130 962 158 820 
E(*S) 212 165 252956 131 047 160 990 
E(*P) 208 431 246 366 129 999 155 365 
E(*D) 206 899 249 410 128 318 157 730 
E(?F’) 228 736 274 385 140 267 170 430 
E(?D’) 229 957 277 332 140712 173 335 
E(?P’) 232 500 276 385 143 532 172 465 
E(?P”) 250 253* 305 400 156 750* 191 350° 
C 204 628 552 1075 
ig 20 40 80 100 

* Estimated. 


Tase III 
Adopted parameters in p* d and p* d configurations 
Ion Ou Ne II 
Configuration 2p*(°P)3d 2p*(?P)3d 
E(?P) 233 470 281 506 
E(?D) 234 425 280 560 
E(?F) 232 880 280 442 
E(*P) 232 554 281 040 
E(*D) 232 726 279 305 
E(*F) 231 430 280 495 
Cop 192 620 
Ca ° ° 


There is an independent check which can be applied to our work. Using the 
transformation matrices we can compute Landé g-values in intermediate coupling 
(8). In LS coupling the g-values are given by the usual formula 
J(J +1)—L(L+1)+S(S+1) 

2J (J +1) ’ 


The IC g-values for p*s and p*s configurations do not differ very greatly from the 
LS values. The IC g-values for the p*p and p*p configurations are given in 
Table IV, together with such experimental values as are available (7). The critical 
case is All 3p*4p, where there are both large deviations from LS coupling and 
adequate comparison data. Inspection of the results for A in Table IV shows 


g(SLJ)=1+ 
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that the introduction of intermediate coupling has made a great improvement in 
the Landé g-values. We find 


X (ove —&1s)* = 0°219 I 
array 


X (ors — 1c)” = 0°015 6 
array 


This provides a remarkable demonstration of the increase in accuracy gained by 
the use of intermediate coupling theory. For the d configurations experi- 
mental g-values are available for only two levels, for which the results are 


Ne 1 2p* 3d 81s Sic Sobs 
*Piie 0-667 0°707 0°70 
"Pyle 1-333 1-293 1-25 
The general improvement of the g-values and the reasonable fit obtained for 
the majority of the energy levels encourages us to proceed to the calculation of 
intermediate coupling line strengths. 


TaBLe IV 
Landé g-factors* for p* p and p* p configurations 
Ol Nell Siu All 
2p* 3p 2p* 3p 3p" 4p 3p* 4p 
2000 2000 1°958 (1°96) 2000 1°670 (1-695) 
0667 0667 0-709 (0°71) 0668 0-985 (0-983) 
1°333 1°333 1°337 (1°33) 1°332 1°179 (1°244) 
0800 o 801 0808 (0°80) o'810 1‘102 (0-918) 
1°200 1°200 1°204 (1°20) 1°207 1°243 (1°241) 
2000 1999 1'989(... ) 1-984 1970 (1987) 
2°667 2°665 2°662 (2°67) 2°658 2639 (2-638) 
1°733 1°733 1°732 (1°73) 1°735 1°724 (1°720) 
1600 1°599 1°598 (1°60) 1°591 1°593 (1°599) 
0:000 0'002 0°005 (0°00) 0'007 0032 (0:000) 
1°200 1201 1*199 (1°20) 1°203 1°179 (1-199) 
1-371 1°372 1*369 (1°37) 1°373 1°330 (1°334) 
1°429 1°429 1-428 (1°43) 1428 1°427 (1°427) 
0667 0667 0°667 (0°67) 0667 0°669 (0677) 
1°333 1°333 1°333 (1°33) 1°334 1°331 (1°332) 
0800 0800 0-801 (0°80) 0801 0805 (0804) 
1°200 1°200 1°200 (1°20) 1°197 1°201 (1°202) 
0°857 0°857 0°857 (0°86) o'861 0°859 (0°857) 
1°143 1°143 1°143 (1°14) 1°143 1°143 (1°140) 
0-667 0667 0°667 (0°67) 0667 0°666( ... ) 
1°333 1°333 1°333 (1°33) 1°333 1°333( +. ) 
* Experimental values in brackets. 


A difficulty was encountered in the A 3p* 4p calculations. An attempt was 
made to use the Landé g-values as a guide to the parameters to be adopted. This 
attempt was not successful. It was found that as certain parameters (particularly 
E(?D)) were changed in such a way as to improve the calculated energy levels the 
agreement of the Landé g-values became worse. There seemed to be no point 
in attempting any further improvement in the calculations. This difficulty 
does not affect our general conclusion that the intermediate coupling calculations 
represent a decided improvement over the LS calculations. It is, however, an 
excellent illustration of the uncertainties which are present in this type of work, 
and it is a warning against over-confidence in the numerical results of the work. 
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The final stage in the calculations is the transformation of the line strengths 
to intermediate coupling. The matrices of the square roots of line strengths 
are written down from the formulae of Condon and Shortley (9) and of Goldberg 
(10). Using the matrices previously obtained the transformation is carried out. 

3. Numerical values of relative line strengths.—The relative line strengths for 
all the arrays are collected in Tables V and VI. The numerical values are 
expressed as multiples of o, the radial integral whose value will be obtained later 
for each array. We have omitted a number of intercombination lines which are 
much too weak to appear in the Revised Multiplet Table (11). Many lines of 
moderate strength have been included even if they do not appear in the R.M.T. 


TABLE V 


Relative line strengths of s—p arrays in Ou, Neu, Sur and A 
(All expressed in units of the radial integral o*) 


Intermediate Coupling 
LS Ou Ne Su 
38-3P 38-3P 45-4p 
4°43 4°61 4°63 
0°57 °'79 0°34 
0°30 0°28 0°27 
1°94 1°67 2°69 
1°70 1°71 1°73 
1°57 1°39 1°30 
1°49 1°55 0°95 


8-00 8-00 8-00 
4°43 4°56 4°71 
1°79 1°90 1°93 
1°57 1°38 1°29 
2°05 1°95 1°96 
1°70 1°72 1°73 
0°16 O'13 oO°10 
0°30 0°28 0°27 


1°37 1°22 1°67 
0°54 
2°16 


5°93 


3°03 
0°95 





Intermediate coupling line strengths 


TABLE V (cont.) 
Intermediate Coupling 

Ou Nell si 
35-3P 38-3P 458-4p 
5°57 5°57 5°26 
3°61 3°52 3°65 
0°39 0°48 0°35 
0°43 0°43 o'74 


4°00 4°00 4°00 
2°00 2°00 2°00 


0064 
0°003 
O'O14 


‘018 


“O09 


2 
3/2 
3/2 
5/2 
3/2 


"002 
003 
‘OO! 


"004 O° 


‘Py 0 dos peer d "004 O'oll 
v2 obs — nie 0°00! 0003 


*S1/2- *Ss2 “a ‘as oct o-oo! 0°004 


* Indicates a line particularly sensitive to adopted parameters. Numerous other lines 
have a small, but finite, strength. All lines with strengths exceeding 0-02 have been included. 


There are certain sum-rules valid in intermediate coupling. These provide 
useful numerical checks on the computations but do not provide any check on the 
physical validity of the transformation matrices. 

On examining the results of the computations for the s—p arrays it is seen that 
intercombination lines are fairly weak. They are particularly weak if the initial 
and final terms are based on different parents. It follows that it is a good approxi- 
mation to neglect interactions between terms based on different parents. 
A second conclusion which may be drawn from Tables V and V1 is that in general 
weak lines are very much more sensitive to departures from LS coupling than are 
strong lines. ‘This conclusion was anticipated: the effect of IC is to redistribute 
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TaBLe VI 


Relative line strengths for O 11 2p*(P)3p—2p?(*P)3d and for 
Ne 11 2p'(*P)3p-2p*(°P)3d 
(All expressed in units of the radial integral o”) 


LS Intermediate Estimated 


Coupling intensities 
Ol Ne Il O1 Nell 
13°33 13°24 g°O1 4 
6-67 “51 7°99 5 


Coupling 


48:00 "90 31°90 
2°40 79 1°30* 
33°60 ‘45 4°88* 


8-40 3°81 7°99 
5°40 “19 4°84 
0-60 “48 O°17 
0-60 ‘47 4°22* 


8°33 . 5°39 
3°33 , I-71 
1°67 . 1°81 
1°67 ‘ 8-24 


27°00 : 14°75 
1500 7°71 
3°00 3°86 


0-60 082 
0°33 , 018 
0-066 0°044 


60-00 60-00 
28-99 
7°08 
11°99 
2°77 
1°76* 
3°83* 
0-20* 
o-o11* 


oOfwueut DAO 


1°04 
1'07* 
0°42 
o-o10* 
o89* 
oll 
0°33 
3°35* 


NN HRS 


12°87 
6-20 
2‘91 
1°34 
1°51 
1°63 


WNHNM RF WWWwHsl 





Intermediate coupling line strengths 
TABLE VI (cont.) 


. Intermediate Estimated 
LS . ; — 

Coupling intensities 
Ou Net Nell 
10°50 3°85 6°33 
1°33 0°34 0°47 
0°83 5°26 1°30 
4°50 1°79 2°19 
4°17 0°84 3°94 
4°50 13°29 5°96 
4°17 9°54 3°67 


Coupling 


2) 
AQuVNHtO | 
— 


i 


20°00 16°57 
13°33 10'II 
6°67 4°66 


pu 


36-00 33°80 

18-90 8-49 

3/2 75° 1°93 

5/2 8-10 13°62 

3/2 9°60 10°79 

1/2 7°5° 3°14 

3/2 0°90 3°59 

32" = 8 4°69 


SOnuUnwtrh uw 


2D, 2D, 2 


3/2 5/2 
5/2 5/2 


nu 


4G, 4 

Ss2~*D5)2 
32 8638/2 
32- =/2 


*P5)2- °F, 2 


3/2 5/2 


5/2 5/2 


2 
*Dz2- Fyre 
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TaB_e VI (cont.) 
Intermediate Estimated 
Coupling intensities 
Ol Nell Ou Nell 
28 mF I 
43 as 4 
12 es aE: 


‘22 


LS 
Coupling 


“49 


oom o'4I 
5/27 5) a «ne °0°78 
5/2 ose 0°87 
es wad 3 or'll 
aa-/8 aed ins 0°08 
yes sis ins 2'12 


*Pija- "Pie eee 0-16 
ya «3/2 oe . o'16 


*P3jo- *P 5/2 eee eee 3°78 
ye- a8 “s ait 2°27 
3/2 = 3/2 see ‘ee 0°95 


*Dije- *P 3/2 eee tee o'12 
a Jor aia 006 
sc 3/2 see eos 0-09 


4P, "Ds /2 oa ves o’51 
sa 8 oes aes 0°23 
3 95/2 eee eee 0°09 
se ase nee a 0°16 
a on oes ove o'12 


* Denotes a line particularly sensitive to adopted parameters. This table includes all 
lines mentioned in the Revised Multiplet Table and all others with strengths not less 
than 0-06. 
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energy among various possible transitions and a weak line will be subject to a 
greater proportional change. In the p-d arrays, however, there are a number of 
strong lines with strengths substantially different from the LS values. The 
multiplet 4P-*P of O11 is of especial interest in this connection. Furthermore 
there are a number of very strong intercombination lines in the p—d arrays of both 
Om and Net. 

There are, unfortunately, no tables of accurate laboratory intensities for any 
of the transitions which we have studied. ‘The available eye-estimates (11) have 
been included in Table VI. In some multiplets (e.g. ?7P—*P in Ou and Ne uand 
4P-*P in O11 the IC strengths give much better agreement with the estimated 


T T T T T ' | q ' 





log S 
2 








L 1 l 
5 6 7 8 
Laboratory Intensity 


Fic. 1.—Relation between the logarithm of the LS coupling strength and the laboratory 
intensity estimates for Ne II 3p-3d. 
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Fic. 2.—Relation between the logarithm of the IC strength and the laboratory intensity 
estimates for NeII 3p-3d. Filled circles denote lines allowed in LS coupling—the same lines 
as in Fig. 1. Crosses denote intercombination lines. 


intensities than that given by LS strengths. In others (e.g. ‘D-*F in Net) 
the agreement is disappointing. The most encouraging feature, however, is the 
behaviour of the intercombination lines in Neu. The theoretical strengths 
obtained for the *D-*F multiplet (amongst others) in Net are excellent. To 
illustrate this Figs. 1 and 2 show the relation between the logarithms of the 
calculated relative strengths and the estimated intensities. In Fig. 2, inter- 
combination lines have been distinguished from lines allowed in LS coupling. 
The agreement of the general relationships between strength and intensity for 
the two groups of lines is very satisfactory and gives added confidence in the 


9 
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results of the calculations. There is a large residual scatter, due to a combination 
of the uncertainties in the calculations and the unknown errors in the intensity 
estimates. 

4. The transition arrays 3p-ns.—The spectra of some early-type stars (for 
example, the B2 supergiant 55 Cygni) contain not only lines from the array 
38-3p of O11 but also some lines of O11 3p—4s. It is of interest to consider the 
line strengths of the latter transition in intermediate coupling. The matrix of the 
square roots of the line strengths, when expressed in the appropriate o as unit, is 
the same for all the arrays 3s—3p, 3p-ns (n>4). ‘This matrix has the special 
property that for any given level of the 3p configuration there is only at most one 
level in ns with any given J to which a transition is permitted in LS coupling. 
This is an exceptional property due to the simplicity of the ps configuration. 
Assume that departures from LS coupling are small for both the 3p and ns 
configurations and consider the order of magnitude of the various contributions 
which arise in the matrix of square roots of line strengths when transformed to 
intermediate coupling. It can be shown that transitions which are forbidden in 
LS coupling are also forbidden in intermediate coupling if contributions of the 
second order of small quantities are neglected. It can also be shown that to the 
same approximation, the arrays 3s—3p, 3p-ns (n>4) have the same strengths 
(expressed in units of the appropriate o”). Thus the calculation of the absolute 
strengths of lines in the arrays 3p-ns (n >4) merely requires an estimate of o?, 
together with the use of the relative strengths obtained above for the 3s—3p array 
and a verification that the departures from LS coupling are small. 

5. Absolute strengths.—The calculations so far discussed have been concerned 
with the relative strengths of lines expressed in multiples of o?, where 


ie as ({-rPeuyper, }.. 1)dr), 


and P(nl), P(n',/—1) are the radial wave functions of the series electron in its 
initial and final states. ‘This formula is an approximation based on the neglect of 
configuration interaction. It would be difficult to obtain a more accurate formula. 
The approximation probably is adequate for our purposes, particularly as there is 
little cancellation of the wave functions in the cases with which we are concerned. 

Two methods are available for the determination of o?._ One method is that 
described by Bates and Damgaard (12), usually known as the Coulomb approxi- 
mation. ‘The second method requires the numerical determination of self- 
consistent field wave functions, with or without exchange and polarizability 
corrections. Use of the Coulomb approximation provides a simple and, as far 
as is known at the present time, a fairly reliable method of calculating o2. We 
shall confine ourselves to this method. The method has been applied to the 
transitions in Ou, Su and At and the results, together with the author’s (13) 
earlier results for Net, are given in Table VII. 


TaB.e VII 

o 
Array (atomic units) 
38-3p 4°2 
3p-3d I'I7 
35-3P 3°3 
3p-3d 0°99 
45-4) 6°5 
45-4 52 
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The relative strengths of the lines given in Tables V and VI can be converted 
to absolute strengths by multiplication of the relative strength by the appropriate 
o*, Transition probabilities and f-values can then be calculated, if required, 
using standard formulae. 

6. Curves of growth for y Pegasi.—One of the principal applications of the 
strengths of large numbers of spectral lines is the construction of curves of growth 
for individual stars. Among the numerous studies of this type the work of Aller 
(14) on y Pegasi is of special interest in connection with our computations on O11. 
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Fic. 3.—Curve of growth for OII 3p—3d in y Pegasi, using LS coupling strengths. 
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Fic. 4.—Curve of growth for OIL 3p—3d in y Pegasi, using IC strengths. Filled circles 
denote lines allowed in LS coupling—the same lines as in Fig. 3. Crosses denote intercombination 
lines. 





Aller measured the equivalent widths of a large number of lines and constructed 
curves of growth for various ions. Examination of his data for O11 shows that 
while the 3s—3p array gives a curve of growth with the observational points 
scattered reasonably closely about the theoretical curve, the 3p—3d array obser- 
vations show a disturbingly large scatter. Part of this scatter is presumably due 
to observational errors and part may be due to errors in the assumed f-values. 
The latter source of scatter can now be investigated using our new f-values. We 
have re-plotted the curves of growth for the 3p~3d array in Om using LS and IC 
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line strengths and the results are shown in Figs. 3and4. The notation and method 
of calculation are the same as used by Aller. It will be seen that there is an 
appreciable reduction of scatter when the IC strengths are used. Unfortunately 
the remaining scatter is rather large and a theoretical curve cannot be well fitted to 
the observations. The writer understands that Aller is obtaining new intensity 
measurements with greater accuracy and further detailed discussion of the bulk 
of the stellar data at the present time does not seem justified. 

The most important changes introduced in IC are in the *P—*P multiplet. 
Data for this multiplet are collected in Table VIII (o? = 1-17 has been used). 


Tas.e VIII 
Multiplet *P—*P of O 11 3p-3d 


Designation A log S Aller 


’ R.M.T. yPeg 
: ae (A) LS IC Intensity E.W. (mA) 


5/2 sB/2 4 169°23 1°09 0°66 blend 
3/2 3/2 4 140°74 o'19 —o'40 és 
1/2 1/2 4121°48 —o-o!l 0-79 28 
5282 4 156°54 0°72 0°32 16 
32/8 4 129°34 0°69 —o-'ol 14 
3/2 5/2 4 153°30 o'72 I‘1g 44 
1/2 3/2 4132°81 0:69 1-05 40 


The improvement in the strengths is very considerable, and gives increased 
confidence in our results. 

It should also be noted that for three intercombination lines observed by Aller 
the theory predicts finite (and not particularly small) intensities. These are 
plotted as crosses on the curve of growth. Considering the uncertainties of the 
calculations these results are quite satisfactory. 

7. Conclusion.—It has been shown that a considerable improvement in the 
strengths of permitted lines can be attained if we perform calculations at the level of 
approximation implied by intermediate coupling theory. The improvement 
obtained would justify much further work, and a wide field is open for routine 
calculations. Investigations of the effect of configuration interaction on the line 
strengths are much to be desired, and additional experimental work must be 
performed to assess the reliability of many of the theoretical calculations. 
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MEAN AREAS AND HELIOGRAPHIC LATITUDES OF 
SUNSPOTS IN THE YEAR 1950 


Royal Greenwich Observatory 
(Communicated by the Astronomer Royal) 


(Received 1954 February 5) 


The following results are in continuation of those given in M.N., 113, 786, 
1953, and are derived from the measurement at Herstmonceux of photographs 
taken at the Royal Observatories of Greenwich and the Cape and at the Kodaikanal 
Observatory, India. The Greenwich photographs were all taken at the new 
site at Herstmonceux Castle. 


Table I gives the mean daily areas of umbrae, whole spots and faculae for 
each synodic rotation of the Sun included in the year 1950; the means for the 
year as a whole are included in Table II, which summarizes the yearly values 
since the last sunspot minimum of 1944. 


TABLE I 
Mean daily areas 
Projected* Corrected for foreshorteningt 
No. Rotation Days — A _ = A” — 
of commenced photo- ao Whole Seite Minden Whol 
rotation 5 i graphed spots spots 
1949-50 
1288 Dec. 20°18 28 266 1691 1838 192 1235 2269 
1289 Jan. 16°52 27 274 1559 1810 206 1186 2220 
1290 6$Feb. 12°86 27 564 3429 2013 422 2607 2541 
1291 Mar. 12°19 28 389 2215 1752 300 1728 2232 
1292 Apr. 8-49 27 514 3300 1781 365 2334 2227 
1293 May 5°74 27 328 1922 1742 250 1480 2150 
1294 June 1-96 27 246 1441 1600 182 1091 1966 
1295 June 29°15 27 264 1543 1467 198 1165 1806 
1296 July 26-36 28 294 1774 1270 212 1260 1600 
1297 Aug. 22°59 27 133 743 ««1184 103 579 =: 1437 
1298 Sept. 18-84 27 144 780 1127 103 567 1361 
1299 6Oct. 16°13 27 148 goo 886 III 685 1038 
1300 6Nov. 12°43 28 182 1112 882 129 806 1060 
* Expressed in millionths of the Sun’s disk. 
+ Expressed in millionths of the Sun’s hemisphere. 








Faculae 


Table III gives for each rotation in the year 1950 the mean daily area of 
the whole spots (corrected for foreshortening) and the mean _ heliographic 
latitude of the spotted areas for both northern and southern hemispheres. 

The mean heliographic latitude of the entire spotted area and the mean 
distance from the equator of all spots are also tabulated. The mean values 
for the year 1950 are included in Table IV. 

Tables II and IV are in continuation of similar tables in Monthly Notices : 
for the years 1874 to 1888, 49, 381, 1889: 1889 to 1902, 63, 465, 1903: 
IQOI to 1914, 76, 402, 1916: 19713 to 1924, 85, 1007, 1925: 1923 to 1933, 
94, 870, 1934: and 1933 to 1945, 110, 501, 1950. 

The rotations in Tables I and III are numbered in continuation of 
Carrington’s series. 
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Tas_e II 


Mean daily areas 


A 


No. of days Projected* Corrected for foreshorteningt 





‘ 


os A. 
c ‘ ¢ + 


Photo- Without Umbrae Whole Faculae Umbrae Whole Faculae 
graphed spots spots spots 
366 157 30 160 284 23 126 344 
365 14 102 560 774 78 429 940 
365 ° 389 2381 1794 291 1817 2188 
365 ° 558 3559 2326 405 2637 2894 
366 ° 419 2618 1849 314 1977 2331 
365 ° 485 2873 2140 356 2129 2597 
365 3 274 1638 1421 203 1222 1754 
* Expressed in millionths of the Sun’s disk. 
+ Expressed in millionths of the Sun’s hemisphere. 





Tas_e III 


Spots north Spots south 
of the equator of the equator Mean 
- “A ~ ¢ A 5 Mean distance 
Mean Mean latitude from 
° ‘ Mean ‘ Mean : : 
No. Rotation dail helio- dail helio- ofentire equator 
aily a aily : 
of commenced graphic graphic spotted of all 
rotation U.T. area latitude area latitude area spots 
1949-50 ; c 2 ; 
1288 Dec. 20°18 1022 13°87 213 6:06 + 10°44 12°53 
1289 Jan. 16°52 960 14°55 226 11°42 + 960 13°95 
1290 Feb. 12°86 1820 12°20 787 14°27 + 4°21 12°82 
1291 Mar. 12°19 1083 18-07 646 12°77 - 6°55 16°09 
1292 Apr. 8-49 1362 13°56 972 14°07 + 2°06 ¥a-py 
1293 May 5°74 1164 II‘O1 315 12°59 5°98 11°35 
1294 June 1°96 521 9°79 570 17°98 4°72 14°07 
1295 June 29°15 808 13°42 357 13°82 + 5:08 13°54 
1296 July 26-36 433 12°46 1121 3°06 11°64 
1297 Aug. 22°59 272 15°32 10°23 + 1°78 12°62 
1298 Sept. 18°84 354 15°45 10°28 + §°81 13°51 
1299 Oct. 16°13 295 13°47 10°75 o'31 11°93 
1300 Nov. 12°43 15°04 12°59 3°70 14°04 








TABLE IV 


Spots north Spots south 
of the equator of the equator Mean 
, ~ s f A ~ Mean distance 
M Mean Mean latitude from 
ean ‘ Mean : : 
; helio- ; helio- ofentire equator 
daily ; daily : 
Days pie graphic asia graphic spotted of all 
photographed latitude latitude area spots 
366 42 19‘00 83 22°81 —8-70 21°53 
Old Cycle 7 4°17 6 7°61 —1'l4 5°72 
366 N 
ew cycle 35 22°18 77 24°02 —9°63 23°45 
365 20°13 309 20°26 —8-92 20°22 
365 20°74 690 += 18-79 +5°73 20°00 
365 16°58 1645 17°86 —4'91 17°38 
366 936 14°92 1041 13°53 —o'06 14°19 
365 14°30 O51 =: 12°14 +2°49 13°33 
365 771 =: 1362 452 = 13°05 +3°77 13°41 
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The chief features of the record for the year 1950 are as follows :— 

(1) The mean daily area of sunspots, corrected for foreshortening (Table II), 
although still large, shows a marked decrease as compared with the previous 
year amounting to 42-6 per cent. The decrease from the maximum year (1947) 
amounts to 53-6 per cent. 

The mean area of faculae for 1950 likewise shows a considerable decrease 
though not so marked as for sunspots. As compared with 1949 the drop in 
mean daily area is 32-5 per cent and as compared with the maximum year 
39°4 per cent. 

The ratio of mean corrected areas of faculae/sunspots is I-44. 

A very large and active group with a mean area of 1853 millionths of the 
Sun’s hemisphere and central meridian passage on February 20 contributed 
to the high area value of 2607 millionths for Rotation No. 1290 (Table I). 

Two other large spots with respective mean areas of 1527 and I0II were 
present during Rotation No. 1292. Thereafter, the mean areas for the synodic 
rotations were appreciably less, especially after Rotation 1297. 

(2) On three days there were no sunspots present, viz. December 21, 22 
and 23, and no faculae were seen on October 22. 

(3) The mean weighted latitude of all spots 13°-4 (Table IV) was almost 
the same as in the previous year, thus representing what is, no doubt, a 
temporary halt in the latitude progression towards the equator. This irregularity 
is shown particularly by the southern spots. The latitudes of the three largest 
spot groups were + 10°-I, +13°-4 and —12°-o. 

Table IV also shows that the ratio of areas of northern spots/southern spots 
for 1950 is I-71. 

(4) The number and distribution, northern and southern hemispheres, of 
spot groups of (a) two days’ duration or longer; (5) one day only; are as follows : 


(a) (6) 

Northern spots 189 33 
Southern spots 126 27 
Total 315 60 


The range in latitude for the longer-lived spots was from 0° to 40°. The latitudes 
of the 1-day spots ranged from 1° to 35°. 

A short descriptive survey of the sunspots of 1950 has been given in Monthly 
Notices, 111, 233, 1951. 

An appended table gives the mean daily areas of sunspots and of faculae 
(corrected for foreshortening and expressed in millionths of the Sun’s hemisphere) 
for each calendar month of 1950. 


Monthly mean daily areas of sunspots and faculae 
1950 
Month Spots Faculae Month Spots Faculae 
Jan. 1416 2278 July 1256 1837 
Feb. 1927 2262 Aug. 1069 1496 
Mar. 1660 2330 Sept. 488 1350 
Apr. 2359 2155 Oct. 571 1220 
May 1578 2230 Nov. 745 842 
June 1048 1931 Dec. 610 1145 


Royal Greenwich Observatory, 
Herstmonceux Castle, Sussex : 
1954 January 29. 





CONTENTS 


Meeting of 1954 January 8 : 
Fellows elected 
Junior Member elected 
Presents announced : 
Award of the Gold Medal to Dr wW. Baade 


I. W. Busbridge and D. W. N. Stibbs, On the intensities of interlocked a 
lines in the Milne-Eddington Model ; ; 


A. B. Hart, Motions in the Sun at the photospheric level. IV. The equatorial 
rotation and possible velocity fields in the photosphere 


T. R. Kaiser, Theory of the meteor height distribution obtained from radio-echo 
observations. I. Shower meteors 


T. R. Kaiser, Theory of the meteor height distribution obtained from radio-echo 
observations. II. Sporadic meteors 


S. Evans, Scale heights and pressures in the upper atmosphere from radio-echo 
observations of meteors 


R. D. Davies, An analysis of bursts of solar radio emission and their association 
with solar and terrestrial phenomena 


A. D. Thackeray, The helium star HD 168476 


Zdenék Kopal, Photometric effects of reflection in close binary systems ... 


R. H. Garstang, Intermediate coupling line strengths 


Royal Greenwich Observatory, Mean areas and heliographic latitudes of 
sunspots in the year 1950 ... 


Printed in London by Taylor & Francis, Ltd., 1954 October 





